TECHNICAL  REPORT 


Defense  Threat  Reduction  Agency 
8725  John  J.  Kingman  Road,  MS 
6201  Fort  Belvoir,  VA  22060-6201 


DTRA-TR-15-80 


Irreversible  Phase- Changes  in 
Nanophase  RE-doped  M2O3  and 
their  Optical  Signatures 

Distribution  Statement  A.  Approved  for  public  release;  distribution  is 
unlimited. 


December  2015 


HDTRA1 -10-1 -0005 
Hergen  Eilers 
Prepared  by: 

Washington  State  University 
1 1 25  NE  Washington  St. 
Pullman,  WA  99164 


DESTRUCTION  NOTICE: 


Destroy  this  report  when  it  is  no  longer  needed. 

Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  THREAT  REDUCTION 
AGENCY,  ATTN:  DTRIAC/  J9STT,  8725  JOHN  J.  KINGMAN  ROAD, 
MS-6201 ,  FT  BELVOIR,  VA  22060-6201 ,  IF  YOUR  ADDRESS 
IS  INCORRECT,  IF  YOU  WISH  IT  DELETED  FROM  THE 
DISTRIBUTION  LIST,  OR  IF  THE  ADDRESSEE  IS  NO 
LONGER  EMPLOYED  BY  YOUR  ORGANIZATION. 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Allison  Lyalikov 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassfied 

SAR 

121 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

703-767-4663 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


UNIT  CONVERSION  TABLE 

U.S.  customary  units  to  and  from  international  units  of  measurement 


U.S.  Customary  Units 


Multiply  by 


Divide  byf 


International  Units 


Length/Area/Volume 

inch  (in) 
foot  (ft) 
yard  (yd) 

mile  (mi,  international) 
mile  (nmi,  nautical,  U.S.) 
barn  (b) 

gallon  (gal,  U.S.  liquid) 
cubic  foot  (ft3) 

Mass/Density 
pound  (lb) 

unified  atomic  mass  unit  (amu) 
pound-mass  per  cubic  foot  (lb  ft  3) 
pound-force  (lbf  avoirdupois) 
Energy/W  ork/Power 
electron  volt  (eV) 
erg 

kiloton  (kt)  (TNT  equivalent) 
British  thermal  unit  (Btu) 

(thermochemical) 
foot-pound-force  (ft  lbf) 
calorie  (cal)  (thermochemical) 
Pressure 
atmosphere  (atm) 
pound  force  per  square  inch  (psi) 
Temperature 
degree  Fahrenheit  (°F) 
degree  Fahrenheit  (°F) 

Radiation 

curie  (Ci)  [activity  of  radionuclides] 
roentgen  (R)  [air  exposure] 
rad  [absorbed  dose] 
rem  [equivalent  and  effective  dose] 


2.54 

x  10-2 

meter  (m) 

3.048 

x  10  1 

meter  (m) 

9.144 

x  10_1 

meter  (m) 

1.609  344 

x  103 

meter  (m) 

1.852 

x  103 

meter  (m) 

1 

00 

b 

X 

square  meter  (m2) 

3.785  412 

x  10  3 

cubic  meter  (m3) 

2.831  685 

x  10-2 

cubic  meter  (m3) 

4.535  924 

x  10  1 

kilogram  (kg) 

1.660  539 

o 

r- — 1 

X 

kilogram  (kg) 

1.601  846 

x  101 

kilogram  per  cubic  meter  (kg  m“3) 

4.448  222 

newton  (N) 

1.602  177 

x  10  19 

joule  (J) 

1 

x  10-7 

joule  (J) 

4.184 

x  1012 

joule  (J) 

1.054  350 

x  103 

joule  (J) 

1.355  818 

joule  (J) 

4.184 

joule  (J) 

1.013  250 

x  105 

pascal  (Pa) 

6.984  757 

x  103 

pascal  (Pa) 

[T(°F)-32]/1.8 

degree  Celsius  (°C) 

[T(°F)  + 459.67]/!. 8 

kelvin  (K) 

3.7 

x  1010 

per  second  (s_1)  [becquerel  (Bq)] 

2.579  760 

x  10^ 

coulomb  per  kilogram  (C  kg1) 

1 

x  10“2 

joule  per  kilogram  (J  kg1)  [gray  (Gy)] 

1 

x  10  2 

joule  per  kilogram  (J  kg-1)  [sievert  (Sv)] 

of  SI  units  may  be  viewed  at  http://www.bipm.org/en/si/. 

^Multiply  the  U.S.  customary  unit  by  the  factor  to  get  the  international  unit.  Divide  the  international  unit  by  the  factor  to  get  the 
U.S.  customary  unit. 


2015-11-16 


ABSTRACT 


Luminescent  rare-earth  doped  metal  oxides  are  evaluated  as  potential  temperature 
sensors  that  are  seeded  into  explosive  fireballs  and  collected  post-detonation  for  analysis. 
The  sensors  consist  of  amorphous  precursors  which,  once  subjected  to  thermal  exposure, 
undergo  irreversible  phase  transitions  such  as  decomposition,  nucleation,  grain  growth, 
and  phase  transformations.  The  extent  and  degree  of  these  phase  transitions  depends  on 
temperature  and  heating  duration.  The  phase  transitions  are  monitored  via  spectral 
changes  in  the  fluorescence  of  the  luminescent  rare-earth  dopants.  The  spectral  emission 
of  these  dopants  is  very  sensitive  to  the  configuration  of  the  oxygen  ions  surrounding  the 
dopant.  As  this  configuration  changes,  the  spectral  emission  properties  change  and  this 
change  can  be  correlated  with  the  temperature  the  materials  were  exposed  to.  The 
correlation  is  based  upon  laboratory-based  reference  measurements.  The  functionality  of 
these  temperature  sensors  was  demonstrated  during  several  explosion  tests.  By  using  two 
different  sensor  materials  combined  with  a  kinetic  analysis,  it  is  possible  to 
simultaneously  extract  information  about  temperature  and  heating  duration.  Initial 
results  using  this  approach  were  demonstrated.  Furthermore,  luminescent  tracking 
particles  were  developed  to  determine  the  origin  of  the  temperature  sensors  during  an 
explosion.  These  tracking  particles  consist  of  fully  crystallized  rare-earth  doped  metal 
oxides.  Using  different  dopants  with  different  spectral  signatures,  provides  a  means  to  tag 
the  temperature  sensors  placed  within  different  locations  of  the  explosive  fireball  and  to 
track  their  origin. 
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I.  SUMMARY 

The  Defense  Threat  Reduction  Agency  (DTRA)  is  developing  the  science  and  technologies 
to  defeat  weapons  of  mass  destruction  such  as  biological  agents.  Developing  such 
weapons  requires  knowledge  of  the  temperature  and  heating  duration  inside  the  fireball 
of  an  explosion.  Current  approaches  rely  on  the  use  of  thermocouples  (TCs)  for  such 
measurements.  However,  TCs  are  limited  in  their  response  time  and  tend  to  get  destroyed 
during  the  explosion.  As  such,  DTRA  is  supporting  the  development  of  alternative  sensors 
that  are  capable  of  measuring  and  recording  temperature  and  heating  duration  under 
extreme  conditions  to  ensure  that  the  temperatures  reached  in  a  post-detonation  fireball 
is  sufficiently  high  to  destroy  biological  agents.  In  order  to  be  relevant  for  this  specific 
application,  DTRA  expects  these  sensors  to  be  of  about  the  same  size  as  biological  agents. 
Furthermore,  DTRA  needs  the  means  to  map  the  spatial  temperature  profile  within  the 
fireball  (i.e.  the  movement  of  the  temperature  sensors  needs  to  be  tracked). 

Our  approach  is  based  upon  monitoring  irreversible  phase  changes  (e.g.  decomposition, 
nucleation,  grain  growth,  and  phase  transformations)  that  occur  during  heating  of 
precursor  of  metal  oxides.  These  phase  changes  are  probed  by  rare-earth  (RE)  dopants 
via  changes  in  their  optical  spectra  (i.e.  spectral  changes  due  to  heating-induced 
rearrangements  of  the  local  environment  of  the  RE  dopants).  Laboratory-based  reference 
measurements  in  which  the  sensor  material  is  heated  for  specific  times  to  various 
temperatures  are  used  for  calibration.  Changes  in  the  spectral  peak  width,  peak  position, 
relative  peak  intensities,  etc.  can  then  be  used  to  determine  the  temperature.  By  using  two 
different  sensor  materials  and  performing  a  kinetic  analysis,  it  is  possible  to  extract 
temperature  and  heating  duration.  In  addition,  fully  crystalline  particles  can  be  added  to 
the  temperature  sensors  as  references  and  for  tracking  purposes. 

Modeling  of  the  heating  process  and  response  of  the  temperature  sensors  has  highlighted 
the  importance  of  a  fast  laboratory-based  heating  technique.  If  the  effective  ramp  time  is 
too  long  compared  to  the  actual  isothermal  heating  duration,  significant  changes  will 
occur  in  the  sensor  material  that  could  lead  to  misinterpretation  of  the  results. 

Our  initial  laboratory-based  reference  measurements  were  performed  using  pyroprobe 
heating  combined  with  TC-based  temperature  measurements.  However,  heating  of  the 
temperature  sensors  was  limited  by  heat  conductivity  and  it  took  several  seconds  for  the 
temperature  sensors  to  reach  their  final  temperature  even  though  the  heating  coil  reached 
its  final  temperature  within  a  few  milliseconds.  Next,  we  used  indirect  heating  via  a  CO2 
laser.  The  temperature  sensors  were  embedded  between  graphite  sheets  and  the  sheets 
were  then  illuminated  by  the  laser.  This  approach  led  to  much  faster  heating  times,  but 
still  not  fast  enough  to  compare  with  conditions  in  an  explosive  fireball.  Finally,  we 
changed  to  direct  laser  heating  using  a  CO2  laser  and  monitoring  the  temperature  via  3- 
color  pyrometry.  This  approached  allowed  us  to  perform  our  reference  measurements 
using  sufficiently  fast  and  high-enough  temperature  rises. 

Different  types  of  our  temperature  sensors  were  tested  during  actual  explosion  tests  at 
the  Naval  Surface  Warfare  Center,  Indian  Head  Explosive  Ordnance  Disposal  Technology 
Division  (NSWC  IHEODTD)  under  various  conditions.  These  tests  were  performed  in 
open  air  and  closed-chamber  settings  using  non-aluminized  and  aluminized  charges. 
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After  each  explosion,  the  debris  including  our  temperature  sensors  was  collected  and 
returned  to  our  lab  for  analysis. 

Early  results  showed  that  precursor  Eu-doped  yttria  (p-Eu:Y203)  and  precursor  Eu-doped 
zirconia  (p-Eu:Zr02)  are  suitable  temperature  sensor  materials.  Because  of  its 
hypersensitive  transition,  europium  is  uniquely  useful  for  monitoring  temperature- 
induced  structural  changes.  Both  sensor  materials  were  successfully  tested  during  actual 
explosion  tests.  The  optical  analysis  of  these  materials  was  based  upon  changes  in  the  full- 
width  half-maximum  (FWHM)  for  Eu:Y203  and  on  changes  in  fluorescence  and 
fluorescence-excitation  peak  position  and  relative  peak  intensities  for  Eu:Zr02. 

The  more  recent  analysis  of  sensor  materials  has  focused  on  using  changes  in  relative 
peak  intensities  instead  of  changes  to  FWHM  and  shifts  in  peak  position.  Combining 
temperature  sensor  materials  with  tracking  sensors,  with  crystalline  (c-)  reference 
materials,  or  mixing  two  different  sensor  material  for  combined  temperature  and  heating 
duration  analysis  can  lead  to  spectral  interferences  from  the  different  emission  sources, 
and  we  found  that  using  the  relative  fluorescence  peak  intensities  as  indicators  is  less 
sensitive  to  spectral  interferences.  Also,  fluorescence  measurements  are  less  time 
consuming  and  require  simpler  experimental  setups  than  fluorescence-excitation 
measurements  -  an  important  consideration  for  the  development  of  potential  portable 
analysis  equipment. 

Time-dependent  measurements  using  p-Eu:Zr02  showed  that  the  material  response  to  a 
heating  event  is  very  fast  (millisecond  timescale  and  possibly  faster)  and,  dependent  on 
the  temperature,  the  changes  in  the  material  can  be  completed  within  a  few  hundred 
milliseconds.  A  mixture  of  p-Eu:Zr02,  c-Ho:Zr02,  and  p-Dy:Y2o3  combined  with  a  kinetic 
analysis  approach  is  currently  being  used  to  extract  temperature  and  heating  time 
information.  We  have  used  this  mixture  in  combination  with  direct  laser  heating  to 
demonstrate  that  we  can  measure  heating  durations  as  short  as  100  ms. 

The  following  materials  were  evaluated  as  potential  tracking  sensors:  Eu:Y203,  Tb:Y203, 
Dy:Y203,  Sm:Y203,  Er,Yb:Y203,  Ho:Y203,  Pr:Y203,  Nd:Y203,  andTm:Y203,  with  several  of 
them  tested  during  explosion  tests.  We  found  that  the  critical  selection  criteria  for 
tracking  particles  is  the  spectral  overlap  between  the  tracking  materials  and  the 
temperature  sensors,  as  the  tracking  particles  are  found  to  have  a  much  higher 
fluorescence  intensity  than  the  temperature  sensors.  Based  upon  a  careful  evaluation, 
Dy:Y203,  Ho:Y203,  Pr:Y203,  Nd:Y203,  and  Tm:Y203  appear  to  be  the  most  promising 
tracking  sensor  materials  in  combination  with  temperature  sensors  such  as  Eu:Y203  and 
Eu:Zr02.  However,  when  combined  with  the  p-Eu:Zr02,  c-Ho:Zr02,  and  p-Dy:Y2o3 
mixture,  Dy:Y203  and  Ho:Y203  cannot  be  used. 

Temperature-induced  phase  changes  in  RE-doped  precursor  metal-oxides  have  been 
successfully  demonstrated  as  temperature  sensors  in  explosive  fireballs.  The  temperature 
readings  from  these  sensors  are  always  higher  than  the  readings  from  TCs,  which  is 
expected  due  to  their  faster  response  time.  The  response  time  of  these  sensor  material  is 
at  least  as  fast  as  a  few  milliseconds  and  possibly  faster.  Extracting  temperature  and 
heating  duration  requires  the  simultaneous  use  of  two  sensor  materials  combined  with  a 
kinetic  analysis.  The  temperature  sensors  can  be  combined  with  tracking  particles  to 
identify  the  location  of  their  origin. 
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II.  INTRODUCTION 

According  to  DTRA,  its  "...  counter-WMD  mission  includes  research  and  development  of 
weapons,  in  particular  Agent  Defeat  weapons.  Such  weapons  utilize  a  combination  of 
explosive,  thermal  and  chemical  kill  mechanisms,  to  create  high  temperatures  of  critical 
duration  to  neutralize  biological  weapon  agents.  When  such  a  weapon  is  deployed,  the 
resultant  turbulence  creates  a  complex  heterogeneous  environment  of  bursting  vessels, 
fragmenting  equipment  and  aerosolized  liquids  within  the  blast  radius.  We  need  to 
characterize  this  spatially  and  temporally  non-uniform  thermal  environment  under 
extreme  conditions  such  as  at  hundreds  ofkilopascals  of  pressure  ( too  KPa),  and  in  very 
transient  atmospheres  that  change  every  microsecond  (1  ys).  Current  technologies 
cannot  provide  the  desired  information,  but  exploration  of  thermo-properties  of 
nanomaterials  could  lead  to  revolutionary  new  ways  to  measure  temperature  in  these 
extreme  environments.”  “Possible  research  directions  can  include  the  creation  and  study 
of  nano-composite  particles  that  change  and  maintain  material  phases  after  a  thermal 
event;  micro-encapsulated  materials  that  when  exposed  to  a  thermal  event  releases  the 
encapsulated  material  for  analysis,  etc.”1 

DTRA’s  objective  for  this  research  topic  was  “to  support  fundamental  research  to  identify 
and  characterize  the  means  and  materials  that  will  measure  the  temperature  and 
perhaps  also  retain  the  complete  thermal  history  of  events  (1  to  180  seconds)  involving 
extreme  conditions  ranging  from  tens  to  hundreds  ofkilopascals  of  pressure  (100  KPa), 
hundreds  of  degrees  of  temperature  ( maximum  700  °C),  and  microsecond  (1  ys)  changes 
of  these  conditions.  “x 

In  response  to  this  challenge,  we  proposed  to  characterize  and  evaluate  temperature- 
induced  irreversible  phase  changes  in  metal-oxide  nanoparticles  as  potential  temperature 
sensors  for  use  in  explosive  fireballs.  The  expected  phase  changes  are  often  relatively 
subtle  and  require  sensitive  means  for  read-out.  As  such,  we  proposed  to  dope  these 
materials  with  rare-earth  ions  such  as  europium  (Eu)  whose  optical  spectra  are  very 
sensitive  to  their  local  environment. 

The  scope  of  this  project  included: 

•  development,  characterization,  and  testing  of  several  temperature  sensor 
materials  and  demonstration  of  their  suitability  for  measuring  and  recording 
temperature  inside  an  explosive  fireball 

•  demonstration  that  these  sensor  materials  respond  at  least  as  fast  as  a  few 
milliseconds  and  possibly  faster 

•  development  of  a  laser-based  heating  and  3-color  pyrometer-based  temperature 
recording  system 

•  initial  demonstration  of  a  kinetic  analysis  approach  based  upon  using  two 
different  temperature  sensor  materials  to  extract  temperature  and  heating 
duration  information 

•  development  and  demonstration  of  tracking  particles  to  monitor  the  origin  of  the 
sensors  during  explosions 
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Our  approach  included: 

•  identifying  the  most  appropriate  synthesis  techniques  for  the  preparation  of 
metal-oxide  precursors  in  a  specific  state 

•  identifying  how  these  precursors  react  (i.e.  undergo  irreversible  phase  transitions) 
under  thermal  exposure 

•  developing  a  heating  technique  that  is  capable  of  reaching  temperatures  with 
heating  rates  comparable  to  those  occurring  in  an  explosion  in  a  controlled 
manner 

•  measuring  the  temperature  with  a  fast  enough  resolution  and  provide  feedback  to 
the  heating  mechanism 

The  results  gained  from  this  investigation  provide  the  required  information  to  develop 
further  the  kinetic  approach  in  combination  with  appropriate  sensor  mixtures  to  extract 
temperature  and  heating  duration.  Besides  the  use  as  a  development  tool  for  c-WMD,  the 
technique  might  find  more  general  applications  during  heating  events  such  as  fire  and 
arson  investigations.  We  also  envision  the  development  of  a  portable  analysis  unit  for  field 
testing. 
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III.  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 
A.  Synthesis 

This  section  summarizes  the  synthesis  techniques  used  for  the  preparation  of  the  sensor 
and  tracking  materials.  The  techniques  are  described  using  specific  example  materials, 
but  are  also  applicable  to  other  materials. 

p-Eu:Y203 

Irregular  sized  nanoparticles  of  precursor  of  Eu-doped  yttria,  p-Eu:Y203,  with  dimensions 
in  the  20-50  nm  range  are  synthesized  by  co-precipitation.  An  aqueous  mixture  of 
ammonium  hydroxide  and  ammonium  bicarbonate  (1:1  mole  ratio,  1.25  M)  is  drop  wise 
added  into  an  aqueous  solution  of  europium  nitrate/yttrium  nitrate  solution  (1:99  mole 
ratio)  until  the  pH  of  the  mixture  solution  reaches  7.5  -  8.0.  The  resulting  white 
precipitate  is  left  to  age  for  3  h.  Subsequently,  it  is  washed  with  deionized  water  and  rinsed 
with  acetone  before  it  is  placed  in  vacuum  to  dry  overnight  followed  by  heat  drying  at  150 
°C  for  12  h.  The  dried  white  cake  is  crushed  into  powder,  and  run  through  a  100  mesh 
sieve.2  The  europium  dopant  may  be  substituted  with  other  rare  earth  elements.  For  the 
preparation  of  tracking  particles,  p-Eu:Y203  is  heated  to  1400  °C  for  30  min. 

p-Eu:Zr02 

Irregular  sized  nanoparticles  of  precursor  of  Eu-doped  zirconium  oxide,  p-Eu:Zr02,  is 
obtained  by  adding  aqueous  solution  of  europium  nitrate/zirconium(IV)  oxychloride 
(1:99  mole  ratio,  0.15  M),  drop-wise  to  an  aqueous  solution  of  ammonium  hydroxide  (2.5 
vol.%)  until  a  pH  of  10.00  is  reached.  A  white  gel  forms  and  is  aged  for  4  h.  Subsequently, 
the  gel  undergoes  multiple  cycles  of  centrifugation  and  redispersion  in  deionized  water. 
The  washed  gel  is  then  rinsed  twice  with  acetone,  dried  in  vacuum,  and  then  dried  at  150 
°C  for  12  h.  Finally,  the  white  dried  cake  is  crushed  into  powder  with  mortar  and  pestles, 
and  run  through  a  100  mesh  sieve.3  The  europium  dopant  may  be  substituted  with  other 
rare  earth  elements. 

REi:Zr02/p-RE2:Y203  Core-shell 

Spherically-shaped  nanoparticles  of  precursor  of  RE-doped  zirconium  oxide,  p-REi:Zr02, 
with  about  300  nm  diameters  are  obtained  via  forced  hydrolysis  and  serve  as  the  core.  RE 
salt/zirconium(IV)  oxychloride  (1:99  mole  ratio,  0.2  M)  dissolved  in  an  isopropanol/ 
water  mixture  (5:1  vol/vol)  containing  hydroxypropyl  cellulose  (1.5  g/L)  is  heated  quickly 
in  a  microwave,  and  then  held  isothermally  in  a  hot  bath  at  75-80  °C  for  30  min.  The 
formed  spherical  nanoparticles  are  stabilized  by  adding  ammonium  hydroxide  solution 
to  neutralize  the  solution’s  pH.  The  nanoparticles  are  then  washed  with  deionized  water, 
rinsed  with  acetone,  dried  in  vacuum,  and  then  dried  at  150  °C  for  12  h.  The  nanoparticles 
are  then  heated  to  1000  °C  for  30  min  to  transform  them  into  metallic  oxides. 

The  precursor  of  RE-doped  yttrium  oxide,  p-RE2:Y203,  is  grown  via  homogeneous 
precipitation  using  urea  in  the  presence  of  REi:Zr02  (5  mg/mL  concentration)  to  yield  the 
REi:Zr02/p-RE2:Y203  core-shell  nanoparticles.  The  total  metallic  salt  concentration  is  5 
mM  and  the  urea  concentration  is  0.2  M.  The  homogeneous  precipitation  is  done  at  90 
°C  for  3  h.  The  obtained  nanoparticles  are  then  washed  with  deionized  water,  rinsed  with 
acetone,  dried  in  vacuum,  and  the  dried  at  150  °C  for  12  h.4 
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B.  Lab-based  heating  and  temperature  measurements 
Furnace 


Initially,  irreversible  phase  transitions  of  the  precursors,  placed  in  crucibles,  were 
induced  by  heating  in  a  box  furnace  (Lindberg  Blue).  This  heating  method  allowed  us  to 
subject  the  precursors  to  well-defined  conditions  and  to  follow  the  phase  changes  over 
time.  The  precursors  which  typically  consist  of  amorphous  carbonates  or  hydroxides 
undergo  a  thermal  decomposition  process  that  leads  to  the  formation  of  oxides.  This 
process  can  occur  over  several  steps.  Depending  on  the  intermediary  products,  nucleation 
may  occur  while  decomposition  is  still  undergoing,  or  nucleation  may  start  only  after  the 
decomposition  process  is  complete.  Under  continued  heating,  the  nucleated  particles  will 
then  start  to  grow  into  larger  structures. 

Pyroprobe 


In  order  to  identify  phase  changes  that  occur  under  much  faster  heating  rates,  we  used  a 
Pyroprobe-iooo  (CDS  Analytical,  Inc.),  see  Figure  1.  This  pyroprobe  is  capable  of  heating 
rates  from  10  °C/s  up  to  20000  °C/s,  temperatures  up  to  1400  °C,  and  heating  durations 
from  0.01  s  up  to  99.99  s.  The  sample  material  can  be  placed  on  a  heating  strip  or  inside 
a  quartz  tube  which  is  then  inserted  into  a  heating  coil.  Placing  our  material  on  the  heating 
strip  did  not  work  well.  The  material  would  typically  come  off  during  the  heating  and  thus 
heating  durations  were  not  well  defined.  The  development  of  gas  products  such  as  C02 
during  the  thermal  decomposition  is  probably  responsible  for  this  effect. 


|J  win! I 


fiiimi,'  riibw 


r  ±—3 


Figure  1.  CDS  Analytical 
pyroprobe  heater,  heating  strip, 
and  heating  coil. 


Time  (s) 

Figure  2.  Actual  pyroprobe  temperature 
profiles. 


While  placing  our  material  into  the  quartz  tubes  and  heating  it  in  the  coil  kept  the  material 
in  place,  we  noticed  that  heating  the  material  at  the  center  of  the  tube  to  the  set 
temperature  often  took  a  significant  amount  of  time.  We  embedded  a  thermocouple 
(Omega  CHAL-005)  with  a  fast  response  time  into  the  sensor  material  and  monitored  the 
temperature  during  the  heating  process.  Even  though  we  used  a  heating  rate  of  20000 
°C/s,  the  actual  heating  process  took  much  longer,  see  Figure  2.  With  a  heating  duration 
of  10  s,  none  of  the  samples  reached  the  set  temperature. 
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Laser  heating 

In  order  to  achieve  faster  heating  rates,  we  used  a  CO2  laser  (Synrad  fioo).  Initially,  we 
enclosed  our  sample  material  between  graphite  sheets  to  hold  the  material  in  place.  Using 
this  approach,  we  improved  our  heating  process  and  achieved  effective  heating  rates  as 
fast  as  400  °C/s,  see  Figure  3. 

To  further  increase  the  heating  rate  while  keeping  the  sample  in  place,  we  are  cold¬ 
pressing  the  sensor  material  into  a  pellet,  place  it  into  a  sample  holder,  and  cover  it  with 
an  IR  transparent  window  such  as  ZnSe.  We  then  use  the  C02  laser  and  a  3-color 
pyrometer  with  home-build  control  unit  and  feedback  loop  to  heat  the  sample  in  a 
controlled  manner.  Using  this  approach,  we  are  able  to  reach  heating  rates  up  to  26600 
°C/s,  see  Figure  3.  At  this  rate,  it  takes  about  30  ms  to  reach  1000  K.  A  fixable  problem 
with  the  controller  is  currently  limiting  the  isothermal  heating  time  and  causes 
fluctuations  in  the  temperature  value. 


Figure  3.  Actual  indirect  (left)  and  direct  (right)  laser  heating  temperature  profiles. 


Pyrometer 

Our  3-color  pyrometer  is  shown  in  Figure  4. 
Fundamentally  a  pyrometer  works  based  on  Planck's 
law  of  blackbody  radiation,  which  states  that  the 
emitted  spectral  radiance  of  a  black  body  at  wavelength 
A  and  temperature  T  is  given  by: 


IU,T)  = 


2  he2 


1 


^5  ehc/AkT  _  1  ’ 


(1) 


where  h  is  Planck's  constant,  c  is  the  speed  of  light,  and 
k  is  Boltzmann's  constant.  While  Equation  1  holds  for  a 
perfect  black  body,  in  reality  most  materials  are 
actually  'grey-bodies'  meaning  that  they  have  a 
wavelength  and  temperature  dependent  emissivity 
s(X,T)  Additionally,  what  we  measure  experimentally  is 

not  the  spectral  radiance,  but  a  photodiode  voltage  which  is  proportional  to  the  spectral 
radiance  integrated  over  a  narrow  wavelength  range,  with  the  proportionality  constant 
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Figure  4.  3-color  pyrometer. 


depending  on  geometry,  optical  properties  of  the  pyrometer,  and  the  quantum  efficiencies 
of  the  photodiodes.  For  simplicity  we  lump  the  various  wavelength  responses  into  one 
function  Si(X).  Taking  into  account  the  emissivity,  wavelength  response  of  the  pyrometer, 
and  spectral  integration  we  find  that  the  voltage  measured  by  one  of  the  photodiodes  at 
temperature  T  is 


Vi(T) 


r 

JAi 


' Ai+A/ 2  2ftc2  S(d)e(d,  T) 

Aj-A/2  d5  ~  1 


dA. 


(2) 


While  Equation  2  describes  the  photodiode  voltage  measured  for  a  sample  at  temperature 
T,  what  we  want  in  practice  is  to  invert  Equation  2  to  find  the  temperature  as  a  function 
of  the  photodiode  voltage. 

In  its  current  form  Equation  2  is  too  computationally  intensive  to  implement  a  simple 
inversion  algorithm,  which  is  necessary  for  the  pyrometer  to  function  quickly.  We 
therefore  make  several  approximations  to  drastically  simplify  the  computation  required. 
First  we  assume  that  XkT  «  he  such  that  Planck's  law  can  be  replaced  by  Wien's 
approximation: 

HT)  =  2 he2A~5e~hc,AkT.  (  3  ) 

This  approximation  is  valid  for  our  wavelength  and  temperature  range,  as  the  maximum 
deviation  between  Equation  l  and  3  is  only  0.22  %  in  our  experimental  regime. 

Next  we  assume  that  for  the  wavelength  range  we  use  the  emissivity  depends  only  on 
temperature,  which  is  a  fairly  good  assumption  for  ceramics  and  the  narrow  wavelength 
range  we  are  using.  With  these  approximations  Equation  2  becomes 

JrAi+ A/2 

S(A)A~5e~hc/AkTdA. 

Ai- A/2  (  4  ) 

Finally  we  assume  that  the  filter  bandwidth  is  negligible  such  that  the  signal  at  each 
detector  is  due  to  a  single  wavelength. 

To  calculate  the  approximate  temperature  we  take  the  ratio  of  the  voltage  at  the  first  and 
second  detector: 


V'l 

Si 

M2\5 

he  A2  -  di 

S2 

\m)  exp 

kT  AxA2 

Where  Si  and  S2  are  the  detector  response  constants  of  the  first  and  second  detector 
respectively.  Rearranging  Equation  5  we  find  the  inverse  temperature  as  a  function  of  the 
voltage  ratio: 

1  ^  k  A1A2  ^ 

T  he  A2  —  A]  n 


O2S2 

s, 

U; 

Equation  6  relates  the  inverse  temperature  to  the  ratio  of  the  detector  voltages  with  the 
assumption  that  the  emissivity  does  not  change  with  wavelength.  This  assumption  is  not 
necessarily  valid  and  therefore  to  increase  our  accuracy  we  move  to  the  next  order  of 
approximation,  which  is  to  include  an  approximate  emissivity  wavelength  dependence. 
One  model  of  the  wavelength  dependent  emissivity  is  an  exponential  expansion  given  by: 

8 


e(d,  T)  =  exp(-ao  -  a\A  -  a2A2  -  . . .)  (7) 

Where  ao  defines  the  magnitude  of  the  emissivity  and  the  higher  order  terms  determine 
the  emissivities  spectral  shape.  Since  we  are  using  a  relatively  narrow  wavelength  range, 
it  is  sufficient  to  only  use  the  first  two  orders,  which  gives  a  wavelength  dependent 
emissivity  of  the  form: 

e(d,  T )  as  e~aoe~a'A 

=  eoe~a,A  (  8  ) 


With  this  correction  we  can  write  the  first  order  single  wavelength  detector  voltage  as 

V,-  =  2 S jhc2eoA~5  exp  {-hc/AjkT  -aid;}.  (  g  ) 

Since  we  now  have  two  unknowns,  T  and  ai,  we  need  to  use  two  separate  ratios  in  order 
to  calculate  the  temperature  from  the  pyrometer.  We  therefore  use  the  ratio  of  detector  1 
and  2,  and  the  ratio  of  detector  1  and  3,  which  are  given  by: 


Yi 

Vl 


r\ 3 


-  Yi  =  h(di 


V3 


/  d2' 

5 

\dl , 

exp 

(A3) 

5  r 

UJ 

1  exp 

he  A2  —  A\ 
kJ  A 1/I2 

he  d3  -  A\ 
kT  d,d3 


+  ai(d2  -  di) 
+  a\(A3  -  di) 


(10) 


(11) 


To  calculate  the  inverse  temperature  we  first  rearrange  Equation  10  to  give  the  inverse 
temperature  dependent  on  m, 


1  k  d,d2 

/ 

In 

rl2S2 

f-f 

T  he  d2  -  A\ 

111 

V 

5, 

\d2  / 

+  «i(di  —  A2)  . 


(12) 


We  then  substitute  Equation  12  into  Equation  11  and  solve  for  ai, 

\5 

1  ri'irs  7  /  a  1 

a  1  = - - —  < - In 


/l3  -  /12  M3  -  ^1 


^13^  3  / 


1 


A2 


A3U2  -  A\) 


In 


r\2S2  (A]_ 
*2 


(13) 


Finally,  we  substitute  Equation  13  into  Equation  12  and  simplify  to  give  the  inverse 
temperature  in  terms  of  the  detector  voltages: 


T 


kA\A2A3 


Mdl-d2)(di-d3)(d2-d3) 


(di-d3)ln 


Vi52/di 
F25i  U2 


+  (d2  -di)ln 


V\S  3 

/dlfl 

V3S 1 

\d3  / 

(14) 


One  of  the  main  assumptions  in  deriving  Equation  5  is  that  the  detectors  only  measure  a 
single  wavelength,  which  is  not  true.  In  order  to  estimate  the  error  in  the  voltage  ratios 
due  to  this  assumption  we  calculate  the  voltage  at  each  detector  for  a  nonzero  filter 
bandwidth  using  Equation  4  with  the  emissivity  assumed  to  be  constant.  Table  1  tabulates 
the  voltage  ratios  for  a  12  nm  and  50  nm  filter,  as  well  as  the  percent  difference  between 
the  multiwavelength  calculations  and  the  single  wavelength  approximation.  From  these 
percent  differences  we  see  that  ratio  errors  are  less  than  1%. 
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Table  l.  Calculated  intensity  ratios,  and  their  percent  differences,  for  different  filter 
bandwidths  at  six  different  temperatures.  The  percent  difference  represents  that 
percentage  that  the  ratio  changes  if  we  make  the  single  wavelength  approximation. 


T  (K) 

Wavelength  (nm) 

Ratio  (Vx/Vi5oo) 

Percent  Difference 

0  nm 

12  nm 

50  nm 

12  nm 

50  nm 

473 

1550 

1.63281 

1.63291 

1.63117 

0.00601 

0.10064 

1600 

2.57262 

2.57262 

2.56631 

-6E-05 

0.24526 

673 

1550 

1.34430 

1.34419 

1-34343 

-0.00814 

0.06469 

1600 

1.76512 

1.76494 

1.76345 

-0.01015 

0.09479 

873 

1550 

1.20990 

1.20126 

1.20961 

-0.71389 

0.02385 

1600 

1.43926 

1.42898 

1.43865 

-0.7145 

0.04251 

1073 

1550 

1.13258 

1.13258 

1-1325 

-0.00053 

0.00800 

1600 

1.26642 

1.26642 

1.26625 

-0.00066 

0.01396 

1273 

1550 

1.08252 

1.08253 

1.08253 

0.00036 

-0.0003 

1600 

1.16024 

1.16024 

1.16028 

-7.9E-05 

-0.00315 

1473 

1550 

1.04758 

1.04764 

1.04748 

0.00472 

0.01009 

1600 

1.08882 

1.0888 

1.08887 

-0.00274 

-0.00353 

For  the  actual  implementation  of  the  pyrometer,  Equation  14  is  problematic  due  to 
spectral  integration  of  the  detectors  and  deviations  between  theory  and  reality.  A  more 
helpful  approach  is  to  use  an  empirical  model  of  the  voltage  ratios  using  a  fit  function  of 

r  =  AeblT,  (  15  ) 

where  r  is  a  given  ratio,  T  is  the  temperature,  and  both  A  and  b  are  fit  parameters.  We  can 
then  invert  Equation  15  to  find  the  inverse  temperature  for  a  given  voltage  ratio 

H"0- 

In  order  to  use  Equation  16  for  the  pyrometer  we  must  first  perform  calibration 
measurements  in  order  to  determine  A  and  b  for  each  ratio  permutation.  We  therefore 
measure  the  emission  from  a  pyroprobe  at  different  temperatures  and  obtain  the  detector 
voltages  as  a  function  of  temperature1,  shown  in  Figure  5.  From  these  voltages  we 
calculate  the  voltage  ratios,  shown  in  Figure  6,  and  fit  them  to  Equation  15,  with  Table  2 
listing  the  fit  parameters  for  each  ratio  permutation. 


1  For  T  <  923  °K  the  detectors  use  a  gain  setting  of  70  db.  For  higher  temperatures  the  gain  has  to  be  set  lower  (30 
dB)  to  avoid  saturation. 
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Temperature  (°K) 

Figure  5.  Detector  voltages  as  a  function 
of  temperature. 


600  650  700  750  800  850  900 


Temperature  (°K) 

Figure  6.  Measured  voltage  ratios  as  a 
function  of  pyroprobe  temperature  and 
fits  to  Equation  15. 


Table  2.  Calibration  parameters  determined  for  pyrometer  using  pyroprobe  heating 
and  predicted  b  value  from  theory. _ 


A 

B(°K) 

b  i  h  (°K) 

Vi/Va 

5.93  ±  0.20 

-330  ±  29 

-309.3 

V1/V3 

5.72  ±  0.15 

-553  ±  23 

-599-5 

V2/V3 

0.967  ±  0.035 

-225  ±  32 

-289.9 

Using  Equation  16  and  the  parameters  for  each  ratio  we  can  compute  three  temperature 
values,  which  can  be  averaged  to  give  a  single  temperature.  To  test  this  method  we  use 
our  data  from  the  pyroprobe  to  calculate  the  pyroprobes  temperature.  Table  3  lists  the 
pyroprobe’s  set  temperature,  pyrometer  voltage  ratios,  and  the  calculate  temperature. 

From  Table  3  we  see  that  the  pyrometer  calculates  the  temperature  far  more  precisely  and 
accurately  at  high  temperatures,  which  is  expected  as  the  signal-to-noise  ratio  is  bigger  at 
higher  temperatures.  Figure  7  shows  the  difference  magnitude  in  the  calculated  and  set 
temperatures  as  a  function  of  set  temperature.  From  Figure  7  we  see  that  as  the  set 
temperature  increases  the  difference  drastically  decreases  to  near  zero  and  also  that  all 
temperature  differences  are  within  uncertainty  of  zero.  The  slight  increase  in  the 
difference  at  higher  temperatures  is  due  to  the  high  T  measurements  requiring  a  different 
detector  gain  setting. 
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Table  3.  Measured  voltage  ratios  at  different  pyroprobe  temperatures  and  calculated 


temperature  from  pyrometer. 


T(°K) 

V1/V2 

V1/V3 

Tcalc  (°K) 

523 

4-5  ±  2.3 

2.22  ±  0.64 

708  ±  574 

573 

3-4  ±  14 

2.08  ±  0.50 

573  ± 126 

623 

3.41  ±  0.44 

2.29  ±  0.20 

601  ±  42 

673 

3.56  ±  0.12 

2.465  ±  0.057 

652  ±  13 

773 

3.871  ±  0.024 

2.803  ±  0.012 

775.3  ±  3.8 

823 

3.968  ±  0.034 

2.912  ±  0.023 

821.0  ±  6.3 

839 

4.004  ±  0.056 

2.964  ±  0.041 

841  ±  11 

873 

4.064  ±  0.030 

3.038  ±  0.022 

874.8  ±  6.4 

923 

4.143  ±  0.013 

3.141  ±  0.010 

922.4  ±  3.0 

1116 

4429  ±  0.017 

3.507  ±  0.013 

1130  ±  19 

1201 

4.511  ±  0.064 

3.622  ±  0.051 

1210  ±  40 

1473 

4.724  ±  0.019 

3.927  ±  0.015 

1470  ±  18 

Set  Temperature 

Figure  7.  Difference  in  calculated  and  set 
temperatures  as  a  function  of  set 
temperature.  As  the  temperature  increases 
the  difference  between  the  set  and  calculated 
temperatures  drastically  decreases  and  all 
differences  are  within  uncertainty  of  zero. 
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C.  Explosion  tests 

During  the  course  of  this  project,  a  four  series  of  explosion  tests  were  conducted  to  test 
the  performance  of  the  sensors.  These  tests  were  performed  at  the  Naval  Surface  Warfare 
Center,  Indian  Head  Explosive  Ordnance  Disposal  Technology  Division  (NSWC 
IHEODTD).  The  tests  were  conducted  either  in  a  large  room  or  inside  a  smaller  steel 
chamber  using  aluminized  and  non-aluminized  charges.  Since  all  of  the  tests  were  slightly 
different,  the  specifics  are  described  in  the  Results  and  Discussion  section. 
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D.  Characterization 

A  variety  of  optical  characterization  techniques  were  used  to  characterize  the  optical 
properties  of  the  sensor  materials  and  to  identify  appropriate  optical  signatures  that 
would  provide  information  about  the 
temperature  that  these  sensors  were 
exposed  to.  The  bulk  of  these 
measurements  was  based  upon 
fluorescence  and  fluorescence-excitation 
measurements.  During  fluorescence 
measurements,  the  sample  material  was 
illuminated  at  a  particular  wavelength 
using  a  laser  or  LED  source  and  the  light 
emitted  by  the  sample  is  monitored  as  a 
function  of  wavelength.  In  fluorescence- 
excitation  measurements,  the 
monitoring  monochromator  was  set  to  a 
specific  wavelength  and  the  excitation 
wavelength  (laser)  was  tuned  over  a 
certain  range.  The  resulting  spectra  were 

then  analyzed  with  respect  to  full-width-half-maximum  (FWHM)  of  peaks,  peak 
positions,  and  relative  peak  intensities.  Figure  8  shows  the  schematic  of  the 
fluorescence/fluorescence-excitation  setup.  The  optical  spectroscopy  measurements  are 
performed  using  a  Continuum  Nd:YAG  laser  operating  at  10  Hz  in  the  3rd  harmonic  at  355 
nm.  This  laser  pumps  a  Continuum  Panther  Optical  Parametric  Oscillator  (OPO).  Using 
several  mirrors,  laser  light  of  about  10  ns  pulse  length  of  the  desired  wavelength  is 
directed  onto  the  sample  material  which  is  placed  in  a  vial,  see  Figure  6  (right).  The 
emitted  fluorescence  is  collected  using  a  2”  diameter  lens  with  a  focal  length  of  15  cm. 
Using  a  2”  diameter  lens  with  a  focal  length  of  50  cm,  the  fluorescence  is  then  focused 
onto  the  entrance  slit  of  an  Acton  2750  monochromator/spectrometer  with  PMT  and  CCD 
array  attached.  The  excitation  spectra  are  acquired  using  a  SpectraHub  (Acton)  interface 
and  in-house  Labview  software.  The  fluorescence  spectra  are  acquired  using  the  CCD 
camera  and  WinSpec  software  (Princeton  Instruments).  The  lifetimes  are  measured  using 
a  PMT  and  an  oscilloscope  (Tektronix  TDS3052B).  The  data  is  further  processed  using 
Origin  Pro  8.5. 

Additional  characterization  techniques  used  during  this  project  include  Fourier 
Transform  Infrared  (FTIR),  x-ray  diffraction  (XRD);  scanning  electron  microscopy 


Figure  8.  Schematic  of  the  fluorescence/ 
fluorescence-excitation  setup. 
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(SEM),  energy  dispersive  x-ray  spectroscopy  (EDS),  transmission  electron  microscopy 
(TEM),  and  differential  scanning  calorimetry  (DSC). 

Fourier-transform-infrared  (FT-IR)  spectra  are  measured  using  a  Cary  680  FTIR 
(Varian/Agilent).  The  nanopowder  morphologies  are  characterized  using  a  scanning 
electron  microscope  (SEM),  NOVA  nanoSEM  230  (FEI),  and  a  transmission  electron 
microscope  (TEM),  Tecnai  G2  20  Twin  TEM  (FEI).  The  SEM  is  equipped  with  X-Max 
Silicon  Drift  Detector  (SDD)  with  an  80  mm2  active  area  (Oxford  Instruments,  PLC.)  for 
energy  dispersive  spectroscopy  (EDS)  measurements.  A  PANalytical  X’Pert  Pro 
diffractometer  (PANalytical  B.V.),  using  Cu-Ka  radiation  (A=i.54i8  A)  and  operated  at  45 
kV  and  40  mA,  is  used  to  measure  x-ray  diffraction  patterns.  Diffraction  pattern  were 
measured  using  a  PIXceD0  detector  (PANalytical  B.V.).  The  x-ray  beam  was  collimated 
using  a  fixed  divergence  slit  (FDS)  with  0.04  radian  Soller  slits,  0.50  divergence  slit,  and 
a  10  mm  mask.  The  x-ray  beam  was  monochromatized  using  a  X’Celerator 
monochromator  (PANalytical  B.V.).  Due  to  the  age  of  the  x-ray  source,  it  also  contained 
W-La  and  Cu-K(3  radiation  that  led  to  small  ghost  peaks  in  the  XRD  spectra. 

Transmission  electron  microscopy  (TEM)  images  are  obtained  using  a  Tecnai  G2  20 
Twin  TEM  (FEI)  or  a  Philips  CM  200  transmission  electron  microscope  (TEM).  ANetzch 
STA  409  PCat,  using  a  ramp  rate  of  20°C/minute  under  Ar  atmosphere,  is  used  to 
measure  differential  scanning  calorimetry  (DSC)  and  thermogravimetric  analysis  (TGA). 
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IV.  RESULTS  AND  DISCUSSIONS 


This  section  is  divided  as  follows:  Section  A  describes  our  modeling  results  for  the  kinetics 
of  temperature-induced  phase  changes;  Section  B  describes  various  temperature  sensors 
we  characterized,  but  focuses  on  Eu:Y203  and  Eu:Zr02,  the  two  most  important  ones; 
Section  C  describes  our  temperature/time  (thermal  impulse)  sensors;  and  Section  D 
describes  the  tracking  particles. 

A.  Phase  Change  Modeling 

General  Model 

We  can  model  the  kinetics  of  thermal  degradation  by  considering  the  degree  of 
conversion,  a,  which  is  found  to  obey  the  differential  equations 

^  =  h(T)f(a), 
at 

where  h(T)  is  an  Arrhenius  equation, 

h(T)  =  Ae\ p  — 

kT  i 

and /(a)  can  generally  be  expressed  as: 

f(a)  =  am(l-a)n\-\n(l-a)f  (19) 

where  m,  n,  andp  depend  on  the  reaction  mechanism.  Substituting  Equations  18  and  19 
into  Equation  17  lets  us  write  the  whole  kinetic  equation  as: 

am{\  -a,)”[-ln(l  —  a)\p.  (20) 


da 

—  =  A  exp 
dt 


mo 


(17) 

(18) 


Moving  all  the  a  and  t  terms  to  opposite  sides  we  can  solve  Equation  20  to  give  a  general 
solution  of: 


da' 

ar'm(l  -<*')"[- ln(l  -  a')F 


E 


dt'. 


(21) 


At  this  point  further  simplification  of  the  general  case  cannot  be  done  as  both  the  integrals 
depend  on  the  specific  functional  forms,  some  of  which  have  analytical  solutions,  and 
others  which  require  numerical  solutions. 


As  an  example  of  a  semi-analytic  solution  (we  leave  T(t)  in  its  general  form),  we  use  the 
case  of  random  nucleation  and  growth,  which  gives  exponent  values  of  m  =  o,  n  =  1,  p  = 
o.  Assuming  that  the  system  begins  with  no  conversion,  i.e.  a0  =  o,  we  can  solve  Equation 
21  to  give  the  degree  of  conversion  as  a  function  of  time  as: 


a(t )  =  1  -  exp 


(22) 


The  integral  in  Equation  22  has  no  general  analytical  solutions  and  typically  requires 
numerical  methods  to  calculate.  However,  the  simplest  case  of  a  temperature  profile,  that 
of  an  isotherm,  gives  an  exact  analytical  solution  of 
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a(t)  =  1  -  exp  <  -A  exp 


(23) 


kT: 


if  -  to) 


Nucleation  and  Growth  Example 

To  better  understand  how  different  temperature-time  profiles  affect  the  degree  of 
conversion,  we  assume  a  nucleation  and  growth  mechanism,  with  the  degree  of 
conversion  as  a  function  of  time  given  by  Equation  22.  For  our  simulations  we  assume 
parameter  values  of  E  =  0.725  eV  and  A  =  io4,  which  are  on  the  same  order  of  magnitude 
as  those  measured  by  Turmanova  et  al.  for  a  different  materials  We  will  first  consider 
the  effect  of  heating  rate  on  the  degree  of  conversion  and  then  consider  an  estimated 
explosion  temperature  profile. 

Effect  of  Heating  Rate 

We  simulate  the  effect  of  the  heating  rate,  g,  on  the  degree  of  conversion  by  assuming  a 
heating  profile  consisting  of  a  linear  ramp  from  273  °K  to  873  °K,  followed  by  a  1  s 
isotherm  and  then  an  exponential  cooling  profile  with  a  time  constant  of  x  =  1  min.  For 
the  ramp  rates  we  use  200  °K/s,  400  °K/s,  600  °K/s,  and  2400  °K/s,  with  the  ramping 
time  varying  such  that  the  final  temperature  is  always  873  °K.  Figure  9  shows  the 
temperature  profile  and  degree  of  conversion  as  a  function  of  time  for  the  600  °K/s 
heating  rate,  with  the  degree  of  conversion  given  both  by  Equation  22  and  by  the 
isothermal  model. 


From  Figure  9  we  see  that  for  the  isothermal  model  the  degree  of  conversion  only  changes 
during  the  isothermal  time,  neglecting  the  effects  of  heating  and  cooling.  By  comparing 
the  full  model  to  the  isothermal  model  we  see  that  the  heating  and  cooling  works  to 
increase  the  degree  of  conversion  further  than  just  for  the  isothermal  period. 


Figure  9.  Temperature  profile  as  a 
function  of  time  plotted  with  the  degree  of 
conversion  using  the  full  model  of 
Equation  22  and  the  isothermal  model, 
which  ignores  heating  and  cooling. 


Tima  (s) 

Figure  10.  Degree  of  conversion 
predicted  by  Equation  23  as  a  function  of 
time  for  different  heating  rates,  as  well  as 
the  isothermal  value.  Note  that  as  the 
heating  rate  decreases  the  full  model 
diverges  more  from  the  isothermal 
model. 
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Next  we  plot  the  degree  of  conversion  from  Equation  22  for  different  ramp  rates  as  shown 
in  Figure  10,  as  well  as  the  full  isotherm  value.  From  Figure  10  we  find  that  as  the  ramp 
rate  increases,  the  final  degree  of  conversion  begins  to  approach  the  isothermal  value. 
This  can  be  understood  as  the  isothermal  calculation  essentially  assumes  an  infinite  ramp 
rate  and  cooling  rate.  This  implies  that  the  heating  time  and  cooling  time  will  increase  the 
degree  of  conversion  and  that  slower  heating/cooling  rates  will  have  a  larger  impact  on 
the  degree  of  conversion. 

Additionally,  Figure  10  shows  that  there  is  an  error  associated  with  the  degree  of 
conversion  when  neglecting  the  heating  rate  and  only  assuming  an  isothermal  profile.  To 
better  understand  this  error  and  how  different  experimental  parameters  affect  it,  we 
define  the  fractional  error  between  the  isothermal  assumption  and  full  model  to  be: 


FE  _  Fuii 


iso 


Full 


(24) 


where  ciFuii  is  the  degree  of  conversion  calculated  by  the  full  model  and  aiso  is  the  degree 
of  conversion  calculated  using  the  isothermal  assumption.  Equation  24  essentially  gives 
us  a  measure  of  how  good,  or  poor,  the  isothermal  assumption  is.  As  an  example,  Figure 
11  shows  the  degree  of  conversion  for  the  isothermal  assumption  and  the  full  model  as  a 
function  of  ramp  rates  and  overlayed  is  the  fractional  error.  From  Figure  11  we  see  that  as 
the  ramp  rate  increases  the  fractional  error  decreases,  which  is  consistent  with  the 
observations  in  Figure  10. 


Heating  Rate  (*K/s) 

Figure  11.  Degree  of  conversion  and 
fractional  error  as  a  function  of  heating 
rate. 


Furthermore,  in  Figure  12  we  plot  the 
fractional  error  as  a  function  of  ramp  rates 
for  different  (a)  isothermal  times  and  (b) 
isothermal  temperatures.  From  Figure  12 
we  find  that  for  long  times  and  high 
temperatures  the  fractional  error  is  small. 
However,  for  short  times,  or  lower 
temperatures  the  fractional  error  can  be 
substantial,  up  to  82%. 

From  these  results  it  is  obvious  that 
achieving  the  highest  possible  ramp  rate  is 
of  utmost  importance  to  obtain  an  accurate 
degree  of  conversion. 
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Heating  Rate  (*K/s)  Heating  Rate  (*K/s) 

Figure  12.  Fractional  error  as  a  function  of  heating  rate  for  different  (left)  isothermal 
times  and  (right)  isothermal  temperature. 

Effect  of  Isothermal  Time 

In  addition  to  considering  the  effect  of  the  heating  rate  on  the  degree  of  conversion,  we 
also  consider  how  changing  the  isothermal  time  affects  the  degree  of  conversion  and  the 
fractional  error  due  to  the  isothermal  approximation.  Figure  13  shows  the  degree  of 


WK) 

Figure  13.  Degree  of  conversion  and 
fractional  error  as  a  function  of 
isothermal  time  for  a  heating  rate  of  200 
°K/s  and  an  isothermal  time  of  100  ms. 


conversion  for  the  full  model  and  isothermal 
model,  as  well  as  the  fractional  error  as  a 
function  of  isothermal  temperature  for  a 
heating  rate  of  200  °K/s  and  an  isothermal 
time  of  100  ms.  From  Figure  13  we  see  that 
for  low  isothermal  temperatures  (T  <  900 
°K)  the  fractional  error  is  in  between  0.7  and 
0.8,  but  decreases  rapidly  as  the 
temperature  exceeds  900  °K.  The  fractional 
error  remains  above  0.1  until  around  1300 
°K,  implying  that  for  isothermal 
temperatures  below  1300  °K  the  isothermal 
approximation  is  poor. 

To  more  fully  get  a  sense  of  the  effect  of  the 
isothermal  time  we  consider  the  fractional 
error  as  a  function  of  isothermal 
temperature  for  different  times  and  heating 
rates,  as  shown  in  Figure  14. 
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Figure  14.  Fractional  error  as  a  function  of  isothermal  temperature  for  different 
isothermal  times,  At,  and  two  different  heating  rates.  As  the  isothermal  time  decreases 
the  fractional  error  increases. 

“Exact”  Solution 

For  a  temperature  profile  consisting  of  a  linear  ramp  from  T0  to  TiS0  at  a  rate  g,  followed 
by  an  isotherm  of  time  dt  and  concluding  with  an  instantaneous  cooling,  we  can  solve 
Equation  19  exactly,  albeit  with  special  functions.  We  begin  with  Equation  22 


a(t )  =  1  -  exp 


-Af" 


exp 


k(T0  +  gt') 

where  the  integral  in  Equation  is  given  by, 


dt'  -  A  exp 


kTi 


Atj; 


(25) 


f 

Jt0 


k(T0  +  gt') 

1 


dt' 


=  -j-  yk(T0  +  gt )  exp 
=  ^  Ut iso  exp 


kT; , 


k(T0  +  gt) 
-Ei 


-Ei 


1 


kT; 


k(T0  +  gt) 
-  k7\)  exp 


-  kTo  exp 


E 

kTo 


E 

kTo. 


+£,y}- 


(26) 


where  the  function  Ei(x)  is  the  exponential  integral  function. 


a(t)  =  1  -  exp 


-u 


kT  iso  exp 


kT, , 


-Ei 


-1 

kTi,, 


-  kTo exp 


E 

'kTo 


A  exp 


kTh 


(27) 


19 


Explosion  Temperature  Profile 


Tima  (ms) 

Figure  15.  Temperature  and  degree  of 
conversion  during  simulated  explosion 
profile. 


To  simulate  the  degree  of  conversion  as  a 
function  of  time  during  an  explosion,  we 
assume  an  explosion  temperature  profile  of, 

T{t)  =  T0  +  tATe~,lT ,  (  2g  ) 

Where  7b  is  the  initial  temperature,  t  is  the 
time,  AT  is  an  amplitude  term,  and  x  is  a 
|  decay  time.  We  also  use  the  same  A  and  E 
0  values  used  in  the  previous  section.  Figure 
15  shows  the  temperature  profile  and  degree 
of  conversion  as  a  function  of  time  during 
the  explosion.  From  Figure  15  we  see  that 
the  majority  of  the  change  in  the  degree  of 
conversion  occurs  during  heating  and 
slightly  into  the  cool  down. 


Temperature  and  Time  for  Isothermal  Nucleation  and  Growth 

The  proposed  goal  for  this  project  is  to  be  able  to  extract  a  temperature  and  time  from 
thermally  induced  changes  in  our  samples.  What  this  means  precisely  for  an  explosion  is 
ambiguous  as  in  explosion  the  temperature  is  constantly  changing.  However,  for 
simplistic  modeling  purposes,  we  can  assume  an  isothermal  temperature-time  profile 
with  the  temperature  instantly  becoming  TiS0  for  a  period  of  time  At,  after  which  it 
instantly  cools  back  to  ambient  temperature.  Assuming  that  we  have  two  materials  in  our 
samples  (with  nucleation  and  growth  parameters  Ai,  Ei  and  A2,  E2 )  we  can  write  the 
isothermal  degree  of  conversion  for  the  two  materials  as: 


ad/)  =  1  -expj 

j-Aj  exp 

El 

kTiSO 

[  (29) 

a2(/)  =  1  -  exp  | 

j-A2  exp 

e2 

kTiS() 

|  (30) 

Performing  some  simple  algebra  on  Equations  29  and  30  we  can  eliminate  the  time 
dependence  and  determine  the  isothermal  temperature  in  terms  of  the  materials 
parameters  and  the  measured  degree  of  conversion,  with  the  inverse  isothermal 
temperature  given  by: 


1  _  k  1  A2  ln(l  -  aj) 

Tiso  E2-E\  _Ajln(l-a2) 


(31) 


Finally,  by  substituting  Equation  31  into  Equation  29  and  rearranging  we  can  find  the 
isothermal  time  as, 


A/  =  -  ln(  1  -  a|)Ai  exp 


E 1 


E2  -  E\ 


In 


A2  ln(l  -  ai) 


A  \  ln(l  -  a2 ) 


(32) 
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Example  Calculation 

To  better  understand  what  happens  when  we  apply  the  isothermal  approximation  to  a 
non-isothermal  temperature  profile  we  start  with  the  temperature  profile  given  by 


Equation  28,  with  the  parameters  varied  to 
give  different  peak  temperatures  and  peak 
widths.  We  then  calculate  asymptotic  values 


0  20  40  60  80  100 

Tima  (ms) 

Figure  16.  Temperature  and  degrees  of 
conversion  as  a  function  of  time  for  a 
temperature  profile  given  by  Equation 
28. 


§  of  ai  and  a2  for  the  temperature  profiles, 
|  from  which  we  use  Equations  31  and  32  to 
1  calculate  the  isothermal  temperature  and 
•s  time.  Figure  16  shows  an  example 
|  temperature,  on,  and  0.2  as  a  function  of 
‘  time,  with  the  isothermal  “pulse”  overlayed 
to  demonstrate  the  profile  calculated  by 
Equations  31  and  32.  Note  that  the 
isothermal  temperature  slightly  less  than 
the  peak  and  the  isothermal  time  is  about 
65%  the  FWHM  of  the  peak.  We  also 
perform  these  calculations  for  different 
peak  temperatures  and  peak  widths,  with 
the  results  tabulated  in  Table  4. 


Table  4.  Calculated  isothermal  temperature  and  time  for  different  explosion 
temperature  profiles  with  different  peak  temperatures  and  peak  widths.  Also 


Peak  T  (°K) 

FWHM  (ms) 

ai 

0t2 

Tiso  (°K) 

At  (ms) 

456.9 

25-35 

9.52E-05 

1.20E-06 

445-5 

15-1 

640.75 

25.35 

0.018114 

0.000826 

619.4 

14.5 

861 

25.35 

0.425829 

0.055041 

823.9 

15.1 

1008 

25-35 

0.90945 

0.29236 

957-7 

15-7 

1008 

36.85 

0.972596 

0.404128 

957-6 

23-5 

1008 

48.95 

0.9918 

0.499249 

957-8 

31-4 

1008 

67 

0.999254 

0.645362 

957-7 

47.1 

1008 

95 

0.99993 

0.748766 

959-0 

61.8 

From  Table  4  we  find  that  using  the  isothermal  approximation  leads  to  the  calculated 
temperature  being  slightly  less  than  the  peak  temperature  (approx.  95%)  and  the 
isothermal  time  being  approximately  65%  the  peak  width. 
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B.  Isothermal  Conversion  Calculation 

In  the  previous  sections  we  considered  different  heating  profiles  and  their  effects  on  the 
degree  of  conversion.  Additionally,  we  considered  how  the  isothermal  assumption 
compared  to  the  results  of  the  full  model  and  found  that  they  diverge  for  low  temperatures 
and  short  isothermal  times.  These  results  suggest  that  a  fast  heating  rate  is  necessary  to 
accurately  use  isothermal  temperature-time  calibration  curves.  In  this  section  we  assume 
an  essentially  infinite  heating  rate  and  consider  how  temperature  and  time  effect  the 
degree  of  conversion,  using  the  same  model  of  nucleation  as  used  above.  We  begin  by 
considering  a  single  phase  system  (similar  to  Eu:Y203),  then  a  two-phase  system  (similar 
to  Eu:Zr02),  and  finally  we  use  experimental  data  for  Eu:Zr02  to  determine  the  A  and  E 
constants  and  compute  the  nucleation  model  based  on  these  parameters. 

Single  Phase 

For  our  single  phase  model  we  assume  that  only  nucleation  and  growth  occur  with  the 
descriptive  differential  equation  given  by  Equation  20  with  m  =  o,  n  =  1,  and  p  =  o.  The 
calculations  are  performed  assuming  a  rate  parameter  of  A  =  10000  and  an  energy  barrier 
of  E  =  0.75  eV.  Figure  17  shows  a  two-dimensional  mapping  of  the  degree  of  conversion 
as  a  function  of  isothermal  time  and  temperature.  From  Figure  17  we  find  that  for  low 
temperatures  the  degree  of  conversion  is  negligible,  even  at  large  isothermal  times,  and 
that  as  the  temperature  increases  the  time  to  full  conversion  shortens.  To  better 
demonstrate  these  effects  we  plot  the  degree  of  conversion  as  a  function  of  temperature 
for  different  isothermal  times  in  Figure  18  (a)  and  the  degree  of  conversion  as  a  function 
of  isothermal  time  for  different  temperatures  in  Figure  18  (b).  The  degree  of  conversion 
is  found  to  be  exponential  in  time  and  a  sigmoid  function  in  temperature.  From 
calculations  using  different  values  for  A  and  E,  we  find  that  E  determines  the  half¬ 
conversion  temperature  of  the  sigmoid  function  and  A  is  the  primary  parameter  for  the 
rate  constant  when  considering  the  isothermal  time. 


Isothermal  Time  (ms) 

Figure  17.  Two-dimensional  mapping  of  the  degree 
of  conversion  as  a  function  of  temperature  and  time. 
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Additionally  from  Figure  17  we  find  that  there  exists  different  combinations  of 
temperature  and  time  which  produce  the  same  degree  of  conversion.  Rearranging 
Equation  23  we  find  that  the  iso-conversion  curves  to  be  defined  by  the  relation: 


ln(At) 

Tiso 


—  -  a)] 


(33) 


where  At  is  the  isothermal  time  and  Tiso  is  the  isothermal  temperature. 


Figure  18.  Degree  of  conversion  as  a  function  of  (a)  temperature  and  (b)  time. 

Two-Phases 

While  Y2O3  is  known  to  only  have  the  one  crystalline  phase,  ZrC>2  has  two  crystalline 
phases:  tetragonal  and  monoclinic.  From  experiments  we  know  that  the  Tetragonal  phase 
forms  first  and  only  after  the  tetragonal  phase  is  formed  does  the  monoclinic  phase  form. 
Additionally,  we  know  that  for  temperatures  below  approximately  400  °C  the  monoclinic 
phase  doesn’t  form.  Based  on  these  facts  and  using  the  same  nucleation  model  for  both 
phases  we  can  write  a  set  of  couple  differential  equations  describing  the  fraction  of 
material  that  is  in  the  tetragonal  phase  ar ,  and  in  the  monoclinic  phase  am : 

^  [1  _  “r®  "  “*<«]  <34) 


f  0  T  <  7\ 

| AMexp{^}aT(tl  T>Tt 

(35) 

where  Ti  is  the  minimum  temperature  required  for  the  formation  of  the  monoclinic  phase. 

With  Equations  34  and  35  we  model  the  tetragonal  and  monoclinic  fractions  as  a  function 
of  isothermal  temperature  and  time.  Figure  19  (a)  shows  a  2D  mapping  of  the  tetragonal 
fraction  as  a  function  of  temperature  and  time,  while  Figure  19  (b)  shows  a  2D  mapping 
of  the  monoclinic  fraction  as  a  function  of  temperature  and  time.  From  Figure  19  we  find 
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that  at  low  temperatures  and  times  the  tetragonal  fraction  plot  is  similar  to  the  one  phase 
model,  but  as  the  time  and/or  temperature  increases  the  tetragonal  fraction  begins  to 
decrease  as  the  monoclinic  fraction  forms,  eventually  dominating  the  sample. 

Intensity  Ratio 

While  we  calculated  the  tetragonal  and  monoclinic  fraction  above,  what  we  measure 
experimentally  is  the  fluorescence  intensity,  and  more  specifically  the  ratio  between 
peaks.  In  order  to  get  a  better  sense  of  how  the  two  phases  affect  the  spectrum  we  calculate 
an  estimated  intensity  ratio  between  the  peak  at  592  nm  and  613  nm  for  Eu:Zr02.  We 
begin  by  calculating  the  experimental  intensity  ratio  from  fully  amorphous  (rQ  =  0.22), 
tetragonal  (rr  =  0.5715),  and  monoclinic  (tm  =  0.0476)  Eu:Zr02. 
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Figure  19.  Two-dimensional  mappings  of  the  fraction  of  material  in  the  (a)  tetragonal 
phase  and  (b)  monoclinic  phase. 

Using  the  experimental  values  and  a  simplistic  model  of  the  intensity  ratios  we  can  write 
the  measured  intensity  ratio  for  a  mixture  of  the  three  phases  as 

rmix  =  (1  -aT-  aM)ra  +  aTrT  +  aMrM  (  36  ) 

With  Equation  36  and  the  fractions  computed  in  the  previous  section  we  determine  the 
intensity  ratio  as  a  function  of  temperature  and  time,  which  is  shown  in  Figure  21.  As  the 
material  becomes  more  rich  in  the  tetragonal  phase  the  ratio  between  592  nm  and  613  nm 
grows,  but  then  as  the  monoclinic  phase  begins  growing  the  ratio  rapidly  decreases  as  the 
material  becomes  rich  in  the  monoclinic  phase,  which  is  characterized  by  a  large  peak  at 
613  nm.  To  get  a  better  sense  of  this  functional  form  we  plot  the  intensity  ratio  as  a 
function  of  temperature  at  different  isothermal  times  in  Figure  22  (a),  and  as  a  function 
of  isothermal  time  at  different  temperature  in  Figure  22  (b).  From  Figure  23  we  find  that 
the  ratio  as  a  function  of  temperature  and  time  goes  through  an  inflection  point  as  the 
sample  changes  from  being  rich  in  the  tetragonal  phase  to  being  rich  in  the  monoclinic 
phase.  The  location  of  the  inflection  point  depends  on  the  isothermal  time  and 
temperature. 
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Modeling  with  Experimentally  Determined  Parameters 

Thus  far  we  have  performed  phase  change  modeling  using  arbitrary  values  which  are  on 
the  order  of  magnitude  of  those  found  in  the  literature.  At  this  point  we  turn  to  looking  at 
experimental  data  to  determine  estimates  of  the  rate  constant  and  energy  barrier  of 
Eu:Zr02.  To  begin  we  set  E=(o.131±o.oio)  eV  and  A=  8.5(±1.3)  x  lcr2  ms-1.  Figure  20 
shows  the  normalized  fluorescence  spectra  of  amorphous,  tetragonal,  and  monoclinic 
Eu:Zr02.  While  all  the  spectra  are  dominated  by  the  Eu3+  5D0  A  ?Fj  transitions,  there  are 
clear  differences  in  peak  positions  and  relative  peak  intensities.  When  amorphous 
Eu:Zr02  is  heated  up,  it  first  transforms  into  the  tetragonal  form  and  then  into  the 
monoclinic  form  for  higher  temperatures  and  longer  heating  times,  see  Figure  21. 


Figure  20.  Normalized  spectra  of 
different  phases  of  Eu:Zr02  using  an 
excitation  wavelength  of  355  nm. 
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Figure  21.  Two-dimensional  mapping  of 
the  intensity  ratio  between  592  nm  and  613 
nm  as  a  function  of  temperature  and  time. 
The  red  area  corresponds  to  being  rich  in 
the  tetragonal  phase,  while  the  blue  area  is 
rich  in  the  monoclinic  phase. 


Figure  22  shows  the  calculated  ratio  between  the  two  phases  as  a  function  of  isothermal 
temperature  for  various  heating  times,  and  as  a  function  of  heating  times  for  various 
isothermal  temperatures.  Initially,  the  ratio  increases,  indicating  the  increasing  amount 
of  the  tetragonal  phase.  However,  once  the  monoclinic  phase  starts  to  form,  the  ratio 
starts  to  decrease.  Figure  23  shows  the  experimental  intensity  ratios  as  a  function  of 
isothermal  temperature  and  isothermal  heating  time  and  fits  to  the  data. 
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Figure  22.  I592/I613  ratio  as  a  function  of  (a)  temperature  and  (b)  time.  As  the 
amorphous  material  changes  into  the  tetragonal  phase  the  ratio  increases.  Once  the 
material  begins  to  convert  into  the  monoclinic  phase  the  ratio  begins  to  decrease. 


Figure  23.  Intensity  ratios  as  a  function  of  isothermal  temperature  and  time. 
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C.  Temperature  Sensors 

Initially,  a  variety  of  potential  sensor  materials  consisting  of  individual  and  core/shell 
rare-earth-doped  metal-oxides  including  precursors  of  Eu:Y203,  Eu:Gd203,  Eu:YGd03, 
Er,Yb:Y203,  Tm,Yb:Y203,  Si02/Eu:Y203,  Eu:Zr02,  Yb203/Er:Y203,  Er203/Yb:Y203, 
Eu:Zr02/Er,Yb:Y203,  and  Tb:Zr02/Eu:Y203  were  synthesized,  characterized,  and 
evaluated.  It  quickly  became  apparent  that  Eu:Y203  and  Eu:Zr02  were  the  most  promising 
sensor  materials  as  they  showed  the  best  sensitivity  with  respect  to  temperature-induced 
changes.  As  such,  they  will  be  described  here  in  more  detail. 

Eu:Y203 


Characterization  o  fcis-siinthesized  material 

Precursor  europium-doped  yttria  (p-Eu:Y203)  is  synthesized  using  co-precipitation, 
resulting  in  a  white  powder.  The  use  of  mixed  precipitants,  ammonium  hydroxide  and 
ammonium  bicarbonate,  leads  to  the  formation  of  nanosized  morphologies,  see  Figure  24 
(a).  Since  the  nanoparticles  tend  to  aggregate,  it  is  not  possible  to  accurately  measure  the 
particle  size.  The  particles  have  irregular  shapes  with  dimensions  in  the  20  -  50  nm  range. 
Nanoparticles  prepared  in  this  manner  are  known  to  favor  the  formation  of  de- 
agglomerable  (soft  agglomerate)  oxide  powders  upon  calcination  6.  X-ray  diffraction 
measurements  show  that  the  as-prepared  precipitate  is  amorphous,  as  seen  from  the 
broad  halo  curve  in  Figure  24  (b).  It  is  worth  noting  that  at  a  pH  greater  than  9.0,  the 
nanosized  particles  are  known  to  form  two-dimensional  microstructures.  However,  this 
formation  was  avoided  in  this  work  by  keeping  the  pH  at  a  lower  value.  According  to 
Huang  et  al.  7,  at  a  lower  pH  the  precipitate  consists  of  Y2(C03)3-nH20,  while  at  higher  pH 
the  Y2(C03)3-nH20  species  transform  into  Y(0H)C03.  The  hydroxyl  group  in  Y(0H)C03 
is  capable  of  forming  oxo-bridges  that  lead  to  the  formation  of  micron-sized  aggregates. 
As  they  spontaneously  assemble  into  micro-aggregates,  the  yttrium  hydroxide 
nanoparticles  tend  to  assemble  more  rapidly  in  the  direction  perpendicular  to  the 
resulting  plane  of  two-dimensional  micro-structures.  Upon  calcination,  Y(0H)C03 
powders  form  hard  agglomerate  7.  However,  because  of  the  lower  pH  used,  the  likelihood 
for  the  presence  of  Y(0H)C03  in  our  samples  is  low. 


Wavenumber  (cm'1) 


Figure  24.  (a)  TEM  image  (scale  bar  is  100  nm),  (b)  XRD  spectrum,  and  (c)  FT-IR 
spectrum  of  coprecipitated  yttrium  carbonate. 


Figure  24  (c)  shows  the  FT-IR  spectrum  of  the  amorphous  precursor.  The  broad 
absorption  peak  at  3,440  cm1  is  due  to  the  O-H  stretching  mode,  indicating  the  presence 
of  hydroxide  ions.  The  absorption  peaks  at  1,520  cm1  and  1,385  cm1  are  due  to  the 
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asymmetric  stretching  of  C-0,  and  the  absorption  peaks  at  1,084  cnr1  and  845  cm-1  are 
due  to  the  symmetric  and  deformation  vibration  of  C-O,  respectively,  indicating  the 
presence  of  carbonate  ions.  The  bands  between  800  cnr1  and  650  cnr1  are  due  to  the 
deformation  of  C032'  and  O-Hv 

Characteristics  of  the  heat-treated  uttria  precursor 

Between  room  temperature  and  700°C,  the  TGA  curve  shows  a  continuous  mass  loss  of 
up  to  about  33  wt%,  see  Figure  25.  This  large  mass  loss  indicates  that  the  yttria  precursor 

is  indeed  a  carbonate  rather  than  a  pure 
hydroxide.8  The  DSC  curve  shows  an 
exothermic  transition  that  has  been  attributed 
to  the  formation  of  Y2O3  following  the 
decomposition  of  yttrium  carbonate  at 
665.7°C  9.  Above  this  temperature,  the  rate  of 
mass  loss  is  significantly  reduced,  indicating 
an  almost  completed  formation  of  the  Y2O3 
phase.  Due  to  prior  drying  (see  section  2),  no 
apparent  endothermic  transitions  in  the  100- 
150°C  range,  corresponding  to  loss  of  water, 
were  observed.  The  overall  decomposition 
steps  for  yttrium  carbonate  precursors  have 
been  reported  previously  ?.  The 
decomposition  begins  with  the  loss  of  water, 
physisorbed,  chemisorbed,  and  molecular  water,  followed  by  the  release  of  C02  resulting 
in  Y20(C03)2.  Due  to  further  release  of  C02,  Y20(C03)2  decomposes  into  Y202C03,  which 
finally  transforms  into  Y203  under  loss  of  C02.  As  such,  the  main  expected  decomposition 
processes  (under  release  of  C02)  is:  Y2(C03)3  — ►  Y20(C03)2  — ►  Y202C03  — ►  Y203. 

The  initial  formation  of  Y203  is  corroborated  by  the  appearance  of  an  FT-IR  absorption 
peak  at  570  cm-1,  which  is  indicative  of  Y-0  bonds,  see  Figure  26  (a)-(e).  The  earliest 
appearances  of  the  Y-0  bond  can  be  seen  for  heating  to  500°C  for  5  minutes,  see  Figure  3 
(b).  The  continuous  loss  of  water  and  decomposition  towards  the  formation  of  Y203  is 
apparent  by  the  relative  increase  in  intensity  of  the  Y-0  peak  at  570  cnr1  compared  to  the 
O-H  and  C-0  peaks.  However,  some  remnants  of  the  peaks  corresponding  to  C-0  and  O- 
H  bonds  are  still  visible  even  after  furnace  heating  to  noo°C  for  180  minutes.  These  peaks 
are  due  to  surface  hydroxyl  and  carbonate  groups  that  form  when  Y203  reacts  with 
ambient  water  and  carbon  dioxide,  respectively 10. 


Temperature  (°C) 

Figure  25.  DSC  and  TGA  curves  of  the 
co-precipitated  carbonates. 
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Figure  26.  FT-IR  spectra  of  co-precipitated  yttria  precursors  heat-treated  at  300°C, 
500°C,  700°C,  900°C,  and  noo°C  for  dwell  times  of  5  minutes,  30  minutes,  and  180 
minutes. 


Figure  27  shows  XRD  spectra  of  yttria  precursors  heat-treated  at  various  temperatures 
and  dwell  times.  Figure  27  (a)  shows  that  at  300°C  the  matrix  remains  amorphous  even 
for  dwell  times  of  180  minutes.  The  influence  of  dwell  time  on  the  crystal  quality  is  best 
demonstrated  in  Figure  27  (b).  As  the  dwell  time  increases  at  500°C,  the  peaks  become 
narrower,  indicating  the  formation  of  crystalline  cubic  Y2O3.  The  three  highest  intensity 
x-ray  diffraction  peaks,  matching  the  crystalline  cubic  Y2O3  phase  (JCPDS#  98-010-7439) 
at  20  =  29. 1°,  48.5°,  and  57.6°,  start  to  appear  for  the  sample  heated  to  500°C  for  30 
minutes.  The  presence  of  the  cubic  crystalline  phase  of  Y2O3  is  apparent  in  the  spectra  of 
samples  heated  to  700°C  and  above,  see  Figure  27  (c)-(e).  Since  Y2O2CO3  is  thermally 
stable  only  up  to  470°C,  the  crystalline  peaks  seen  for  heating  to  500°C  for  30  minutes 
are  likely  to  correspond  to  Y2O3.11  After  a  dwell  time  of  180  minutes  the  peaks  can  be 
unambiguously  matched  to  the  cubic  Y2O3  phase. 
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Figure  27.  XRD  patterns  of  co-precipitated  yttria  precursors  heat-treated  at  300°C, 
500°C,  700°C,  900°C,  and  noo°C  for  dwell  times  of  5  minutes,  30  minutes,  and  180 
minutes. 


Figure  28  shows  the  TEM  images  of  the  heated  samples.  The  images  show  that  the 
samples  heated  to  500°C  have  a  perforated  appearance,  which  may  be  attributed  to  the 
carbonate  decomposition  resulting  in  the  release  of  carbon  dioxide.  This  decomposition 
results  in  smaller  particle  sizes  for  samples  heated  to  700°C  for  dwell  times  of  5  minutes, 
30  minutes,  and  180  minutes.  Above  700°C,  these  crystalline  particles  grow  in  size  with 
increasing  temperature  —  up  to  70  -  80  nm  in  diameter  for  heating  to  noo°C  for  180 
minutes.  The  particle  size  of  this  sample,  as  determined  from  the  TEM  image  is 
comparable  to  the  grain  size  determined  from  the  XRD  spectrum,  indicating  complete 
crystallization.  These  particles  have  polyhedral  geometry  and  are  loosely  agglomerated. 
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Figure  28.  TEM  images  of  co-precipitated  carbonates,  heat-treated  at  various 


temperatures  and  dwell  times.  Scale  bar  is  100  nm. 


Figure  29  (a)  shows  the  grain  size  (determined  from  XRD  data)  vs.  dwell  time  for  various 
heating  temperatures.  The  trend  in  the  crystalline  growth  is  characteristic  for  nanosized 
particles.12  Initially,  the  grains  grow  rapidly  within  a  short  period  of  annealing  time.  This 
period  is  followed  by  a  period  of  substantially  slower  growth  rate,  and  finally  a  period  in 
which  the  growth  levels  off.  During  the  last  stage,  the  grain  size  approaches  its  limiting 
grain  size  for  a  given  temperature.  No  crystalline  domain  is  observed  at  300°C.  At  all 
higher  temperatures,  the  grain  growth  slows  down  considerably  after  about  two  hours. 
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Figure  29.  (a)  Grain  size  vs.  dwell  times  for  yttrium  carbonates,  heat-treated  at  various 
temperatures  and  dwell  times,  (b)  Plot  of  Ln(K)  versus  l/T  (K-i). 

Fitting  of  the  above  data  using  the  conventional  grain  growth  equation:12 


Dn  —  Dq  =  Kt  (37) 

with  D  being  the  average  grain  size  after  annealing,  Do  being  the  initial  grain  size,  n  being 
the  grain  growth  exponent,  t  being  the  annealing  time,  and  K  being  a  constant 
proportional  to  exp(-Q/RT)  with  Q  being  the  activation  energy,  R  being  the  gas  constant, 
and  T  being  the  calcination  temperature,  yielded  different  values  for  the  grain  growth 
exponent  n  for  the  various  calcinations  temperatures. 


Fitting  our  data  to  a  structural  relaxation  equation  yielded  more  consistent  results  12. 
While  all  data  points  for  the  900  and  noo°C  isothermals  could  be  fitted  to  this  equation, 
for  the  500  and  700°C  isothermals  the  data  points  for  dwell  times  of  5  minutes  and  30 
minutes  could  not  be  fitted  and  were  left  out.  It  may  be  that  a  different  process  occurs  for 
the  shorter  dwell  times  at  lower  temperatures.  The  relaxation  equation: 


P-Pp  _  !  _  e-Ktm 
Dm~D  0 


(38) 


with  D  being  the  grain  size,  Do  being  the  initial  grain  size,  Dm  being  the  limiting  grain  size, 
t  being  the  annealing  time,  m  being  the  relaxation  order  (=1  for  linear  relaxation),  and  K 
being  the  rate  constant  proportional  to  exp(-Q/RT)  with  Q  being  the  activation  energy,  R 
being  the  gas  constant,  and  T  being  the  calcination  temperature,  describes  crystalline 
growth  due  to  gradual  ordering  of  the  bond  lengths  and  bond  angles  at  the  nanoparticles’ 
interface  12.  This  model  has  been  successfully  used  to  describe  isothermal  growth  of  sol- 
gel  synthesized  nanoparticles.  Figure  29  (b)  shows  the  plot  of  the  log  of  the  rate  constant 
vs.  inverse  temperature.  Fitting  of  this  data  set  provides  the  activation  energy  for  the  grain 
growth  process.  Table  5  shows  the  fitted  parameters  for  each  temperature.  The  fitted 
activation  energy  of  50.3  KJ/mol  is  on  par  with  values  reported  in  the  literature  for 
nanosized  particles.12 
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Table  5.  Fitting  parameters  for  relaxation  equation. 


Temperature  (°C) 

K 

Dm 

m 

x2 

500 

0.017 

14.1 

0.90 

0.009 

700 

0.099 

20.0 

0.55 

0.075 

900 

0.280 

46.0 

0.28 

0.033 

1100 

0.509 

98.0 

0.18 

4-32 

Optical  properties 

Eu3+  dopants  can  be  used  to  monitor  the  structural  properties  of  their  host  and  to  gain 
information  about  how  these  properties  change  in  response  to  stimuli.^  For  example,  if 
the  europium  ion  is  located  in  a  site  with  inversion  symmetry,  electric  dipole  transitions 
are  forbidden.1^  *5  As  the  site  symmetry  is  lowered,  electric  dipole  transitions  become 
allowed  and  the  optical  spectra  change  from  being  dominated  by  magnetic  dipole 
transitions  to  being  dominated  by  electric  dipole  transitions.  Also,  when  the  host  is 
amorphous,  a  variety  of  sites  are  available  for  the  europium  ions,  leading  to 
inhomogeneous  broadening  of  the  optical  spectra.  As  the  material  changes  from  an 
amorphous  phase  to  a  crystalline  phase,  the  spectra  become  narrower,  indicating  the 
more  homogeneous  nature  of  the  host  material. 


Excitation  Wavelength  (nm) 

Figure  30.  Excitation  spectra  of  Eu-doped  samples  calcined  at  various  temperatures 
for  5  minutes  (top),  30  minutes  (middle),  and  180  minutes  (bottom),  observed  at  611 
nm. 

Figure  30  shows  the  excitation  spectra  of  europium-doped  nanoparticles,  monitored  at 
611  nm  and  heated  to  various  temperatures  (300°C,  500°C,  700°C,  900°C,  and  noo°C) 
for  dwell  times  of  5  minutes,  30  minutes,  and  180  minutes.  The  excitation  spectra  of  the 
samples  heated  to  300°C  for  all  dwell  times  show  a  broad  peak  centered  at  about  578.3 
nm.  As  the  temperature  increases,  this  broad  peak  decreases  rapidly  and  instead,  a  sharp 
peak  centered  at  about  580.8  nm  (?Fo— >5Do)  and  some  weak  double  peaks  at  about  568 
nm  and  570  nm  appear.  The  broad  excitation  peak  centered  at  about  578.5  nm  disappears 
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with  increasing  temperature.  This  disappearance  is  due  to  the  decomposition  of  the 
amorphous  yttrium  carbonate  phase.  As  seen  in  Figure  30,  the  full-width-half-max 
(FWHM)  of  the  narrow  peak  at  580.8  nm,  representing  the  formation  of  the  crystalline 
Y2O3  phase,  continues  to  decrease  with  increasing  temperature.16  No  obvious  distinctions 
are  apparent  in  the  excitation  spectra  of  samples  heated  for  different  dwell  times. 


Figure  31.  Fluorescence  spectra  of  Eu-doped  samples  calcined  at  various  temperatures 
for  5  minutes,  30  minutes,  and  180  minutes,  and  excited  within  excitation-peaks 
observed  in  Figure  30. 

The  fluorescence  spectra  of  the  samples  heated  to  300°C,  500°C,  700°C,  900°C,  and 
noo°C  for  dwell  times  of  5  minutes,  30  minutes,  and  180  minutes  are  shown  in  Figure 
31.  The  samples  heated  to  300°C  were  excited  at  577  nm,  in  the  broad  excitation  peak  for 
the  amorphous  phase.  This  excitation  resulted  in  broad  emission  lines  centered  at  about 
620  nm.  The  samples  heated  to  500°C,  700°C,  900°C,  and  iioo°C  were  excited  at  580.8 
nm,  leading  to  emission  centered  at  about  611  nm,  typical  for  crystalline  Y203:Eu.  The 
peaks  between  585  nm  and  605  nm  are  due  to  the  sDo— >7Fi  transition,  the  peaks  between 
610  nm  and  635  nm  are  due  to  the  sDo— >7F2  transition,  and  the  peaks  around  650  nm  are 
due  to  the  5D0— >7F3  transition.17  While  the  FWHM  of  the  excitation  peak  of  the  crystalline 
Y203  phase  continues  to  decrease  with  temperature,  the  FWHM  of  the  HDo— >7F2 
fluorescence  of  the  crystalline  Y203  phase  does  not  show  any  apparent  decrease.  This 
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difference  is  due  to  the  fact  that  for  the  excitation  spectra,  monitored  at  611  nm,  the 
amorphous  phase  and  the  crystalline  phase  both  emit,  but  the  amorphous  phase 
disappears  with  increasing  temperature,  reducing  the  FWHM.  For  the  fluorescence 
spectra  however,  580.8  nm  excites  only  the  crystalline  phase  and  not  the  amorphous 
phase.  Similar  to  the  excitation  spectra,  no  obvious  distinctions  are  apparent  in  the 
fluorescence  spectra  of  samples  heated  for  different  dwell  times. 


Temperature  (°C) 

Figure  32.  Fluorescence  lifetime  of  Eu-doped  samples, 

heat  treated  at  various  temperatures  and  dwell  times. 

Figure  32  shows  the  fluorescence  lifetimes  of  the  various  europium-doped  samples.  In 
each  case,  a  single  exponential  was  used  to  fit  the  data.  The  traces  for  the  three  furnace- 
heated  samples  (5  minutes,  30  minutes,  and  180  minutes)  are  very  close  to  each  other. 
Samples  heated  to  low  temperatures  start  out  at  about  1-1.2  ms.  As  the  temperature  of  the 
samples  is  increased,  the  fluorescence  lifetime  increases  and  reaches  a  maximum  of  about 
2.4  ms  at  700°C.  Samples  heated  to  higher  temperatures  show  a  decrease  in  fluorescence 
lifetime.  As  already  seen  before,  the  samples  undergo  a  phase  transition  at  around  665°C 
and  convert  into  a  crystalline  cubic  oxide  phase.  As  more  and  more  of  the  sample  material 
undergoes  this  transition,  the  fluorescence  lifetime  settles  at  a  value  of  about  1.2  ms, 
which  is  typical  for  cubic  Eu:Y2C>3.  The  fluorescence  lifetimes  observed  at  lower 
temperatures  reflect  the  various  decomposition  products  present  in  the  material.  The 
dwell  time  at  each  temperature  setting  does  not  appear  to  affect  the  fluorescence  lifetime 
very  much. 

Puroprobe-heated  YoO  r.Eu 

In  the  lab,  samples  were  heated  via  a  CDS  Analytical  Pyroprobe  1000.  The  nanopowders 
were  placed  in  a  quartz  tube,  which  was  sealed  at  both  ends.  Subsequently,  the  quartz 
tube  was  placed  inside  a  coil  heater  and  heated  to  various  temperatures  and  dwell  times 
using  a  ramp  rate  of  20000°C/s.  However,  these  setting  are  based  upon  heating  of  the 
coil  heater,  not  the  material  inside  the  quartz  tube.  Reference  measurements  were 
performed  using  a  thermocouple  placed  inside  the  quartz  tube.  Figure  33  shows  the 
results. 
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As  can  be  seen  from  the  data  in  Figure  10,  for  short  heating  times  of  l  s,  the  maximum 
temperature  reached  was  less  than  50%  of  the  set  temperature,  while  for  longer  heating 
times,  the  maximum  temperature  is  10-25%  below  the  set  temperature.  Also,  the  observed 
ramp  rate  is  significantly  slower  that  the  set  ramp  rate.  For  the  further  analysis,  the 
measured  temperatures  were  used. 


Figure  33.  Reference  measurements  of  the  pyroprobe  heater. 


Figure  34  shows  a  TEM  image  of  the  co-precipitated  and  dried  precursor.  The  particles 
are  about  20  -  50  nm  in  diameter  and  appear  to  be  agglomerated.  Calcination  of  this 
precursor  is  known  to  lead  to  the  formation  of  Y203:Eu  powders  characterized  by  soft 
agglomeration.6  Figure  35  shows  the  DSC  and  TGA  curves  of  the  precursor.  Between  room 
temperature  and  700°C,  the  TGA  curve  shows  a  continuous  mass  loss  of  up  to  about  33 
wt%.  This  large  mass  loss  indicates  that  the  yttria  precursor  is  indeed  a  carbonate  rather 
than  a  pure  hydroxide.8  According  to  Zhuang  et  al., 7  the  precipitate  is  expected  to  consists 
ofY2(C03)3-nH20.  The  DSC  curve  shows  an  exothermic  transition  at  665. 7°C,  which  has 
been  attributed  to  the  formation  of  crystalline  Y2O3  following  the  decomposition  of 
yttrium  carbonate. 9  Above  this  temperature,  the  rate  of  mass  loss  is  significantly  reduced, 
indicating  an  almost  completed  formation  of  the  Y203  phase.  Because  of  prior  drying  (see 
section  2),  no  apparent  endothermic  transitions  in  the  100-150°C  range,  corresponding 
to  loss  of  water,  are  observed. 
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Figure  34.  TEM  image  of  co-precipitated 
precursor.  Scale  bar  =  50  nm. 


Temp  (°C) 

Figure  35.  DSC  and  TGA  curves  of  the 
co-precipitated  carbonates. 


The  overall  decomposition  steps  have 
been  reported  for  yttrium  carbonate 
precursors.?  The  decomposition  begins 
with  the  loss  of  water,  physisorbed, 
chemisorbed,  and  molecular  water, 
followed  by  the  release  of  CO2  resulting 
in  ¥20(003)2.  Due  to  further  release  of 
CO2,  Y20(C03)2  decomposes  into 
Y2O2CO3,  which  finally  transforms  into 
Y2O3  under  further  loss  of  C02.  As  such, 
the  main  expected  decomposition 
process  (under  release  of  CO2)  is: 
Y2(C03)3  -►  Y20(C03)2  Y2O2CO3 
Y2O3.  According  to  Hussein,  Y2(C03)3 
transforms  into  Y20(C03)2  and  Y2O2CO3 
at  375°C  and  400°C,  respectively.18 

Figure  36  shows  the  FT-IR  spectra  of  the 
precursor  and  samples  heated  to  various 
temperatures  for  dwell  times  of  1  s,  50  s, 
and  100  s.  The  precursor  spectrum, 
shown  in  Figure  36  (a),  is  consistent  with 
the  presence  of  hydroxyl  and  carbonate 
groups.  The  O-H  stretching  mode  is 
responsible  for  the  broad  peak  at  3440 
cnr1  and  the  asymmetric  stretching  of  C- 
O  is  responsible  for  the  peaks  at  1520  cnr 
1  and  1385  cnr1.  The  C-0  symmetric  and 
deformation  vibrations  absorb  at  1084 
cm-1  and  845  cm-1.  Absorption  peaks 
between  800  cnr1  and  650  cnr1  are  due  to 
deformation  vibrations  of  CO32-  and  OH\ 


Figure  36  (b)-(f)  show  the  FT-IR  spectra  of  the  precursor  after  it  has  been  heated  for  dwell 
times  of  1  s,  50  s,  and  100  s  to  various  temperatures.  All  spectra  indicate  the  presence  of 
hydroxyl  and  carbonate  groups,  which  form  at  the  surface  when  the  material  reacts  with 
ambient  water  and  carbon  dioxide.  The  presence  of  Y-0  (570  cm-1)  shows  in  the  spectra 
of  samples  heated  for  50  s  and  100  s  to  720°C  and  850°C.  This  observation  is  consistent 
with  the  DSC  measurement  indicating  the  formation  of  Y2O3  at  about  665°C. 
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Figure  36.  FT-IR  spectra  of  (a)  precursor  and  (b)-(f)  precursor  heated  to  various 
temperatures  for  1  s,  50  s,  and  100  s. 

Figure  37  shows  the  XRD  spectra  of  the  precursor  and  the  heated  precursors.  A  crystalline 
signature  first  appears  for  the  sample  heated  to  640°C  for  50  s.  The  diffraction  peaks  are 
consistent  with  cubic  Y2O3  (JCPDS  No.  98-010-7439).  Table  6  summarizes  the  grain  size 
of  the  samples  heated  to  640°C,  720°C,  and  850°C  for  dwell  times  of  50  s  and  100  s.  This 
data  indicates  only  small  changes  in  grain  size  heating  temperature  and  dwell  time. 
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Figure  37.  XRD  of  (a)  precursor  and  (b)-(f)  precursor  heated  to  various  temperatures 
for  1  s,  50  s,  and  100  s. 
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Table  6.  Grain  size  of  samples  heated  to  640°C,  720°C,  and  850°C  for  dwell  times  of 


50  s  and  100  s. 


Time  (s) 

Temperature  (°C) 

640 

720 

850 

50 

11.5  nm 

14.1  nm 

14.5  nm 

too 

13.0  nm 

14.5  nm 

18.2  nm 

While  the  grain  growth  kinetics  are  often  described  by  the  conventional  grain  growth 
equation  or  the  relaxation  equation,12  these  equations  are  unable  to  describe  any  growth 
occurring  at  the  much  shorter  time  scales  described  in  this  work.  This  inability  is  due  to 
the  multi-step  decomposition  process  that  takes  place  as  described  in  Equation  33. 

Figure  38  shows  the  TEM 
images  of  samples  heated  for 
dwell  times  of  1  s,  50  s,  and  100 
s  at  various  temperatures.  The 
samples  heated  to  640  °C  for  50 
s  and  100  s  have  a  perforated 
appearance,  indicating 

decomposition  and  release  of 
CO2  gas.  The  TEM  images  of 
the  samples  heated  to  720°C 
for  100  s  and  850°C  for  50  s 
and  too  s  show  the  crystalline 
lattice.  Due  to  their  tendency  to 
agglomerate,  measurement  of 
the  particle  size  from  the  TEM 
images  is  not  possible.  Also, 
because  the  material  has  mixed 
amorphous  and  crystalline 
phases,  the  crystallite  size 
derived  from  XRD  spectra  does 
not  necessarily  correspond  to 
the  particle  size  seen  in  the 
TEM  images. 

Figure  39  shows  the  excitation 
spectra  and  fluorescence 
spectra  of  the  two  extreme 
cases,  amorphous  Y2(C03)3:Eu  Figure  38.  TEM  images  of  the  co-precipitated 
precursor  and  fully  crystalline  carbonates  heat  treated  at  various  temperatures  and 
Y203:Eu.  The  excitation  peaks  dwell  times.  Scale  bar  is  50  nm. 
are  due  to  the  pure  electronic 

7F0— >5D0  transition  in  the  amorphous  and  crystalline  samples,  respectively.^  The 
fluorescence  peaks  are  due  to  the  sDo— >7Fi  (585  nm  -  600  nm)  and  5Do— >7F2  (610  nm  - 
630  nm)  transitions,  respectively. ^  The  excitation  peaks  for  the  amorphous  and 
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crystalline  samples  are  well  separated,  while  the  fluorescence  peaks  for  the  amorphous 
and  crystalline  samples  have  significant  overlap. 


Wavelength  (nm)  Wavelength  (nm) 

Figure  39.  Excitation  spectra  (left)  and  fluorescence  spectra  (right)  of  amorphous  and 
fully  crystalline  samples  (left).  Only  minimal  overlap  occurs  for  the  excitation  spectra, 
while  significant  overlap  occurs  for  the  fluorescence  spectra  of  these  two  samples. 

Under  the  influence  of  temperature  and  time,  the  broad  excitation  peak  first  shifts  toward 
shorter  wavelength  while  increasing  in  width.  As  high  enough  temperatures  are  reached 
for  extended  dwell  times,  the  crystalline  peak  appears  and  increases  in  strength,  while  the 
amorphous  peak  collapses  and  disappears.  The  crystalline  peak  width  continues  to 
decrease  with  increased  heating  temperature.16 
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Figure  40.  Excitation  spectra  (left),  monitored  at  611  nm,  of  yttrium  carbonate 
precursor  heated  for  1  s  to  various  temperatures.  Peak  wavelength  and  peak  width  as  a 
function  of  heating  temperature,  (right)  based  on  Gaussian  fit. 


Figure  40  shows  the  excitation  spectra,  monitored  at  611  nm,  of  yttrium  carbonate 
precursors  heated  for  1  s  to  130°C,  250°C,  340°C,  38o°C,  and  440°C.  The  spectra  show 
one  broad  absorption  peak  that  shifts  to  shorter  wavelength  and  increases  in  width  as  the 
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temperature  is  increased.  A  single  Gaussian  peak  was  fitted  to  the  excitation  peak,  and  its 
peak  position  and  peak  width  are  also  shown  in  Figure  40. 

As  the  original  amorphous  Y2(C03)3  decomposes  with  increasing  temperature  and  starts 
forming  Y20(C03)2  and  Y202C03,  each  of  these  compounds  contributes  to  the  excitation 
spectra,  increasing  the  width  of  the  peak.  The  decomposition  products  Y20(C03)2  and 
Y202C03  must  have  higher  crystal  field  strength  than  Y2(C03)3  to  account  for  the  shift  to 
higher  energy  (shorter  wavelength). 


Figure  41  shows  the  excitation  spectra  and  shifts  when  monitored  at  620  nm  instead  of 
611  nm.  In  this  case,  the  peak  shift  is  about  35%  larger  than  when  monitored  at  611  nm. 
However,  the  change  in  peak  width  is  about  the  same. 
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Figure  41.  Excitation  spectra  (left),  monitored  at  620  nm,  of  yttrium  carbonate 
precursor  heated  for  1  s  to  various  temperatures.  Peak  wavelength  and  peak  width  as  a 
function  of  heating  temperature,  (right)  based  on  Gaussian  fit. 

The  larger  peak  shift  is  probably  due  to  the  overlap  of  crystalline  and  amorphous 
fluorescence  components  at  611  nm  while  monitoring  almost  pure  amorphous 
components  at  620  nm.  The  crystalline  peak  component  is  not  expected  to  shift  as  much 
as  the  amorphous  peak  components. 

Figure  42  shows  how  the  excitation  spectra  of  the  pyroprobe-heated  yttrium  carbonate 
precursor  evolve  as  a  function  of  temperature  and  time  while  monitored  at  611  nm  and 
620  nm.  The  general  trend  starts  out  with  a  broad  peak  caused  by  the  amorphous  phase, 
which  increases  in  width  with  increasing  temperature  until  its  intensity  decreases  and  a 
narrow  peak  due  to  a  crystalline  phase  appears.  The  transition  from  amorphous  to 
crystalline  phase  occurs  earlier  in  time  and  temperature  for  monitoring  at  611  nm  than 
monitoring  at  620  nm.  This  difference  is  because  at  620  nm  the  observed  signal  is  almost 
exclusively  due  to  the  amorphous  phase,  while  at  611  nm  the  observed  signal  is  due  the 
combination  of  amorphous  and  crystalline  phase. 
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Figure  42.  Excitation  spectra  of  yttrium  carbonate  precursor  heated  via  pyroprobe  for 
dwell  times  of  25  s,  50  s,  75  s,  and  100  s  to  various  temperatures. 
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Figure  43  summarizes  the  shift  in  peak 
wavelength  and  peak  width  of  the  spectra 
monitored  at  620  nm.  The  initial  shift  in 
peak  position  with  increasing 
temperature  is  toward  shorter 
wavelength,  followed  by  a  shift  toward 
longer  wavelength.  The  shift  in  peak 
width  is  increasing  with  increasing 
temperature.  Dwell-time  does  not 
appear  to  have  much  of  an  influence  on 
the  observed  shift. 

The  data  shown  in  Figure  42  for  a  dwell-  l7.  ci  •<+  •  Teniperature(  C) 

£  •  7  -xt  Figure  43.  Shift  m  peak  position  (left  scale) 

time  ot  1  s,  are  consistent  with  the  *5  7  n  .  c  .• 

’  ,  and  peak  width  (right  scale)  of  excitation 

measurements  shown  in  Figure  43.  .  ..  \  .  ^ _  r  ..  ■ 

j  ^  ^  ^  t  ^  j  spectra  monitored  at  620  nm  ot  yttrium 

These  data  represent  the  shut  toward  ,  .  ,  . 

-i  ^  -i  ..  ^  carbonate  precursors  heated  to  various 

shorter  wavelength  for  temperatures  up  .  .  r  n  nn .  •  f  „ 

^  °  0  -T  temperatures  for  dwell  times  of  25  s,  50  s,  75 

to  about  440°C,  as  seen  in  Figure  18  with  u  0  ' 0 

higher  resolution.  The  minimum  peak  ’ 

position,  corresponding  to  the  highest  energy,  occurs  at  a  temperature  when  most  of  the 
sample  is  expected  to  be  made  up  of  Y2O2CO3.  As  such,  this  material  appears  to  have  the 
highest  crystal  field  strength  of  the  various  materials  present. 
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Figure  44.  Fluorescence  spectra,  excited  at  579  nm  (left)  and  at  580.8  nm  (right),  of 
yttrium  carbonate  pre-cursors  heated  to  various  temperatures  for  1  s. 

Figure  44  shows  the  fluorescence  spectra  of  yttrium  carbonate  precursors  heated  for  1  s 
to  various  temperatures  for  excitation  at  579  nm  and  580.8  nm,  respectively.  As  the 
heating  temperature  increases,  the  emission  peak  between  about  610  nm  and  625  nm 
increases  in  width.  When  excited  at  579  nm,  the  fluorescence  appears  to  have  two  peaks, 
centered  at  about  613  nm  and  617  nm.  With  increasing  temperature,  the  intensity  of  the 
former  decreases,  while  the  intensity  of  the  latter  increases.  Also,  three  peaks,  which  are 
due  to  the  sDo— >7Fi  transition  are  located  between  585  nm  and  600  nm.  The  fact  that  the 
7Fi  level  splits  into  three  levels  indicates  the  low  symmetry  of  the  Eu  site.  When  excited  at 
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580.8  nm,  a  main  peak  is  located  at  about  615  nm.  Only  when  the  temperature  reaches 
440°C  does  the  intensity  at  longer  wavelength  increase.  As  opposed  to  excitation  at  579 
nm,  the  spectra  for  excitation  at  580.8  nm  show  one  main  peak  with  two  shoulders 
between  585  nm  and  600  nm. 
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Figure  45.  Fluorescence  spectra,  excited  at  579  nm  (left)  and  at  580.8  nm  (right),  of 
yttrium  carbonate  pre-cursors  heated  to  various  temperatures  for  too  s. 


The  fluorescence  spectra  for  heating  of  yttrium  carbonate  precursors  for  100  s  to  various 

temperatures  are  shown  for  excitation  at  579 
nm  and  at  580.8  nm  in  Figure  45.  In  this  case, 
the  heating  temperatures  are  higher  than  for 
the  spectra  shown  in  Figure  11.  Heating  to 
higher  temperatures  highlights  the  difference 
between  the  samples  excited  at  579  nm  and 
those  excited  at  580.8  nm.  The  higher 
temperatures  lead  to  the  formation  of 
crystalline  Y2O3,  which  is  characterized  by  a 
strong  emission  peak  centered  at  about  611  nm 
and  which  has  its  main  excitation  at  580.8  nm. 
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Figure  46.  Fluorescence  spectrum  of 
yttrium  carbonate  pre-cursors  heated 
for  1  s  to  440°C  and  its  deconvoluted 
spectral  fit. 


Figure  46  shows  the  fluorescence  spectrum  of 
the  sample  heated  for  1  s  to  440°C  and  its 
deconvoluted  spectral  fit.  Because  of  the  low 
symmetry  indicated  by  the  three  peaks  for  the 
5Do— >7Fi  transition,  five  peaks  are  fitted  to  the 
5Do^7F2  transition. 


The  fluorescence  lifetimes  of  the  heated  samples  are  shown  in  Figure  47  for  excitation  at 
532  nm  and  monitoring  at  611  nm,  and  for  excitation  at  532  nm  and  monitoring  at  620 
nm.  The  samples  monitored  at  611  nm  reach  fluorescence  lifetimes  of  about  2  ms  when 
heated  to  about  850°C,  while  reaching  only  about  1.4  ms  when  heated  to  850°C  but 
monitored  at  620  nm.  This  difference  reflects  the  presence  of  crystalline  vs.  amorphous 
phase  for  the  respective  observation  wavelengths.  At  611  nm,  a  large  contribution  of  the 
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signal  is  due  to  the  crystalline  phase,  while  at  620  nm  a  large  contribution  is  due  to  the 
amorphous  phase.  As  expected,  the  crystalline  phase  has  a  longer  fluorescence  lifetime 
than  the  amorphous  phase. 


Temperature  (°C)  Temperature  (°C) 

Figure  47.  Fluorescence  lifetimes,  excited  at  532  nm  and  observed  at  611  nm 
(left)  and  620  nm  (right),  of  yttrium  carbonate  precursors. 


Explosion  tests 


The  explosion  tests  were  performed  in  a  450  cm  x  450  cm  steel  plated  blast  proof  chamber 
at  the  Sigmund  Jacobs  Detonation  Science  Facility  located  at  the  Indian  Head  Division  of 
the  Naval  Surface  Warfare  Center  in  Indian  Head,  MD,  by  Dr.  James  Lightstone  and  his 
colleagues,  who  provided  the  following  information.  Figure  48  shows  the  experimental 
setup.  The  sensor  particles,  placed  in  small  pieces  of  paper  tissue,  were  hung  -21.5  cm 
from  the  center  of  a  charge  containing  aluminum  and  HMX.  A  plastic  Imhoff  cone  with  a 
10  cm  maximum  diameter  was  used  to  collect  the  particles.  An  R-type  fast-thermocouple 
(Omega  P13R-002),  housed  in  a  small  metal  cylinder  and  placed  at  approximately  the 
same  distance  from  the  center  of  the  charge  as  the  particle  sample,  was  used  to  record  the 
temperature. 


Sample  and  Thermocouple 
located  ~8.5”  from  charge 
center. 

Charge  (C):  1”  right  cylinder 
(~20g)  of  plastic  bonded 
explosive  formulation 
containing  HMX  and  Al 
particles. 


Image  of  typical 
aluminized 
detonating 


Figure  48.  Setup  for  explosion  testing 
thermocouple  (right). 


Time  (ms) 

(left)  and  temperature  measured  via 
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The  fluorescence  and  fluorescence-excitation  spectra  of  amorphous  and  fully-crystalline 
Eu:Y2C>3  are  shown  in  Figure  49. 


Figure  49.  The  as-prepared  sample  is  amorphous  and  displays  broad  spectral  features, 
while  the  heated  sample  is  crystalline  and  displays  narrow  spectral  features. 


Reference  samples  were  heated  using  a  pyroprobe  for  1  s  to  various  temperatures  and 
their  spectra  were  then  compared  with  those  of  the  explosion  sample,  see  Figure  50.  This 
comparison  indicates  a  temperature  of  about  250  °C  (523  K)  for  the  explosion  sample. 
This  reading  compares  with  a  peak  temperature  of  about  390  K  measured  by  the 
thermocouple,  corresponding  to  a  34%  higher  reading.  Since  thermocouples  have  a 
relatively  large  mass  compared  to  the  small  particulate  sensors,  and  an  explosion 
represents  a  very  fast  heating  event,  it  is  fully  expected  that  the  Eu:Y203  temperature 
sensors  yield  higher  reading  than  the  relatively  slow  thermocouple. 


Excitation  Wavelength  (nm) 

Figure  50.  Comparing  excitation  spectra  of  Eu:Y203  subjected  to  explosive  heating  with 
those  of  Eu:Y203  heated  via  a  pyroprobe  (left)  and  temperature  dependence  of  peak 
wavelength  and  peak  width. 
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Another  explosion  test  focused  on  identifying  any  dependencies  for  samples  collected  at 
different  distances  from  the  charge  location.  The  sample  material  was  hung  about  8.5” 
away  from  an  explosive  charge.  A  thermocouple  was  also  placed  about  8.5”  away  from  the 
charge.  Several  collection  bins  were  placed  at  different  distances  from  the  sample 
material.  Bins  labeled  “a”,  “b”,  and  “c”  collected  materials  at  0-18”,  12-30”,  and  24-40”, 
respectively.  Figure  51  shows  the  experimental  setup  and  the  temperatures  measured  by 
the  thermocouple.  Six  different  tests  were  conducted. 


Sample  # 

Peak  Temp.  (K) 

1 

436 

2 

442 

3 

397 

4 

506 

5 

556 

6 

615 

Figure  51.  Explosion  test  setup  (left)  and  temperatures  measured  via  thermocouple 
(right). 

The  temporal  evolution  of  the  temperature  measured  by  the  thermocouple  is  shown  in 
Figure  52.  The  maximum  is  reached  after  about  0.1  s.  However,  the  thermocouple  is 
known  to  not  be  able  to  accurately  track  the  temporal  evolution  of  an  explosive  fireball, 
which  is  significantly  faster. 


Figure  52.  Thermocouple  response  during  the  explosions. 
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Figure  53  shows  photos  of  the  six  post-detonation  sensor  materials.  Vials  “a”,  “b”,  and  “c” 
correspond  to  the  different  collection  distances.  The  sensor  material  samples  contain 
various  impurities  from  the  explosion.  However,  in  all  six  cases,  samples  “b”  appear  to 
have  the  least  amount  of  impurities. 


Figure  53.  Photos  of  post-detonation  sensor  materials  for  six  explosion  test. 


Table  7Table  7  lists  the  temperatures  determined  by  the  nanophase  sensors.  Because  of 
the  inhomogeneity  within  each  sample,  three  small  sample  amounts  were  extracted  from 
each  vial  and  measured  separately.  These  samples  are  labeled  “A”,  “B”,  and  “C”.  In 
addition,  Table  7Table  7  lists  the  average  temperature  as  well  as  the  temperature 
measured  by  the  thermocouple. 
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Table  7.  Comparison  of  temperatures  measured  via  particulate  sensors  based  upon 


FWHM  al 

611  nm  and  thermocouples  (TC). 

Temperature  (K) 

Average  Temp. 

(K) 

TC  Temperature  (K) 

Sample 

A 

B 

C 

la 

803 

721 

717 

747±49 

436 

ib 

633 

652 

611 

632±21 

1C 

751 

622 

684 

686±65 

2a 

622 

541 

687 

6i7±73 

442 

2b 

651 

585 

613 

6i6±33 

2C 

690 

624 

673 

662±34 

3a 

551 

597 

541 

563±30 

397 

3b 

536 

660 

582 

593±63 

3C 

608 

542 

553 

568±35 

4a 

686 

653 

630 

656±28 

506 

4b 

627 

544 

546 

572±47 

4c 

803 

745 

627 

725±90 

5» 

780 

* 

864 

822±59 

556 

5b 

653 

762 

691 

702±55 

5c 

752 

789 

777 

773±19 

6a 

824 

863 

844±28 

615 

6b 

598 

586 

623 

6o2±i9 

6c 

657 

690 

802 

7i6±76 

*  Materia' 

[  is  fully  crystallized,  i.e.  temperature  was  too  high  1 

:or  this  sensor  material. 

Except  for  sample  3,  the  “b”  samples  always  read  the  lowest  temperature.  This  may  be 
explained  with  the  expected  temperature  profile  within  an  explosive  fireball.  Also,  in  all 
cases  the  maximum  temperature  recorded  by  the  thermocouple  is  lower  than  the 
temperature  recorded  by  the  nanosensors.  The  slow  response  time  of  thermocouples 
prevents  them  from  following  a  fast  temperature  profile  such  as  is  expected  within  an 
explosive  fireball.  As  is  apparent  from  the  data  in  Table  7,  the  nanosensor  materials  have 
a  much  faster  response  time.  The  variation  in  measured  temperature  reflects  the 
inhomogeneous  temperature  distribution  within  an  explosive  fireball.  The  difference 
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between  the  two  techniques  ranges  anywhere  from  17%  to  58%,  with  the  higher 
temperature  measurements  showing  less  of  a  difference. 


A  third  explosion  test  was  conducted  inside  a  closed  chamber,  see  Figure  54.  Also  shown 
are  photos  of  the  debris  collected  from  these  tests. 


Figure  54.  Closed  chamber  configuration  with  locations  of  three  thermocouples  (left) 
and  photos  of  debris  collected  after  the  explosion  (right). 


In  our  previous  measurements  for  the  open-air  tests,  we  used  the  FWHM  of  the 
amorphous  peak  to  determine  the  temperature,  see  Figure  55.  The  temperature  inside  the 
closed  chamber  was  so  high  that  all  the  temperature  sensor  material  crystallized.  As  such, 
we  used  the  FWHM  of  the  crystalline  peak  to  determine  the  temperature,  see  Figure  55. 


TC  -  Temperature  /  K  Excitation  wavelength  /  nm 


Figure  55.  Temperature-dependent  FWHM  of  amorphous  peak  (left)  and 
temperature-dependent  crystalline  peaks  (right). 


Using  our  temperature  sensors,  we  determined  a  temperature  of  about  1020  K.  The  three 
thermocouples  measured  395  K,  750  K,  and  1060  K,  respectively. 
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Eu:Zr02 

Characterization  of  as-sunthesized  material 

Nanophase  Zr02:Eu  precursor  is  prepared  by  precipitation  as  described  by  Su  et  al.20  Zr02 
has  three  polymorphs  at  ambient  pressure,  monoclinic  phase  (m-phase)  at  low 
temperatures,  tetragonal  phase  (t-phase)  above  about  1200  °C,  and  cubic  phase  (c-phase) 
at  temperatures  above  about  2370  °C.21  Because  the  tetragonal  and  cubic  phases  are  of 
high  technological  interest,  many  approaches  have  been  reported  in  attempting  to 
stabilize  either  of  these  phases  at  ambient  conditions.22  Most  of  these  approaches  include 
doping  with  trivalent  cations  leading  to  the  formation  of  oxygen  vacancies.  These 
vacancies  are  believed  to  be  responsible  for  the  stabilization  of  the  tetragonal  or  cubic 
phase.23  Due  to  surface  energy  effects,  oxygen  vacancies  can  be  stabilized  in  undoped 
isolated  nanoparticles  of  less  than  10  nm  diameter  and  in  aggregated  nanoparticles  of  less 
than  33  nm  diameter.  24  The  tetragonal  phase  can  also  be  stabilized  by  the  presence  of 
anions  such  as  OH.24  The  metastable  Zr02  phases  may  also  be  stabilized  by  addition  of 
solutes,  e.g.,  Y203,  MgO,  and  Ce203.25-27 

The  optical  properties  of  rare-earth-doped  Zr02  have  been  characterized  by  several 
authors.  Of  particular  interest  is  the  ability  to  use  rare-earth  dopants  such  as  Eu3+  to 
monitor  structural  changes.28^0  For  example,  the  ratio  of  the  electric-dipole  allowed 
sDo— >7F2  transition  to  the  magnetic-dipole  allowed  5D0— >7Fi  transition  can  be  used  to 
extract  changes  about  the  symmetry  of  the  Eu3+  site.  The  lower  the  site  symmetry,  the 
higher  the  intensity  of  the  sDo— >7F2  transition  and  the  higher  the  ratio.  Also,  the  splitting 
of  the  Eu3+  energy  levels  will  be  different  in  monoclinic,  tetragonal,  and  cubic  symmetry. 
Higher  symmetry  leads  to  degeneration,  while  lower  symmetry  removes  the  degeneracy 
and  leads  to  the  observation  of  additional  spectral  lines. 

Processing  completeness  requires  that  Zr02  nanoparticle  precursors,  prepared  via  wet- 
chemistry  techniques,  undergo  a  subsequent  calcination  step  to  remove  hydrates, 
carbonates,  ammonium,  etc.,  and  to  induce  crystallization.  The  decomposition  and 
crystallization  processes  lead  to  variations  in  the  chemical  surrounding  of  the  embedded 
europium  ions  in  the  Zr02  precursor  matrix.  These  structural  changes  affect  the  optical 
properties  of  Eu3+  as  the  matrix  undergoes  various  irreversible  phase  transformation 
from  amorphous  to  metastable  t-phase  and  to  the  stable  m-phase  with  increasing 
temperatures.31  According  to  Ray  et  al,  the  drop-wise  addition  of  Zr0Cl2-8H20  aqueous 
solution  to  NH4OH  aqueous  solution  yields  Zr(OH)4;32 

ZrOCl2  +  2NH4OH  +  H20  Zr(OH)4  +  2NH4CI 

The  precipitated  powder  is  composed  of  monodispersed  nanosized  particles.  Due  to  their 
tendency  to  aggregate,  precise  measurement  of  the  particle  size  is  ambiguous. 

The  TEM  image  shown  in  Figure  56  (a)  suggests  a  particle  diameter  of  about  2-5  nm.  The 
XRD  spectrum,  shown  in  Figure  56  (b),  shows  broad  halo  peaks,  indicating  that  the  as- 
prepared  precipitated  powder  is  amorphous.  The  FTIR  spectrum  of  the  Zr02  precursor, 
shown  in  Figure  56  (c),  is  consistent  with  the  presence  of  lattice  water.  The  broad  peak  at 
3440  cm-1  and  the  narrow  peak  at  1635  cm-1  are  caused  by  the  symmetric  and  asymmetric 
H-O-H  stretching  and  deformation  vibrations,  respectively. 33  The  shoulder  at  about  1580 
cm1  and  the  peak  at  1350  cm1  are  attributed  to  carbonates  that  are  produced  during 
exposure  to  C02  in  the  atmosphere.34 
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2  Theta  Wavenumber  (cm1) 

Figure  56.  (a)  TEM  image,  (b)  XRD  spectrum,  and  (c)  FTIR  spectrum  of  as-prepared 
Zr(OH)4. 


Furnace  heating. 

Figure  57  (a)  shows  the  XRD  spectra  of  the  Zr(OH)4  precursor  when  heat-treated  for  3 
hours  at  temperatures  between  300  °C  and  900  °C.  Heating  to  300  °C  for  3  hours  does 
not  lead  to  crystallization  of  the  precursor.  When  heated  to  400  °C  for  3  hours,  some  XRD 
peaks  become  visible,  which  can  be  matched  to  a  mixed  ZrC>2  t-phase  and  m-phase. 


2  Theta 


Figure  57.  (a,  b)  XRD  and  FTIR  spectra  of  Zr(OH)4  heat-treated  for  3  hours  at 
temperatures  between  300  °C  and  900  °C.  The  gray  and  black  histograms  mark  the 
positions  of  monoclinic  (PDF  No.  072-0597)  and  tetragonal  peaks  (PDF  No.  088-1007), 
respectively. 


Figure  57  (b)  shows  Zr-0  vibrational  bands  at  495  cm-1,  575  cm-1,  and  745  cm-1  (indicated 
by  dashed  lines)  for  samples  heated  to  400  °C  and  above,  confirming  the  presence  of  the 
crystalline  ZrCX  m-phase  and  t-phase. 35, 36  Figure  57  (b)  also  shows  that  even  with  heat 
treatment  at  900  °C,  the  presence  of  OH  and  CO3  groups,  due  to  reactions  with  the 
atmosphere,  can  still  be  seen.  Following  calcination,  Zr(OH)4  has  been  proposed  to 
undergo  continuous  dehydroxylation;37 
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Zr(OH)4  -»•  ZrOi.5(OH)  +  1.5H2O 
ZrOl.sCOH)  — >  ZrOz  +  0.5H2O 

The  crystallite  size  and  the  percent  composition  of  the  t-phase  with  respect  to  the  m- 
phase,  can  be  estimated  from  the  XRD  spectra,  see  Table  8.  As  the  heating  temperature 
increases  from  400  °C  to  900  °C,  the  crystallite  size  increases  from  8  nm  to  33  nm  for  the 
t-phase,  and  from  11  nm  to  25  nm  for  the  m-phase. 


Table  8.  Characteristics  of  the  furnace-heated  ZrCX  crystalline  phases. 


Temperature  (°C) 

t-phase  /  m-phase 
Crystallite  size  (nm) 

t-phase/ 
m-phase  (%) 

400 

8/11 

45/55 

500 

13/16 

28/72 

700 

18/25 

14/86 

900 

33/25 

11/89 

At  400  °C,  the  crystallite  sizes  are  in  the  range  for  the  t-phase  to  m-phase  transformation 
for  isolated  crystallites,  10  nm.24  The  crystallite  sizes  are  slightly  smaller  for  the  t-phase 
which  may  suggest  phase  transformation  of  the  t-phase  to  the  m-phase  as  it  grows  in  size. 
Indeed,  the  amount  of  metastable  t-phase  decreases  with  increasing  temperatures.  As  the 
heat  treatment  temperature  increases  from  400  °C  to  900  °C,  the  percent  of  t-phase 
composition  decreases  from  45%  to  11%  (Table  1).  The  non-vanishing  t-phase  indicates 
that,  in  addition  to  its  role  as  an  optical  probe,  the  europium  dopant  might  act  as  stabilizer 
of  this  phase.38 

Figure  58  shows  the  fluorescence 
spectra,  excited  at  536  nm,  of  Zr02:Eu 
precursors  furnace-heated  for  3  hours  to 
various  temperatures.  The  peaks  at  592 
nm  and  607  nm  can  be  assigned  to  the 
sDo— >7Fi  and  sDo— >7F2  transitions  of  the 
t-phase.39  The  remaining  peaks  belong  to 
the  m-phase,  which  has  a  lower 
symmetry  than  the  t-phase  and  can  lead 
to  splitting  of  the  transitions.  At  300  °C 
only  emission  from  amorphous 
components  are  observed.  As  the 
temperature  increases,  a  mix  of  m-phase 
and  t-phase  forms,  just  as  observed  in  the 
XRD  spectra.  At  400  °C  the  percentage  of 
t-phase  is  significantly  higher  than  at 
900  °C,  confirming  that  the  t-phase  is 
metastable  and  converts  into  the  m- 
phase  as  the  temperature  increases. 


Wavelength  (nm) 

Figure  58.  Fluorescence  spectra  of 
Eu:Zr02  precursor,  furnace-heated  for  3 
hours  to  various  temperatures. 
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.Rapid  puroprobe  heating  of_Zr[_ OH)d 

Figure  59  (a)  and  (b)  show  the  XRD  and  FT-IR  spectra  of  Zr(OH)4  when  heat  treated  for 
10  s  to  temperatures  between  420  °C  and  835  °C.  The  samples  heated  to  420  °C  and  495 
°C  do  not  show  any  crystalline  Zr02  signatures,  while  samples  heated  to  higher 
temperatures  show  clear  crystalline  signatures.  Just  as  in  the  case  of  the  furnace -heated 
samples,  the  presence  of  the  mixed  t-phase  and  m-phase  leads  to  Zr-0  vibrational  peaks 
495  cnr1,  575  cm-1,  and  745  cm’1,  (indicated  by  dashed  lines  in  Figure  5b). 


Figure  59.  (a,  b)  XRD  and  FTIR  spectra  of  Zr(OH)4  heat-treated  at  various 
temperatures,  300-900  °C  for  10  s.  The  gray  and  black  histograms  mark  the  positions 
of  monoclinic  (PDF  No.  072-0597)  and  tetragonal  peaks  (PDF  No.  088-1007), 
respectively. 


The  XRD  spectra  of  Zr(OH)4  heat  treated  for  10  s  in  the  pyroprobe  differ  from  those  heat 
treated  for  3  hours  in  the  furnace  (compare  Figure  57  a  and  Figure  59  a).  Both  heating 
procedures  result  in  a  mixture  of  m-phase  and  t-phase.  However,  for  the  pyroprobe- 
heated  samples,  the  peaks  indicate  a  stronger  presence  of  the  metastable  t-phase  and  for 
the  furnace-heated  samples,  the  peaks  indicate  a  stronger  presence  of  the  m-phase.  In 
comparison  with  the  furnace-heated  samples  that  are  ramped  up  at  a  rate  of  25  °C/min, 
the  rapidly  heated  samples  (nominal  ramping  rate  of  the  pyroprobe  heated  sample  is 
20000  °C/s)  may  not  have  sufficient  time  to  completely  dehydroxylate  the  Zr(OH)4 
compound  and  to  achieve  large  crystallites.  For  the  matching  crystallite  sizes  of  the 
pyroprobe  and  the  furnace  heated  samples,  the  pyroprobe  heated  sample  clearly  shows  a 
larger  amount  of  the  metastable  t-phase.  The  remaining  OH  groups  may  have  stabilized 
the  t-phase  in  the  pyroprobe-heated  samples  (compare  Table  8  and  Table  9).  In  contrast 
to  the  3-hours  furnace-heated  samples,  the  crystallite  sizes  of  the  t-phase  are  larger  than 
the  ones  of  the  m-phase.  This  observation  suggests  that  for  rapid  heating  the  metastable 
t-phase  is  favored  over  the  m-phase. 
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Table  9.  Characteristics  of  the  pyroprobe-heated  ZrQ2  crystalline  phases. 


Temperature  (°C) 

t-phase/m-phase 
Crystallite  size  (nm) 

t-phase/ 
m-phase  (%) 

570 

10/6 

68/32 

710 

11/6 

68/32 

835 

14/11 

64/36 

TEM  images  confirm  the  significant  role  of  temperature  on  the  ZrC>2  crystallite  growth, 
even  at  heating  times  as  short  as  10  s  (Figure  60).  The  observed  particle  sizes  are  on  par 
with  those  estimated  from  the  XRD  spectra.  It  is  worth  noting  that  the  particle  sizes  do 
not  necessarily  correspond  to  the  crystallite  size  due  to  its  polycrystalline  nature. 


Figure  60.  TEM  micrographs  of  Zr(OH)4  calcined  at  various  temperatures  for  10 
seconds:  (a)  420  °C  (b)  495  °C  (c)  570  °C,  (d)  710  °C,  and  (e)  835  °C.  Scale  bar  is  20  nm. 

Figure  61  shows  the  fluorescence  spectra,  excited  at  536  nm,  of  Zr02:Eu  precursors 
pyroprobe-heated  for  10  s  to  various  temperatures.  At  the  two  lowest  temperatures,  only 
amorphous  components  are  present  and  emit  light.  As  the  temperature  reaches  570  °C, 
crystalline  emission  peaks  appear,  which  show  the  presence  of  the  t-phase  and  m-phase. 
Compared  to  the  furnace-heated  samples,  the  t-phase  has  a  much  stronger  presence  in 
these  samples.  Again,  the  peaks  at  592  nm  and  607  nm  can  be  assigned  to  the  sDo— >7Fi 
and  5Do— >7F2  transitions  of  the  t-phase,  and  the  remaining  peaks  to  the  m-phase.  As  the 
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temperature  reaches  835  °C,  the  signal  for  the  m-phase  increases  compared  to  those  at 
lower  temperatures. 


Besides  Eu-doped  Y2O3,  we  recently  also 
presented  some  initial  results  for  Eu- 
doped  ZrC>2.3  The  synthesis, 
morphological  characterization,  and 
pyroprobe  heating  of  the  Eu-doped 
Zr(OH)4  nanoparticles  is  described 
elsewhere.3  Briefly,  an  aqueous  solution 
of  Zr0Cl2-8-H20  and  Eu(N03)3-5-H20  is 
added  to  an  aqueous  solution  of 
ammonium  hydroxide  to  reach  a  pH  of 
10.  The  white  gel  that  forms  is  aged  for  2 
h  and  is  then  several  times  centrifuged 
and  redispersed  in  ionized  water. 
Subsequently,  the  material  is  rinsed  in 
acetone,  dried  for  12  h  at  80  °C,  and  then 
crushed  into  a  powder.  Using  this 
approach,  Euo.01Zro.99O2  was  synthesized 
and  used  for  the  following  characterization. 


Figure  61.  Fluorescence  spectra  of  Eu:Zr02 
precursor,  pyroprobe-heated  for  10  s  to 
various  temperatures. 


We  previously  determined  that  the  Euo.oiZro.99(OH)4  nanoparticles  crystallize  into 
tetragonal  and  monoclinic  phases  upon  heating.3  Figure  62  (left)  shows  the  fluorescence 
spectra  of  Euo.oiZro.99(OH)4  samples  that  have  been  heated  for  three  hours  in  a  furnace  at 
673  K  and  at  1123  K,  respectively. 


Wavelength  (nm)  Wavelength  (nm) 

Figure  62.  Fluorescence  spectra  of  Euo.oiZro.99(OH)4  heated  for  three  hours  to  673 
K  and  1123  K,  respectively  (left),  and  heated  for  10  s  to  693  K  and  1108  K,  respectively 
(right). 


While  the  spectra  are  normalized  with  respect  to  each  other,  it  should  be  kept  in  mind 
that  the  emission  intensity  of  the  sample  heated  to  1123  K  is  significantly  stronger  than 
the  emission  intensity  of  the  sample  heated  to  673  K.  These  two  spectra  show  some  of  the 
main  features  of  Eu  ions  embedded  in  Zr02.  The  peaks  at  about  592  nm  and  at  about  606 
nm  correspond  to  the  sD0  — ►  7Fi  and  sD0  — ►  7F2  transitions  of  Eu  ions  (D2d  site  symmetry) 
in  tetragonal  Zr02  (P42/n  m  c),  respectively.40  The  remaining  peaks  at  about  597  nm  and 
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610  nm  through  635  nm  correspond  to  the  sDo  — ►  ?Fi  and  sDo  —>  7F2  transitions  of  Eu  ions 
(Ci  site  symmetry)  in  monoclinic  Zr02  (Pi  21/c  i).4° 

Shown  in  Figure  62  (right)  are  the  spectra  of  Euo.oiZr0.99(OH)4  samples  heated  for  10  s  in 
the  pyroprobe  heater  to  693  K  and  1108  K,  respectively.  The  spectral  features  of  the 
sample  heated  to  693  K  are  relatively  broad,  indicating  that  most  of  the  sample  is  still 
amorphous.  The  crystalline  spectral  features  become  apparent  somewhere  between  770 
K  and  840  K.3  Compared  to  the  furnace-heated  samples,  the  spectrum  of  the  pyroprobe- 
heated  sample  (1108  K)  shows  a  stronger  signal  for  the  emission  peaks  of  the  tetragonal 
phase  compared  to  the  monoclinic  phase.3 

Figure  63  shows  the  ratio  of  the  intensities  of  the  peaks  at  606  nm  and  614  nm  for  a  range 
of  temperatures  for  the  3  h  furnace-heated  sample  and  the  10  s  pyroprobe-heated  sample. 
For  the  furnace-heated  samples  the  ratios  show  a  linear  dependency  with  temperature, 
providing  a  means  to  determine  temperature  based  upon  the  ratio  of  tetragonal  to 
monoclinic  composition.  However,  no  ratio/temperature  relationship  is  apparent  for 
short  10  s  pyroprobe  heating,  see  Figure  63  (right).  Thus,  alternative  temperature 
indicators  need  to  be  determined. 


Temperature  (K) 


Figure  63.  Tetragonal/monoclinic  intensity  ratios  for  3  h  furnace  heating  (left)  and 
10  s  pyroprobe  heating  (right). 

Our  investigation  of  heated  Euo.o2Yi.98(C03)3  indicates  that  the  peak  width  and  peak 
position  of  excitation  spectra  are  good  temperature  indicators. 41.  42  Thus,  excitation 
spectra  of  heated  Euo.oiZro.99(OH)4  are  measured  between  510  nm  and  590  nm  for 
samples  heated  for  10  s,  and  between  510  nm  and  600  nm  for  samples  heated  1  s  to  various 
temperatures,  while  monitored  at  606  nm  (tetragonal  phase)  and  614  nm  (monoclinic 
phase),  respectively,  see  Figure  64  and  Figure  65.  The  samples  monitored  at  606  nm 
shows  peaks  between  520  nm  and  550  nm  corresponding  to  the  "F0  — >  sDi  transitions, 
and  peaks  between  575  nm  and  590  nm  corresponding  to  the  ?F0  — ►  sDo  transitions, 
respectively.  Figure  64  (right)  which  shows  the  excitation  spectra  monitored  at  614  nm, 
shows  some  of  the  same  peaks  as  shown  in  Figure  64  (left)  which  are  monitored  at  606 
nm,  indicating  that  the  emission  at  614  nm  may  not  be  purely  due  to  the  monoclinic  phase 
but  may  also  contain  contribution  from  the  tetragonal  phase. 
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Figure  64.  Excitation  spectra  of  Euo.oiZro.99(OH)4  heated  for  10  s  to  various 
temperatures  and  monitored  at  606  nm  -  tetragonal  phase  -(left)  and  614  nm  - 
monoclinic  phase  (right). 


Figure  65  shows  the  excitation  spectra  of  samples  heated  for  1  s.  These  spectra  are  very 
similar  to  the  spectra  of  samples  heated  for  10  s  and  shown  in  Figure  64,  indicating  that 
even  for  heating  times  as  short  as  1  s  some  of  the  material  crystallizes.  However,  the  peaks 
of  the  samples  monitored  at  614  nm,  see  Figure  65  (right)  remain  relatively  broad 
compared  to  the  samples  heated  for  10  s  and  the  samples  heated  to  1  s  but  monitored  at 
606  nm.  This  observation  could  be  explained  if  the  tetragonal  phase  starts  to  crystallize 
before  the  monoclinic  phase.  The  higher  the  temperature  and  the  longer  the  heating  time, 
the  more  intense  the  peaks  between  520  nm  and  540  nm  compared  to  the  peaks  between 
575  nm  and  590  nm.  This  observation  is  more  pronounced  for  the  excitation  spectra 
monitored  at  606  nm  (tetragonal)  than  the  ones  monitored  at  614  nm  (monoclinic),  and 
is  consistent  with  the  fact  that  the  5D0  — ►  ?Fo  transition  is  forbidden  in  D2d  symmetry 
(tetragonal)  but  electric  dipole  allowed  in  Ci  symmetry  (monoclinic).  However,  unlike  for 
heated  Euo.o2Yi.98(C0.3)3,  the  excitation  spectra  of  heated  Euo.oiZro.99(OH)4  do  not  readily 
show  any  apparent  changes  in  peak  position  or  peak  width  that  could  be  used  as  a 
temperature  indicator.  Therefore,  the  fluorescence  lifetime  and  fluorescence  spectra  are 
evaluated  as  potential  temperature  indicators. 
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Figure  65.  Excitation  spectra  of  Euo.oiZro.99(OH)4  heated  for  1  s  to  various 
temperatures  and  monitored  at  606  nm  -  tetragonal  phase  -  (left)  and  monitored  at  614 
nm  -  monoclinic  phase  -(right). 


58 


Figure  66  shows  the  fluorescence  lifetimes  of  Eu-doped  Zr02  heated  for  10  s  to  various 
temperatures,  excited  at  536  nm,  and  monitored  at  606  nm  and  614  nm,  respectively.  The 
decay  curves  of  the  samples  monitored  at  614  nm  (monoclinic)  show  a  very  short  lifetime 
component  as  well  as  a  longer  lifetime  component.  These  curves  are  fitted  with  a  single 
exponential  for  all  temperatures  starting  at  0.3  ms.  The  lifetime  data  for  the  samples 
monitored  at  606  nm  (tetragonal)  are  fitted  with  a  single  exponential  for  temperatures 
below  700  K  (amorphous  phase).  Above  700  K,  the  data  are  fitted  with  a  double 
exponential.  One  of  the  exponentials  is  the  same  as  the  monoclinic  lifetime,  indicating 
that  a  mixed  phase  is  observed  at  606  nm.  The  second  exponential  (tetragonal)  is  the  one 
shown  in  Figure  66  (left).  As  expected,  at  temperatures  below  about  700  K,  when  the 
material  is  still  amorphous,  the  observed  lifetimes  are  about  the  same  for  both 
observation  wavelengths.  As  the  temperature  is  further  increased,  the  tetragonal  phase 
displays  a  jump  in  fluorescence  lifetime  from  about  1  ms  to  more  than  4  ms.  Subsequently, 
the  lifetime  increases  a  little  more  before  it  drops  again  to  about  3  ms  at  about  1400  K. 
The  monoclinic  phase  on  the  other  hand  displays  a  different  temperature  dependence. 
The  lifetime  continues  to  increase  more  slowly  to  a  maximum  of  about  1.6  ms,  at  which  it 
levels  out.  Only  at  temperatures  above  1200  K  does  the  lifetime  decrease  again  to  about 
1.5  ms  at  about  1300  K.  The  increase  and  subsequent  decrease  of  fluorescent  lifetime  with 
temperature  is  not  unusual  in  nanoparticles.  A  similar  behavior  is  also  seen  in  Eu:Y203,2 
and  the  concept  of  this  effect  was  described  by  Meltzer  et  al.43  Equation  35  describes  the 
dependence  of  the  radiative  lifetime  of  an  electronic  transition: 


1.5X104  Ig 
f(ED)  [i(n2+2)]2n 


(39) 


where  f(ED)  is  the  oscillator  strength,  n  is  the  index  of  refraction,  and  X  is  the  wavelength 
in  vacuum.  Because  the  nanoparticles  are  much  smaller  than  the  wavelength  of  light,  their 
immediate  surrounding  environment  needs  to  be  considered  to  determine  an  effective 
index  of  refraction.  In  our  case,  the  nanoparticles  start  out  as  Zr(OH)4  and  decompose 
into  Zr02  as  they  are  heated.  During  the  decomposition  process  the  size,  shape,  and 
surrounding  of  the  particles  change,  resulting  in  changes  in  the  radiative  lifetime.  A  more 
detailed  description  of  this  effect  can  be  found  in  the  literature.2’  43  The  fact  that  the  shape 
of  the  fluorescence  lifetime  data  for  Eu  ions  in  the  tetragonal  phase  is  different  from  the 
shape  of  the  fluorescence  lifetime  data  of  Eu  ions  in  the  monoclinic  phase  indicates  that 
their  morphological  development  must  be  different  during  the  heating  process. 
Considering  the  shape  of  the  fluorescence  lifetime  curves,  the  data  monitored  at  614  nm 
could  provide  some  indication  for  temperature  in  the  range  between  about  500  K  and 
about  1100  K. 
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Figure  66.  Fluorescence  lifetime  of  Euo.oiZro.99(OH)4  heated  for  10  s,  excited  at  536 
nm,  and  monitored  at  606  nm  (left)  and  614  nm  (right),  respectively. 


Figure  67  (left)  shows  a  small  section  with  two  peaks  of  the  fluorescence  spectra  of 
Euo.oiZr0.99(OH)4  heated  for  10  s  and  1  s  to  various  temperatures.  Both  of  these  peaks  have 
the  same  fluorescence  lifetime,  indicating  that  they  are  both  resulting  from  the  same  level 
of  Eu  ions  in  the  tetragonal  phase.  It  is  apparent  that  the  shorter  wavelength  peak 
decreases  in  intensity  as  the  heating  temperature  is  increased.  Figure  67  (right)  shows  the 
ratio  of  the  emission  peak  intensities  as  a  function  of  temperature.  As  can  be  seen,  a  linear 
relationship  exists  between  peak  ratio  and  temperature  in  the  range  between  about  800 
K  and  about  1300  K.  At  lower  temperatures  the  material  is  still  amorphous  and  the  peaks 
are  not  well  characterized.  At  temperatures  above  about  1300  K,  the  curve  flattens  out. 
The  data  from  the  1  s  heating  tests  overlaps  well  with  the  data  from  the  10  s  heating  tests. 


Wavelength  (nm) 

Figure  67.  Fluorescence  spectra  of  Euo.oi! 
temperatures  (left)  and  peak  ratio  (right). 


Temperature  (K) 

0.99(011)4  heated  for  10  s  and  1  s  to  various 


The  changing  intensity  peak  ratio  with  temperature  is  most  likely  due  to  the  continued 
crystallization  process  of  the  tetragonal  phase  over  the  observed  temperature  range.  The 
transition  from  an  amorphous  phase  to  a  crystalline  phase  with  higher  symmetry  leads  to 
the  reduction  of  energy  level  splittings.  This  process  starts  at  the  crystallization 
temperature  and  ends  when  all  material  has  been  crystallized. 
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D.  Temperature/time  (Thermal  impulse)  Sensors 

In  order  to  simultaneously  extract  information  about  heating  temperature  and  heating 
duration,  we  rely  on  a  kinetic  analysis  of  our  sensors.  The  theoretical  background  and 
some  modeling  information  was  provided  in  section  A.  This  approach  requires  that  we 
use  two  different  sensor  materials  with  sufficiently  different  kinetic  parameters.  Ideally, 
these  two  sensor  materials  are  physically  combined  in  such  a  manner  that  they  always 
measure  and  record  in  the  same  location.  We  prepared  different  types  of  core/shell 
particles  and  evaluated  different  sensing  approaches. 

Initial  evaluation  of  various  compounds 


For  our  first  type  of  core/shell  sensor,  we  combine  the  well-known  Zr02  and  Y2O3 
matrices  that  work  well  for  crystallization-based  changes.  However,  in  order  to  probe 
temperature-induced  changes  in  these  two  materials,  we  have  to  use  two  different  rare- 
earth  dopants,  and  selected  Er,Yb:Zr02/Eu:Y203  for  evaluation.  The  Er,Yb  combination 
yields  upconverted  light  in  the  ref  and  green  spectral  regions  when  excited  with  infrared 
980  nm  light.  Depending  on  the  Er/Yb  molar  ratio,  the  ratio  of  the  red  and  green  emission 
intensity  shows  a  linear  temperature  dependence  for  furnace  heating,  see  Figure  68. 
However,  when  heated  for  short  periods  of  time  via  a  pyroprobe,  the  temperature 
dependence  was  no  longer  apparent.  Furthermore,  it  appears  that  the  shell  prevents  the 
growths  of  Zr02,  limiting  its  functionality. 


Temperature  /  °C 


Figure  68.  Temperature  dependence  up  upconverted  red  and  green  emission 
intensities  for  furnace  heating  (left)  and  pyroprobe  heating  (right). 


Next,  we  evaluated  potential  temperature-induced  diffusion  and  resulting  upconversion 
in  Yb203/Er:Y203.  Ytterbium  is  a  well-known  activator  for  upconversion  and  if  enough  Yb 
diffuses  from  the  core  into  the  shell,  the  shell  should  show  enhance  upconversion.  Again, 
we  saw  this  effect  occurring  for  furnace-heated  samples,  but  not  for  pyroprobe-heated 
samples. 


Finally,  we  attempted  to  use  an  intensity  reference.  In  particular,  we  used 
Tb:Zr02/Eu:Y203,  in  which  the  core  was  fully  crystalline  and  the  shell  was  a  temperature 
sensitive  precursor.  During  actual  heating  events,  the  core  would  no  longer  change  and 
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could  be  used  as  a  reference.  Figure  69  shows  the  spectra  and  intensity  ratios  for 
pyroprobe-heated  samples.  These  results  demonstrate  the  promising  nature  of  using  an 
intensity  reference. 


Figure  69.  Temperature  dependent  fluorescence  spectra  (left),  corresponding  peak 
intensity  ratios  (center),  and  corresponding  integrated  intensities  (right). 


Next,  we  evaluate  c-Tb:Zr02/p-Eu:Y203.  In  order  to  measure  two  unknowns,  temperature 
(T)  and  time  (t),  the  sensor  needs  to  contain  two  materials  with  sufficiently  different 
decomposition  characteristics.  Also,  due  to  potential  temperature  gradients  within  the 
fireball,  the  two  sensor  materials  need  to  be  able  to  measure  and  record  in  the  same 
location.  Furthermore,  since  the  decomposition  characteristics  are  evaluated  by  the 
optical  emission  from  rare-earth  dopants  in  the  sensor  materials,  appropriate  optical 
indicators  with  non-interfering  optical  spectra  need  to  be  identified. 

We  previously  concluded  that  using  sensor 
materials  in  a  core/shell  configuration  is  the  best 
approach  to  ensure  that  they  measure  and  record 
temperature  in  the  same  location  within  the 
fireball.  Multiple  types  of  core/shell  systems  based 
upon  a  variety  of  approaches  including 
crystallization,  diffusion,  and  intensity  were 
evaluated.  Based  upon  those  measurements,  a 
core/ shell/ shell  system  with  a  built-in  fluorescence 
intensity  reference  was  selected  as  the  most 
promising  approach. 

Figure  70  shows  a  simplified  design  version  that  is 
used  to  evaluate  this  approach.  These  core/shell 
nanoparticles  are  composed  of  a  fully  crystalline  Tbs+-doped  Zr02  core  (calcined  at  1000 
°C  for  30  minutes),  and  a  shell  consisting  of  unheated  amorphous  Eu:Y203  precursors. 
The  shell  was  prepared  by  growing  three  layers  of  Eu3+-doped  Y203  precursor  (p-Eu:Y203) 
around  the  crystalline  Tb:Zr02  core,  resulting  in  Tb:Zr02/p-Eu:Y203  core/shell 
nanoparticles.  These  core/shell  nanoparticles  were  heated  to  various  temperatures  for 
nominal  heating  times  of  10  s  (actual  time  at  maximum  temperature  was  about  1  s)  via 
pyroprobe  heating.  X-ray  diffraction  (XRD)  spectra  show  that  for  high-enough 
temperatures  the  p-Eu:Y203  shell  undergoes  an  amorphous  to  crystalline  phase  changes. 


Figure  70.  Schematic  of  c- 
Tb:Zr02/p-Eu:Y203  core/shell 
design. 
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Figure  71.  Fluorescence  spectra  of  Figure  72.  Integrated  Eu3+  5D0— >?F2 
Tb:Zr02/p-Eu:Y2C>3  core/shell  fluorescence  intensity  as  a  function  of 

nanoparticles  heated  to  various  temperature, 
temperatures  for  10  s. 

Figure  71  shows  the  previously  measured  fluorescence  spectra  of  these  samples  for 
excitation  at  488  nm.  Since  the  Tb:Zr02  core  is  already  fully  crystallized,  its  fluorescence 
intensity  is  not  expected  to  change  due  to  further  heating  for  short  timeframes.  The 
fluorescence  spectra  were  normalized  to  the  intensity  of  the  Tb3+  5D4— >7F5  transition  at 
543  nm.  The  integrated  fluorescence  intensities  of  the  Eu3+  5D0— >7F2  transition  from  606 
nm  to  635  nm  are  then  measured  and  plotted  as  a  function  of  temperature,  see  Figure  72. 
This  integrated  fluorescence  intensity  monotonically  increases  with  temperature  until 
about  800  °C  and  then  slightly  decreases. 

Continuing  with  this  work,  we  are  evaluating  the  effect  of  time  on  the  optical  spectra  of 
Tb:Zr02/p-Eu:Y203  core/shell  nanoparticles.  About  4-7  mg  of  the  core/shell 
nanoparticles  is  encapsulated  in  a  0.0028"  thick  pyrolitic  graphite  sheet  (Panasonic, 
Corp.)  and  heated  using  a  100  Watt  CO2  laser,  Firestar  series  fioo  (Synrad,  Inc.),  to  600 
°C,  700  °C,  800  °C,  and  900  °C.  The  setpoint  temperatures  are  maintained  for  holding 
times  between  2  s  and  60  s.  A  high  speed  PID  temperature  controller,  HA400  (RKC 
Instruments  Inc.),  is  used  to  achieve  and  maintain  an  isothermal  at  the  setpoint 
temperature.  A  type  K  thermocouple,  CHAL005  (Omega  Engineering,  Inc.)  is  used  to 
monitor  the  heating  temperature.  The  heating  is  performed  in  a  vacuum  chamber  with 
reduced  atmosphere  of  about  (5-10)  Torr  to  minimize  oxidation  of  the  graphite  sheet. 
The  heating  rate  is  200  °C/s.  Figure  73  shows  typical  temperature  profiles  of  the  laser- 
based  heating  method  employed  in  this  work. 
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Figure  73.  Typical  temperature  profiles  of 
the  laser-based  heating  method  employed 
in  this  work. 


Time  /  s 

Figure  74.  Relative  fluorescence 
intensity,  I611/I543,  versus  time  for  c- 
Tb:Zr02/p-Eu:Y203  core/shell 

nanoparticles  heated  to  700  °C,  800  °C, 
and  900  °C  for  2  s  -  60  s. 


In  agreement  with  our  previous  data,  crystallization  of  the  Y2O3  phase  occurs  for  heating 
temperatures  of  700  °C,  800  °C,  and  900  °C,  as  indicated  by  the  appearance  of  the  Eu3+ 
5Do^7F2  emission  peak  at  611  nm.2  Heating  to  600  °C  is  not  sufficient  to  initiate 
crystallization,  even  for  holding  times  of  60  s.  Figure  74  shows  the  plots  of  I611/I543. 

The  relative  intensity  of  the  Eu3+  fluorescence,  and  in  particular  the  intensity  of  the 
hypersensitive  5D0— >7F2  transition,  depends  on  various  parameters.  Since  the  4f  electrons 
are  shielded  by  outer  electrons,  interaction  with  the  lattice  is  minimized.  As  such,  it  is 
mostly  the  strength  of  the  electric  field  and  the  symmetry  of  the  ligands  that  determine 
the  optical  properties  of  the  material.44  Sharma  et  al.  have  reported  the  effect  of  solvent, 
host  precursor,  dopant  concentration,  and  crystallite  size  on  the  fluorescence  properties 
of  Eu3+  doped  yttria.s0  In  the  present  case,  changes  in  the  intensities  can  be  ascribed  to 
the  amorphous  to  crystalline  phase  transition,  and  subsequently  nucleation  and  growth 
of  the  crystalline  phase.  In  addition,  an  initial  analysis  of  the  asymmetry  R  of  the  Eu  sites 
via  comparison  of  the  electric  dipole  allowed  sDo— >7F2  transition  and  the  magnetic  dipole 
allowed  5D0— >7Fi  transition  is  shown  in  Figure  22.  These  differences  indicate  that  effects 
due  to  asymmetry  cannot  be  ruled  out. 
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asymmetry  ratio  for  various  heating  fluorescence  intensity,  I611/I543,  versus 
temperatures.  Higher  ratios  indicate  a  time  for  Tb:Zr02/p-Eu:Y203  core/shell 
larger  degree  of  asymmetry.  45  nanoparticles  that  were  heated  to  700 

°C,  800  °C,  and  900  °C  for  2  s  -  60  s 
fitted  with  the  kinetic  model. 

Figure  74  shows  that  for  the  2  s  heating  the  fluorescence  intensity  is  directly  proportional 
with  temperature.  The  value  is  highest  for  the  900  °C  and  lowest  for  the  700  °C  heating 
temperature.  Likewise,  the  equilibrium  fluorescence  intensity  of  the  samples  heated  to 
800  °C  is  higher  than  those  that  were  heated  to  700  °C.  However,  for  the  samples  that 
were  heated  to  900  °C,  the  equilibrium  fluorescence  intensity  is  on  in  par  with  that  of  the 
samples  heated  to  700  °C. 

It  is  believed  that  the  increasing  trend  in  the  fluorescence  intensity  of  the  5D0— >7F2 
transition  in  the  samples  heated  to  700  °C,  800  °C,  and  900  °C  between  2  s  and  60  s 
correspond  to  nuclei  formation  of  the  Eu:Y203  phase,  and  that  the  fluorescence  intensity 
is  directly  proportional  to  the  number  of  nuclei  formed.  At  900  °C  heating  temperature, 
the  nucleation  stage  has  quickly  reached  completion  within  the  2  s  timeframe,  and  the 
kinetic  mechanism  is  subsequently  dominated  by  crystallite  growth.  At  this  stage,  there  is 
little  untransformed  material  for  nuclei  to  form  in.  Instead,  each  of  the  formed  nuclei  fuse 
into  each  other  to  form  larger  crystals.  It  was  reported  that  the  fluorescence  intensity  of 
the  sDo — >-7F2  transition  of  the  Eu3+  ion  is  inversely  proportional  with  crystallite  size. 3°.  46 
As  such,  the  lower  fluorescence  intensity  in  the  samples  heated  to  900  °C  might  reflect 
significantly  larger  crystallite  grain  sizes  compared  to  those  heated  to  700  °C  and  800  °C. 
Alternatively,  the  reduction  in  the  intensities  may  be  due  to  the  asymmetry  effect.  As 
shown  in  Figure  75,  the  sample  heated  to  900  °C  has  a  lower  equilibrium  R  value  than  the 
samples  heated  to  700  °C  and  800  °C.  A  lower  R  value  indicates  higher  symmetry  and 
lower  fluorescence  intensity  due  to  the  electric  dipole  nature  of  this  transition.45 

For  the  samples  heated  to  700  °C  and  800  °C,  where  the  amorphous  to  crystalline 
nucleation  mechanism  dominates,  the  trend  in  the  change  in  the  relative  intensity  of  the 
Eu3+  fluorescence  as  a  function  of  time  has  the  characteristics  of  the  well-known  Johnson- 
Mehl-Avrami-Kolmogorov  (JMAK)  model,  which  has  the  following  general  form:47-so 

<p(t )  =  1  —  e_fct” 
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where,  k  and  n  are  the  Avrami  constants,  and  0(t)  is  the  time-dependent  fraction  of 
transformation.  The  major  characteristic  of  the  JMAK  model  is  the  s-shaped  profile  that 
indicates  a  low  rate  of  transformation  at  the  beginning  and  at  the  end  of  the 
transformation,  but  rapid  transformation  in  between.  The  equation  was  used  to  describe 
the  formation  of  spherical  nuclei  having  a  “(3”  phase  in  a  matrix  having  an  “a”  phase.  The 
above  equation  has  also  been  successfully  applied  to  describe  the  grain  growth  kinetics  of 
nanocrystalline  materials.2’  12  By  assuming  that  the  fluorescence  intensity,  I,  of  the 
sDo— >7F2  transition  of  the  Eu3+  ion  for  the  samples  that  were  heated  to  700  and  800  °C  is 
proportional  to  the  fraction  of  transformed  phase,  we  write  Equation  36  as  follows: 

±^L=l_e-fctn  (41) 

/00-/0 

where,  I  is  the  instantaneous  fluorescence  intensity,  I0  is  the  initial  intensity  at  the  start 
of  nucleation,  which  is  taken  to  be  zero,  and  I=o  is  the  equilibrium  (limiting)  intensity  at 
that  temperature. 

At  the  onset  of  nucleation  the  relative  fluorescence  intensity  at  611  nm  is  close  to  zero,  see 
Figure  76.  Over  time,  the  fluorescence  intensity  evolves  towards  a  maximum  equilibrium 
(limiting)  value.  In  between  the  nucleation  and  the  equilibrium  stage,  the  fluorescence 
intensity  increases  rapidly.  Figure  76  shows  that  equation  2  fits  well  to  the  plot  of  relative 
intensity  versus  temperature  of  the  samples  that  were  heated  to  700  °C  and  800  °C 
between  2  s  and  60  s.  The  former  yields  k7ooc=o.o8o2,  n=i.03,  and  Lo=o.54  with 
R2=o.998,  and  the  latter  yields  ksooc=o.i9  and  n=i.os  with  R2=o.98.  It  can  be  seen  that 
ksooc  >  k7ooc  which  is  consistent  with  increased  nucleation  rate  at  higher  temperature. 

Work  is  underway  to  verify  the  suitability  of  the  kinetic  model,  and  to  generate  additional 
data  in  the  vicinity  of  the  700  °C  and  800  °C  heating  temperatures  to  yield  a  temperature- 
dependent  Arrhenius  plot,  which  is  needed  to  obtain  the  complete  time-temperature 
dependent  kinetic  equation.  In  addition,  the  effects  of  asymmetry  and  grain  size  need  to 
be  further  evaluated.  Finally,  additional  indicators  such  as  excitation  spectra  need  to  be 
evaluated  to  include  lower  temperature  regimes  in  the  analysis. 

As  part  of  this  evaluation,  we  also  characterized  p-Tb:Zr02/p-Eu:Y2C>3  and  p-EmZrCb/p- 
Tb:Y203.  The  spectra  shown  in  Figure  71  indicate  that  the  emission  intensity  of  the  Eu3+ 
5Do— <-7Fi  transition  (about  580  -  600  nm  range)  appears  to  be  constant.  This  observation 
is  not  completely  surprising,  as  this  particular  transition  is  magnetic  dipole  allowed  and 
as  such  not  sensitive  to  symmetry  changes  in  the  host  material.  It  does  however  indicate 
that  it  might  be  possible  to  use  this  particular  transition  as  a  built-in  intensity  reference 
which  would  eliminate  the  need  for  a  separate  reference  material. 

In  order  to  evaluate  the  potential  of  using  the  Eu  5D0— >7Fi  transition  as  a  built-in 
reference,  we  synthesized  p-Tb:Zr02/p-Eu:Y203  and  p-Eu:Zr02/p-Tb:Y203  core/shell 
nanoparticles  and  evaluated  them  for  a  large  (T,  t)  space  using  the  newly  operational 
laser-heating  setup  with  the  goal  to  perform  a  kinetic  analysis  on  these  sensors. 
Temperature  and  times  of  the  (T,t)  space  include  (200,  400,  600,  800,  and  950)  °C  and 
(2, 5, 10,  20, 40,  60)  s,  respectively.  Excitation  wavelength  of  488  nm  and  533.6  nm  were 
used.  488  nm  excites  Tb  and  Eu,  while  533.6  nm  only  excites  Eu.  Figure  77  shows  an  XRD 
spectra  and  two  fluorescence  spectra  of  p-Eu:Zr02/p-Tb:Y203  core/shell  nanoparticles 
excited  with  488  nm  and  533.6  nm,  respectively. 
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Figure  77.  XRD  spectra  (left),  and  fluorescence  spectra  of  p-Eu:Zr02/p-Tb:Y2C>3 
core/shell  nanoparticles  excited  with  488  nm  (center)  and  533.6  nm  (right), 
respectively. 

The  fluorescence  spectra  can  be  correlated  with  the  formations  of  the  tetragonal  Zr02  (t- 
ZrCX)  and  cubic  Y2O3  (C-Y2O3)  phases.  XRD  shows  that  the  C-Y2O3  phase  is  formed  at  800 
°C  and  higher  temperatures,  while  the  t-ZrCX  phase  is  formed  at  600  °C  and  higher 
temperatures.  Similarly,  the  fluorescence  spectra  show  that  Tb  is  located  in  a  crystalline 
t-Zr02  phase  starting  at  600  °C,  while  Eu  is  located  in  a  C-Y2O3  phase  starting  at  800  °C. 
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Figure  78.  Fluorescence  excited  with  488  nm  (left)  and  XRD  (right)  of  p-EmZrCX/p- 
Tb:Y2C>3  core/shell  nanoparticles  heated  to  600  °C  for  10  s  compared  with  those  of  a 
fully  crystalline  Tb:ZrC>2  sample. 


The  optical  and  XRD  spectra  of  the  heated  samples  can  then  be  compared  with  reference 
samples.  Figure  78  shows  the  data  for  p-Eu:Zr02/p-Tb:Y203  core/shell  nanoparticles 
heated  to  600  °C  for  10  s.  Figure  79  shows  the  fluorescence  spectra  and  XRD  spectrum  of 
p-Eu:ZrC>2/p-Tb:Y203  core/shell  nanoparticles  heated  to  800  °C  for  10  s  and  reference 
samples.  The  spectra  indicate  that  diffusion  across  the  interface  occurs  and  the  spectra 
can  be  assigned  to  the  Eu:t-Zr02  and  Eu:Y2C>3  phases. 
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Figure  79.  Fluorescence  spectra  (left)  and  XRD  spectrum  (right)  of  p-Eu:Zr02/p- 
Tb:Y203  core/shell  nanoparticles  heated  to  800  °C  for  10  s  and  reference  samples.  The 
spectra  can  be  assigned  to  the  Eu:t-Zr02  and  Eu:Y2C>3  phases. 

Using  the  fluorescence  spectra,  five  distinct  phases  can  be  identified  within  the  evaluated 
(T,t)  space.  However,  in  some  cases  (e.g.  regions  IV  and  V)  the  differences  are  not  as 
conclusive  as  in  other  cases.  Figure  80  shows  the  fluorescence  spectra,  exciting  Tb  and  Eu 
using  488  nm,  belonging  to  the  different  regions,  as  well  the  resulting  phase  diagram. 
Regions  I,  II,  III,  IV,  and  V  correspond  to  the  following  heat  treatments:  200  °C  and  400 
°C  for  (10,  20,  30,  45,  and  60)  s;  600  °C  for  (10,  20,  30,  45,  and  60)  s;  800  °C  for  (2  and 
5)  s;  800  °C  for  (10,  20,  30,  45,  and  60)  s  and  950  °C  for  (2  and  5)  s;  and  950  °C  for  (10, 
20,  30,  45,  and  60)  s,  respectively. 
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Figure  80.  Distinct  phases  in  (T,t)  space  and  resulting  phase  diagram  for  488  nm 
excitation. 
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Figure  81  shows  the  spectra,  exciting  Eu  using  533.6  nm,  belonging  to  the  different 
regions,  as  well  the  resulting  phase  diagram.  These  assignments  are  somewhat  less 
ambiguous  than  the  ones  resulting  from  excitation  at  488  nm.  Regions  I,  II,  III,  IV,  and 
V  correspond  to  the  following  heat  treatments:  200  °C  and  400  °C  for  (10, 20, 30, 45,  and 
60)  s;  600  °C  for  (2,  5, 10,  20,  30,  45,  and  60)  s;  800  °C  and  950  °C  for  (2  and  5)  s;  800 
°C  for  (10,  20,  30,  45,  and  60)  s;  and  950  °C  for  (10,  20,  30,  45,  and  60)  s,  respectively. 
Depending  on  the  specific  excitation  wavelength  chosen,  some  emission  from  the 
crystalline  Tb  compound  can  overlap  with  the  emission  from  the  precursor  Eu  compound. 
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Figure  81.  Distinct  phases  in  (T,t)  space  and  resulting  phase  diagram  for  533.6  nm 
excitation. 


p-Eu:ZrC>2  +  c-Ho:Zr02  +  p-Dy:Y2(>3 

To  minimize  spectral  overlap  issues,  the  following  mixture  has  been  chosen  as  the  most 
promising  compound,  even  though  it  requires  two  measurements  using  two  different 
excitation  wavelength:  p-Eu:Zr02  +  c-Ho:ZrC>2  +  p-Dy:Y203. 

First,  we  evaluate  a  mixture  of  c-Ho:ZrC>2  and  p-Dy:Y203.  The  emission  spectra  of  the 
unheated  materials  for  excitation  at  355  nm  are  shown  in  Figure  82  and  Figure  83.  Figure 
82  shows  the  spectra  of  the  individual  components.  While  there  is  some  spectral  overlap, 
significant  portions  of  the  peaks  do  not  suffer  from  overlap.  Figure  83  shows  the  spectra 
of  combined  mixtures  for  different  ratios.  Even  though  p-Dy:Y203  is  a  precursor  material, 
it  emission  intensity  is  very  strong  compared  to  the  crystalline  c-Ho:Zr02.  Figure  84 
shows  the  emission  spectra  of  heated  (950  °C)  c-Ho:Zr02  +  p-Dy:Y203  for  several  ratios. 
The  spectra  are  dominated  by  the  Dy  emission,  indicating  that  the  ratio  needs  to  be 
adjusted. 
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Figure  82.  Individual 
spectra  of  unheated  c- 
Ho:Zr02  and  p-Dy:Y203. 


Wavelength  (nm) 

Figure  83.  Spectra  of 
unheated  mixtures  of  c- 
Ho:Zr02  and  p-Dy:Y203. 
The  Ho  emission  is  much 
weaker  than  the  Dy 
emission. 


Figure  84.  Spectra  of 
heated  c-Ho:Zr02+p- 
Dy:Y203.  The  Dy  peak 
dominates  the  emission  for 
all  ratios  tested. 


Figure  85  shows  the  individual  spectra  of  unheated  Eu:Zr02  and  c-Ho:Zr02/p-Dy:Y203 
core/shell  particles  excited  at  355  nm,  and  Figure  86  shows  the  spectra  of  the  unheated 
mixture  for  various  ratios.  Figure  87  shows  the  spectra  of  heated  (950  °C)  mixtures  for 
various  ratios.  Again,  the  Dy  emission  dominates  the  spectra. 


Wavelength  (nm) 

Figure  85.  Individual 
spectra  of  unheated 
Eu:Zr02  and  c-Ho:Zr02/p- 
Dy:Y203. 
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Figure  86.  Spectra  of 
different  mixtures  of 
unheated  Eu:Zr02  and  c- 
Ho:Zr02/p-Dy:Y203. 
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Figure  87.  Spectra  of 
heated  mixture  of  Eu:Zr02 
and  c-Ho:Zr02/p-Dy:Y203. 
The  Dy  emission 
dominates  the  heated 
spectrum. 


However,  the  overall  shape  of  the  spectra  depends  on  the  excitation  wavelength.  Figure 
88  shows  the  spectra  of  heated  Eu:Zr02  +  c-Ho:Zr02/p-Dy:Y203  for  four  different 
excitation  wavelengths,  355  nm,  365  nm,  488  nm,  and  532  nm,  and  Figure  89  shows  a 
direct  comparison  of  the  emission  spectra  of  Eu:Zr02  +  c-Ho:Zr02  +  p-Dy:Y203  in  a  1:1:1 
ratio.  Considering  the  relative  emission  intensities  of  the  different  compounds  and 
potential  corresponding  spectral  overlap,  one  possible  combination  of  excitation 
wavelengths  is  355  nm  and  532  nm. 
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Figure  88.  Spectra  of  heated  mixtures  of  Eu:Zr02+c-Ho:Zr02/p-Dy:Y203  for 
excitation  at  355  nm  (top  left),  365  nm  (top  right),  488  nm  (bottom  left),  and  532  nm 
(bottom  right). 
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Figure  89.  Spectra  of  an  equal  ratio  of  Eu:Zr02+c-Ho:Zr02+p-Dy:Y203  for  different 
excitation  wavelengths.  For  excitation  using  532  nm  the  Dy  peak  is  unobserved  and  a 
residual  pump  peak  is  observed  even  with  the  notch  filter. 
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We  performed  an  initial  evaluation  of  the  thermal  impulse  sensor  combination  p-Eu:ZrC>2 
+  c-Ho:ZrC>2  +  p-Dy:Y203.  Using  indirect  laser  heating,  sample  material  was  heated  to 
various  temperatures  for  specific  isothermal  heating  durations.  The  fluorescence  spectra 
of  these  samples  were  then  measured  using  pulsed  355  nm  laser  pulses,  see  Figure  90. 
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Figure  90.  Fluorescence  spectra  of  p-Eu:Zr02  +  c-Ho:Zr02  +  p-Dy:Y203  heated  to 
400  °C  (top  left),  500  °C  (top  center),  600  °C  (top  right),  700  °C  (bottom  left),  800  °C 
(bottom  center),  and  900  °C  (bottom  right)  for  various  times. 


In  addition  to  using  Ho  peaks  for  computing  intensity  ratios,  we  also  calculate  ratios  using 
internal  peaks  (i.e.  different  peaks  from  the  same  material).  For  Dy:Y203  we  use  the  ratio 
between  the  intensity  at  573  nm  and  578  nm.  These  wavelengths  are  chosen  because 
amorphous  Dy:Y2C>3  is  found  to  have  a  broad  peak  centered  at  578  nm  and  crystalline 
Dy:Y203  is  found  to  have  a  sharp  peak  at  573  nm.  The  ratio  between  the  two  therefore 
corresponds  to  a  direct  measure  of  the  crystallinity  of  the  Dy:Y203.  Figure  91  shows  the 
573/578  ratio  as  a  function  of  isothermal  time  with  fits  to  a  simple  exponential.  From 
Figure  91  we  find  that  there  is  a  drastic  change  in  the  ratios  behavior  as  a  function  of  time 
between  700  C  and  800  C.  This  corresponds  to  the  temperature  range  in  which 
crystallization  can  occur. 
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Figure  91.  Ratio  of  fluorescence  intensities  at  573  nm  and  578  nm  as  a  function  of 
time. 

While  Dy:Y203  has  one  very  clear  and  clean  intensity  ratio  to  consider  for  measuring 
phase  changes,  Eu:Zr02  has  multiple  peaks  which  change  with  temperature  and  time 
leading  to  several  different  intensity  ratios  we  can  use  for  extracting  temperature  and 
time.  The  main  peaks  of  Eu:Zr02  we  can  use  for  this  purpose  are  at  592  nm,  606  nm,  613 
nm,  631  nm,  660  nm,  and  714  nm.  Additionally,  the  intensity  at  620  nm  drastically 
changes  during  crystallization  as  the  613  nm  peaks  FWHM  collapses.  Figure  92  shows  the 
592/613  ratio  as  a  function  of  time  for  various  temperatures. 
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Figure  92.  592/613  Ratio  as  a  function  of  time  for  different  temperatures. 

After  obtaining  the  ratio  curves  as  a  function  of  time  and  temperature  we  first  attempted 
to  fit  the  fluorescence  ratios  to  the  isothermal  degree  of  conversion  function  and 
discovered  that  the  single  phase  model  fails  to  adequately  describe  the  observed 
phenomenon.  Therefore  in  order  to  use  the  calibration  data  to  determine  temperature 
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and  time  we  instead  use  a  phenomenological  model.  For  our  phenomenological  model  we 
define  the  ratios  as  a  function  of  time  t,  to  follow  an  exponential  function: 

R(T,  t)  =  Rm(T)  -  [Rm(T)  -  J?0(O]e~flCr)t  (  42  ) 

where  g(T)  is  the  exponential  rate  as  a  function  of  temperature,  Rcc(T)  is  the  asymptotic 
ratio  as  a  function  of  temperature,  and  Ro(T)  is  the  ratio  offset  as  a  function  of 
temperature.  The  ratio  offset  is  related  to  the  effects  of  slow  heating  rates. 

Using  Equation  42,  we  fit  the  573/541  ratio  and  592/613  ratio  as  a  function  of  time  for 
each  temperature  tested.  Figure  93  shows  the  fit  parameters  as  a  function  of  temperature 
for  the  573/541  ratio  and  Figure  94  shows  the  fit  parameters  as  a  function  of  temperature 
for  the  592/613  ratio.  The  parameters  for  592/613  behave  as  expected,  with  each 
parameter  increasing  with  temperature.  However,  for  the  573/541  ratio  the  exponential 
rate  behaves  opposite  to  what  is  expected  with  the  rate  decreasing  with  temperature.  The 
nature  of  this  result  is  currently  unknown. 

While  the  basic  behavior  of  the  fit  parameters  as  a  function  of  time  is  as  expected,  the 
precise  functional  form  is  found  to  be  different  than  expected.  The  parameters  (aside  from 
the  573/541  rate)  are  found  to  follow  a  modified  Arrhenius  function  given  by: 

/CO  =  fo  +  A  exp  j-  (~)^j  ( 43 ) 

where  the  stretch  parameter  [3  is  found  to  be  greater  than  1.  For  fitting  the  573/541  rate 
we  use  a  simple  power  function. 
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Figure  93.  Ratio  versus  time  fit  parameters  as  a  function  of  temperature  for  the 
573/541  ratio. 
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Figure  94.  Ratio  versus  time  fit  parameters  plotted  as  a  function  of  temperature  for 
the  592/613  ratio.  All  three  parameters  are  found  to  follow  a  modified  Arrhenius  fit 
function. 


Using  the  parameter  curves  as  a  function  of  temperature  -  determined  from  fitting  the 
data  shown  in  Figure  93  and  Figure  94  -  and  using  Equation  42,  we  can  determine  a 
functional  mapping  of  the  ratios  as  a  function  of  temperature  and  time.  Figure  95  shows 
the  ratio  maps  as  a  function  of  temperature  and  time  for  the  573/541  ratio  and  592/613 
ratio. 
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Figure  95.  Ratio  maps  as  a  function  of  time  and  temperature  for  two  573/541  (left) 
and  592/613  (right)  intensity  ratios. 
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The  current  measurements  have  relatively  large  uncertainties.  A  variety  of  factors  are 
responsible: 

1.  The  laser  profile  is  spatially  and  possibly  temporally  inhomogeneous,  showing 
various  hot  spots.  These  hot  spots  lead  to  inhomogeneous  heating  of  the  samples 
during  direct  laser  heating. 

2.  Due  to  the  inhomogeneous  laser  profile,  we  are  currently  working  with  an  indirect 
laser  heating  method.  This  method  relies  on  heat-conductivity  (graphite  sheets)  to 
spread  the  heat  throughout  the  sample.  The  finite  heat  conductivity  leads  to 
temperature  gradients  across  the  graphite  sheet  and  the  sample  material. 

3.  The  indirect  laser  heating  method  also  limits  the  heating  rate  to  300  °C  to  400  °C 
(i.e.  a  sample  heated  to  800  °C  for  100  ms  also  experienced  2  s  of  ramp-up  heating). 

4.  There  is  some  uncertainty  due  to  sample  handling  during  the  heating  procedure. 

5.  Shot-to-shot  noise  during  spectral  measurements. 

The  root  causes  of  these  uncertainties  will  be  addressed  next  to  improve  the  accuracy  of 
the  combined  temperature  and  heating  duration  measurements.  Meanwhile,  an  initial 
test  using  five  blind  samples  yielded  the  following  results,  see  Table  10.  Addressing  the 
above  listed  sources  of  uncertainty,  the  results  are  expected  to  significantly  improve. 


Table  10.  Analysis  result  for  five  blind  samples. 


Heating  temperature 
and  duration 

Measured  heating 
temperature  and  duration. 

Sample 

T  (K) 

t  (ms) 

T  (K) 

t  (ms) 

1 

673 

1003 

623  -  723 

100  -  6000 

2 

1173 

114 

1181  - 1228 

350  -  790 

3 

873 

303 

888  -  958 

100  -  400 

4 

unheated 

unheated 

5 

1073 

605 

1077  - 1159 

13  - 185 
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E.  Tracking  Sensors 

Since  heating  in  an  explosive  environment  is  expected  to  be  spatially  inhomogeneous,  it 
is  important  to  know  the  provenance  of  the  temperature  sensors.  As  such,  we  prepared 
and  evaluated  different  tracking  particles  that  can  be  mix  in  with  the  temperature  sensors 
to  provide  information  about  their  original  location.  These  tracking  particles  consist  of 
fully  crystallized  rare-earth  doped  yttria,  including  Eu:Y203,  Tb:Y203,  Dy:Y203,  Sm:Y203, 
Er,Yb:Y2C>3,  Ho:Y203,  Pr:Y2C>3,  Nd:Y203,  and  Tm:Y2C>3.  After  preparing  these  materials, 
their  fluorescence  spectra  were  measured  and  compared  with  those  of  the  temperature 
sensor  materials  to  make  sure  the  spectral  emission  does  not  interfere,  see  Figure  96. 


Figure  96.  Fluorescence  spectra  of  potential  tracking  particles  for  excitation  with  365 
nm  (left)  and  488  nm  (right). 


As  can  be  seen  from  Figure  96,  excitation  with  different  wavelengths  leads  to  significantly 
different  fluorescence  spectra.  Based  upon  our  analysis,  we  selected  Dy:Y203,  Sm:Y203, 


and  Er,Yb:Y203  for  initial  explosion  testing. 
The  layout  for  the  explosion  testing  is  shown 
in  Figure  97.  A  charge  was  placed  in  the 
center  of  the  room  and  1.5  g  each  of  the  three 
tracking  particles  was  wrapped  in  a  paper 
bag  and  hung  from  the  ceiling  at  distances  of 
2’,  4’,  and  6’,  respectively.  Post-detonation, 
debris  samples  were  collected  in  the  areas 
marked  a  through  i.  Also  shown  in  the  figure 
are  the  tracking  particles  detected  in  those 
specific  area.  No  tracking  particles  were 
detected  in  area  c.  Dy:Y203  was  detected  in 
all  other  areas,  which  is  probably  due  to  its 
very  strong  fluorescence  compared  to  the 
other  samples.  Sm:Y203  was  detected  only  in 
areas  a  and  d.  Its  fluorescence  intensity  is 
significantly  weaker  than  that  of  Dy:Y203. 
Er,Yb:Y203  was  detected  in  five  of  the  areas. 


Bombproof  Door 

Figure  97.  Layout  for  explosion  test  of 
tracking  particles. 
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Bombproof  Door 

Figure  98.  Layout  for  explosion  test  of 
tracking  particles  combined  with  p- 
Eu:Y203/c-Tb:Y203  temperature  sensors. 
We  were  unable  to  measure  a  response  from 
this  sensor  material. 


Bombproof  Door 


Figure  99.  Layout  for  explosion  test  of 
tracking  particles  combined  with  p- 
Eu:Zr02  temperature  sensors.  The 
measured  temperature  ranges  (°C)  are 
indicated  in  the  layout. 


The  quasi-static  temperature,  measured  via  a  thermocouple,  for  shot  2  at  a  distance  of  6’ 
was  430  °C.  This  corresponds  to  a  difference  of  about  17%  to  24%  between  the 
thermocouple  and  the  temperature  sensor  measurements,  with  the  temperature  sensors 
measuring  the  higher  values. 
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V.  CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  this  investigation: 

1.  Temperature  Sensors 

Modeling  of  the  phase  transitions  confirms  that  the  effect  of  time  is  significantly 
less  than  the  effect  of  temperature  and  that  therefore  calibration  measurements  at 
longer  times  can  be  used  to  extract  temperature  information.  The  following 
materials  have  been  successfully  demonstrated: 

a.  p-Eu:Y203 

This  material  works  well  as  a  temperature  sensors  in  explosive  fireballs. 
Various  spectral  signatures  can  be  used  to  determine  the  temperature. 
Below  the  crystallization  temperature,  the  FWHM  and  peak  wavelength  of 
the  fluorescence-excitation  peak  of  the  ?Fo  5Do  transition  are  good 
indicators  for  temperature.  Above  the  crystallization  temperature,  the 
FWHM  of  the  sDo  ->  ?F0  fluorescence  transition  is  a  good  indicator  for 
temperature.  In  addition,  the  intensity  ratio  for  the  sDo  ->  ?F2  vs.  sDo  ->  ?Fi 
fluorescence  peaks  can  be  used  as  an  indicator.  Furthermore,  peak 
intensities  with  respect  to  reference  materials  can  be  used  to  determine 
temperature. 

b.  p-Eu:Zr02 

This  material  also  works  very  well  as  a  temperature  sensor  in  explosive 
fireballs.  High  spectral  resolution  measurements  show  that  the  peak 
intensity  ratio  and  peak  position  of  the  split  sD0  ->  ?Fi  fluorescence  can  be 
used  as  indicators  for  temperature.  In  addition,  the  ratio  of  the  tetragonal 
to  monoclinic  phases,  as  well  as  the  intensity  ratio  compared  with  reference 
materials  can  be  used  to  determine  temperature. 

In  general,  fluorescence  measurements  are  easier  and  faster  to  perform  than 
fluorescence-excitation  measurements,  an  important  parameter  for  future 
portable  equipment  suitable  for  field  operations. 

2.  Thermal  impulse  sensors 

Thermal  impulse  sensors  consist  of  a  mixture  of  two  sensor  materials  to  be  able  to 
simultaneously  determine  temperature  and  time.  The  following  material  has  been 
successfully  demonstrated  under  laboratory  conditions: 

a.  p-Eu:Zr02  +  c-Ho:Zr02  +  p-Dy:Y203 

This  “cocktail”  consist  of  the  already  demonstrated  host  materials  Zr02  and 
Y203.  While  Zr02  is  still  doped  with  Eu,  the  dopant  for  Y203  had  to  be 
changed  to  avoid  spectral  overlap.  As  such,  Dy  was  chosen  as  it  shows  strong 
emission  lines  and  temperature-induced  spectral  changes.  In  addition,  c- 
Ho:Zr02  was  added  as  a  reference.  Using  this  material  combination,  we 
were  able  to  extract  temperatures  between  400  °C  and  900  °C  for  heating 
durations  as  short  as  100  ms.  Shorter  heating  durations  are  possible  but 
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require  improved  reference  measurements.  Also,  further  optimizations  of 
the  sensor  composition  are  expected. 

3.  Tracking  sensors 

Of  the  currently  tested  tracking  sensors,  Dy:Y203  appears  to  be  the  most  promising 
one  as  is  shows  very  strong  emission  intensities.  However,  this  material  cannot  be 
used  in  combination  with  the  current  thermal  impulse  sensor  material 
combinations  since  it  contains  the  same  material.  Er,Yb:Y203  is  the  second  best 
material  and  does  not  interfere  with  any  of  the  other  materials.  Sm:Y203  has  the 
weakest  emission  intensity  of  the  tested  materials.  Also,  depending  on  the  selection 
of  the  thermosensor  host  material,  it  may  interfere  with  Eu  from  the  thermal 
sensors. 


We  have  successfully  demonstrated  that  small  luminescent  particles  can  be  used  as 
temperature  sensors  during  explosions.  As  expected,  the  readings  from  these  sensors 
typically  show  higher  temperatures  than  the  thermocouples.  Also,  our  sensors  show  a 
range  of  temperatures,  reflecting  the  range  of  temperatures  within  an  explosive  fireball, 
in  particular  when  reactive  metals  such  as  aluminum  are  present.  Initial  tests  of  thermal 
impulse  sensors  and  tracking  sensors  also  look  promising. 
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VI.  RECOMMENDATIONS 


Given  the  demonstrated  functionality  of  the  temperature  sensors  and  the  progress  made 
today,  it  is  clear  that  the  approach  works.  However,  further  refinements  can  be  made  to 
the  overall  design  and  functionality  of  the  various  sensors.  In  particular,  the  following  is 
recommended: 

1.  Refine  laboratory-based  short-pulse  laser  heating  system. 

2.  Optimize  thermal  impulse  (TI)  sensors  for  heating  durations  as  short  as  l  ms. 

3.  Optimize  tracking  sensors  to  minimize  spectral  interference  with  temperature 
sensors. 

4.  Develop  portable  analysis  equipment. 
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IX.  APPENDIX 


A.  Alternative  Temperature  Sensors 
Eu:La203  -  based  materials 

Monodispersed  spherical  p-Eu:La203  nanoparticles  are  obtained  by  heating  an  aqueous 
mixture  of  La(N03)3.6H20,  Eu(N03)3-6H20,  citric  acid,  and  urea  with  salt/citric 
acid/urea  ratio  of  1/0.25/40  in  a  90  °C  thermostated  oil  bath  for  3  h.  The  total  metallic 
salt  concentration  is  0.005  M.  Likewise,  monodispersed  spherical  nanoparticles  of  p- 
Tb:Y203  are  obtained  via  homogeneous  precipitation.  An  aqueous  mixture  of 
Yb(N03)3-6H20,  Tb(N03)3.6H20,  and  urea  with  salt/urea  ratio  of  1/40  is  heated  under 
continuous  stirring  in  a  90  °C  thermostated  oil  bath.  The  total  metallic  salt  concentration 
is  0.020  M. 

Core/shell  p-Eu:La203/p-Tb:Y203  nanoparticles  are  prepared  by  repeating  the 
homogeneous  precipitation  of  p-Tb:Y203  three  times  in  the  presence  of  p-Eu:La203 
nanoparticles.  The  salt  concentration  of  the  metallic  salt  is  0.005  M  and  the  salt  to  urea 
ratio  is  1/40.  Likewise,  the  core/shell  p-Tb:Y203/p-Eu:La203  nanoparticles  are  prepared 
by  repeating  the  homogeneous  precipitation  of  p-Eu:La203  in  the  presence  of  citric  acid 
three  times  in  the  presence  of  p-Tb:Y203  nanoparticles.  The  salt  concentration  of  the 
metallic  salt  is  0.005  M  and  the  salt/citric  acid/urea  to  urea  ratio  is  1/0.25/40. 

Figure  100  shows  SEM  images  of  p-Eu:La203  prepared  by  homogeneous  precipitation 
using  variable  amounts  of  urea.  The  anisotropy  of  the  resulting  structures  increases  with 
increasing  urea  concentration.  A  salt/urea  ratio  of  1/1,  results  in  the  formation  of  2-D 
sheets,  while  salt/urea  ratios  of  1/40  and  1/200  result  in  the  formation  of  l-D  structures 
with  similar  morphologies. 


Figure  100.  SEM  micrographs  of  p-Eu:La203  synthesized  with  (a,  d)  1/1,  (b,  e)  1/40,  and 
(c,  f)  1/200  salt/urea  ratios.  In  all  cases  the  salt  concentration  is  20  mM. 

Next,  we  perform  a  series  of  homogeneous  precipitation  of  p-Eu:La203  in  the  presence  of 
citric  acid  as  a  capping  agent  at  various  salt  to  citric  acid  ratios,  but  keep  the  salt/urea 
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ratio  fixed  at  1/40.  Figure  101  shows  spherical  nanostructures  formed  at  a  1/0.05 
salt/citric  acid  ratio,  and  a  bimodal  nanoparticle  distribution  at  a  1/0.50  salt/citric  acid 
ratio.  When  the  salt/citric  acid  ratio  is  1/2,  the  solution  remains  transparent,  indicating 
no  precipitation  has  taken  place,  after  3  h  of  heating  at  90  °C  (SEM  micrographs  not 
shown). 


(a)  Kb) 


Figure  101.  SEM  micrographs  of  p-Eu:La203  synthesized  with  (a) 
1/0.05,  (b)  1/0.10,  (c)  1/0.25,  and  (d)  1/0.50  salt/citric  acid  ratio.  In  all 
cases  the  salt/urea  ratio  is  1/40.  The  salt  concentration  is  20  mM. 


By  comparing  Figure  100  (b)  and  Figure  100  (e)  with  Figure  101,  we  see  that  the  addition 
of  citric  acid  reduces  the  formation  of  anisotropic  structures.  In  the  presence  of  citric  acid, 
lanthanum-citric  acid  complexes  form.  The  coordinated  citric  acid  molecules  suppress  the 
preferential  growth  direction  during  precipitation.51  In  the  case  of  the  1/0.25  salt/citric 
acid  ratio  the  nanoparticles  are  more  or  less  monodispersed  and  segregated,  see  Figure 
101  (c).  However,  micrometer-sized  aggregates  can  be  seen  as  well.  It  is  likely  that  the 
large  particles  form  as  a  result  of  Ostwald  ripening  with  adjacent  nanoparticles  fusing  into 
larger  particles.  To  further  improve  the  monodispersity  of  the  spherical  nanoparticles,  we 
vary  the  salt  concentration.  Figure  102  shows  the  quality  of  the  nanoparticles  at  various 
salt  concentrations.  The  salt/citric  acid/urea  ratio  is  fixed  at  1/0.25/40.  The  extent  of 
aggregation  increases  with  increasing  salt  concentration.  Relatively  monodisperse 
spherical  particles  with  an  average  diameter  of  (375  ±  60)  nm  form  at  a  salt  concentration 
of  5  mM,  see  Figure  102  (a).  All  further  discussion  is  limited  to  p-Eu:La203  nanoparticles 
synthesized  with  a  salt/citric  acid/urea  ratio  of  1/0.25/40  and  a  salt  concentration  of  5 
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mM.  In  addition,  precursor  Tb:Y203  (p-Tb:Y20s)  is  also  prepared  by  means  of 
homogeneous  precipitation,  but  without  citric  acid. 


Figure  102.  SEM  micrographs  of  p-Eu:La203  synthesized  using  various  salt 
concentrations:  (a)  5  mM  (b)  its  corresponding  size  distributions,  (c)  20  mM,  and 
(d)  80  mM.  The  salt/citric  acid/urea  ratio  is  fixed  at  1/0.25/40. 


Figure  103  shows  a  simultaneous  TGA/DSC  measurement  of  p-Eu:La203  and  p-Eu:Y203. 
The  TGA  graphs  show  a  mass  loss  of  34%  and  35%  for  p-Eu:La203  and  p-Tb:Y203, 
respectively.  These  losses  can  be  ascribed  to  the  endothermic  decomposition  of 
La0HC03-H20  and  YOHCO3H2O.  MOHCO3H2O  compounds  undergo  multi-stage 
decomposition  towards  M2O3,  which  proceeds  as  follows  ;52>  53 


2  MOHCO3H2O  -»•  M2O2CO3  +  3  H2O  +  CO2  -»•  M2O3  +  CO2 
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(a)  (b) 


Figure  103.  DSC/TGA  graph  of  (a)  p-Eu:La203  and  (b)  p-Tb:Y203  nanoparticles. 


Table  11  compares  the  calculated  mass  loss  with  those  derived  from  TG  measurement 
while  assuming  that  the  starting  compound  is  M0HC03-H20.  The  experimental  results 
indicate  good  agreement  with  the  proposed  decomposition  channel. 


Table  11.  Calculated  and  measured  mass  loss  assuming  that  the  starting  compound  is 
MQHC03-H20. _ 


Assumed 

Calculated 

Measured 

Assumed 

Calculated 

Measured 

compound 

mass  /  % 

mass  /  % 

compound 

mass  /  % 

mass  /  % 

La0HC03-H20 

100 

100 

YOHCO3  -H20 

100 

100 

La202C03 

79 

74 

y2o2co3 

73 

74 

La203 

70 

66 

y2o3 

6l 

65 

The  p-Eu:La203  nanoparticles  are  heated  to  temperatures  between  300  °C  and  900  °C  for 
3  h.  The  corresponding  x-ray  diffraction  patterns  show  crystalline  phase  formation  in  the 
samples  heat  treated  >  500  °C/3  h.  The  crystalline  phases  of  the  samples  heat  treated  at 
500  °C/3  h,  600  °C/3  h,  and  700  °C  for  3  h,  can  be  assigned  to  m-La202C03  (PDF  00- 
048-1113),  h-La(OH)3  (PDF  01-083-2034),  andh-La203  (PDF  01-074-1144),  respectively, 
see  Figure  104  (a).  All  of  the  carbonate  phases  of  samples  heat  treated  at  900  °C/3  h  are 
expected  to  have  transformed  into  the  h-La203  phase.  However,  this  h-La203  phase 
readily  hydrolyzes  into  the  La(OH)3  phase  at  ambient  conditions  in  air.54  As  a  result,  and 
depending  on  its  humidity  exposure  history,  varying  degrees  of  mixed  La203/La(0H)3 
phases  may  be  observed. 

The  susceptibility  of  the  h-La203  phase  to  moisture  absorption  is  widely  reported  in  the 
literature.54-56  Annealing  may  be  used  to  transform  the  hydrated  La203  species  back  into 
the  La203  phase.  Also,  a  protective  layer  may  be  incorporated  to  minimize  hydration  in 
air.s6  Figure  104  (a)  shows  that  the  phase  of  the  sample  heat  treated  at  500  °C/3  h  is 
mainly  m-La202C03.  Samples  heat  treated  at  >  600  °C/3  h  consist  of  mixtures  of  the  m- 
La202C03,  h-La(OH)3,  and  La203  phases.  The  majority  of  the  sample  heat  treated  at  600 
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°C/3  h  is  m-L^CXCOs,  while  the  sample  heat  treated  at  900  °C/3  h  consists  of  the  h- 
La(0H)3  and  La203  phases.  The  sample  heat  treated  at  700  °C/3  h  contains  all  three 
phases. 
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1  1  1  T  1  ’  1  1  1  1 

10  20  30  40  50  60  70 


2  Theta  /  degree 


540  560  540  €00  620  640 
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Figure  104.  (a,  c)  XRD  and  (b,  d)  fluorescence  of  (a,  b)  p-Eu:La203  and  (c,  d)  p- 
Tb:Y2C>3  heated  to:  (Uncalcined,  300,  500,  600,  700,  900)  °C/3  h.  Markers  in  (a)  and 
(c)  correspond  to:  m-La202C03  (black),  h-La(OH)3  (magenta),  h-La203  (blue),  and  c- 
Y2O3  phase  (gray). 


Figure  104  (b)  shows  the  fluorescence  spectra  of  the  p-Eu:La203  samples  with  peaks 
corresponding  to  the  Eu3+  5D0  — >  ?Fj  (J=o-3)  transitions,  and  markers  for  the  positions  of 
the  Eu3+  free  ion  transitions.  The  emission  spectra  reflect  the  site  symmetry  of  the  Eu3+ 
ions  in  the  various  phases.  The  inhomogeneous  peak  broadening  reflects  the  amorphous 
nature  of  the  p-Eu:La203  samples  in  which  the  Eu3+  ions  are  surrounded  by  disordered 
ligands.  Following  heat  treatments  at  300  °C/3  h  and  500  °C/3  h,  the  color  of  the  p- 
Eu:La203  samples  change  from  white  to  black  due  to  the  decomposition  of  the  p-Eu:La203 
and  the  formation  of  carbon  ash.  No  emission  peaks  characteristics  of  Eu3+  ions  are 
observed  for  the  sample  heat  treated  at  500  °C/3  h  because  the  ash  absorbs  most  of  the 
excitation  and  fluorescence  light.  As  suggested  by  the  formation  of  narrow  diffraction 
peaks  in  their  XRD  spectra,  the  samples  heat  treated  at  or  above  6oo°C/3  h  have 
transformed  into  a  crystalline  phases  such  as  Eu:La202C03,  Eu:La(OH)3,  and  Eu:La203.52 
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The  spectra  of  the  sample  heat  treated  at  900  °C/3  h  shows  strong  emission  in  the  sDo  — > 
7Fi  transitions  compared  to  the  sD0  — ►  7F2  transitions,  indicating  a  high  degree  of  spherical 
site  symmetry.  According  to  the  powder  diffraction  file,  the  site  symmetry  of  the  La3+ 
cations  in  a  h-La203  phase  is  C3  and  the  site  symmetry  of  the  La3+  cations  in  a  h-La(OH)3 
phase  is  C3I1.  By  assuming  that  the  Eu3+  cations  replace  the  La3+  host  cations  and  possess 
the  same  site  symmetries,  the  number  of  energy  levels  of  the  7Fj  may  be  derived  in  terms 
of  the  irreducible  representations,  see  Table  12.57 


Table  12.  Irreducible  representations  for  7Fj  levels  in  La2C>3  and  La(OH)3  . 


J 

C3  site 

C3h  site 

0 

Ti 

Ti 

1 

n  +  r2  + 

Ti  +  r5,6 

2 

Ti  +  2r2  +  2T3 

Ti  +  r2,3  +  T5,6 

Figure  104  (c)  and  Figure  104  (d)  show  the  x-ray  diffraction  patterns  and  fluorescence 
spectra  of  p-Tb:Y203,  respectively,  following  heat  treatments  between  300  °C/3  h  and  900 
°C/3  h.  Figure  104  (c)  shows  that  the  crystalline  Y203  phase  is  observed  only  for  the 
samples  heat  treated  at  700  °C/3  h  and  900  °C/3  h.  Their  XRD  patterns  can  be  matched 
with  the  cubic  Y2O3  phase  (PDF  01-083-0927).  No  diffraction  peaks  other  than  those  due 
to  the  cubic  Y2O3  phase  are  observed.  According  to  the  powder  diffraction  file,  the  site 
symmetries  of  the  Y3+  cations  are  C2  and  S6=C3i.  In  rare-earth  doped  Y2O3,  the  rare  earth 
ion  is  known  to  substitute  for  the  Y3+  cations.s8  By  assuming  that  the  Tbs+  cations  replace 
the  Y3+  host  cations  and  maintain  the  same  site  symmetries,  the  number  of  the  7Fj  energy 
levels  may  be  derived  in  terms  of  the  irreducible  representations,  see  Table  13.57 


Table  13.  Irreducible  representations  for  7Fj  levels  in  Y2O3. 


j 

C2  site 

C3i  site 

5 

5Ti  +  6T2 

3Ti  +  4T2  +  4I3 

4 

5Ti  +  4r2 

3Bi  +  3Bs  +  31^3 

3 

3Bi  +  4B2 

3Ti  +  2T2  +  2^ 

Figure  105  shows  the  SEM  micrographs  of  polycrystalline  Eu:La203  and  Tb:Y203 
nanoparticles  heat  treated  between  600  °C/3  h  and  900  °C/3  h.  The  Tb:Y203 
nanoparticles  are  segregated  and  remain  spherical  when  heated.  However,  due  to  necking 
and  sintering,  the  Eu:La203  nanoparticles  gradually  form  hard  aggregates  as  the 
temperature  increases  and  they  lose  their  spherical  morphology  when  heat  treated  at  900 
°C/3  h.  This  difference  is  due  to  the  lower  crystallization  temperature  of  the  p-La203 
phase  compared  to  the  p-Y203  phase,  see  Figure  104  (a)  and  Figure  104  (c). 
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Figure  105.  SEM  micrographs  of  the  polycrystalline  Eu:La203  (a,  b,  c)  and 
polycrystalline  Tb:Y2C>3  (d,  e,  f)  nanoparticles  heat  treated  at  (a,  d)  600  °C/3  h,  (b,  e) 
700  °C/3  h,  and  (c,  f)  900  °C/3  h. 

Oxide  nanoparticle  precursors  synthesized  via  homogeneous  precipitation  can  be  grown 
over  core  templates^8-  59  To  assess  the  versatility  of  the  modified  homogeneous 
precipitation  method  we  perform  our  synthesis  procedure  in  the  presence  of  Si02 
nanoparticles  (5  mg/mL)  using  salt/citric  acid/urea  ratios  of  1/0.25/40  and  1/0/40. 
Figure  106  compares  the  results  of  the  homogeneous  precipitation  of  a  p-Eu:La203  shell 
in  the  presence  and  in  the  absence  of  citric  acid.  Figure  106  (a)  shows  that  following  the 
homogeneous  precipitation  in  the  presence  of  citric  acid  some  of  the  Si02  nanoparticles 
exhibit  rougher  surface  morphologies  than  others.  In  the  absence  of  citric  acid,  a  mixture 
of  one-dimensional  p-La203  and  spherical  (Si02)  nanoparticles  is  obtained  (see  Figure 
106).  Figure  102  shows  that  increasing  the  concentration  of  the  La(N03)3-6H20  salt 
results  in  particle  aggregation.  In  order  to  ensure  that  all  of  the  core  nanoparticles  are 
covered  with  a  p-La203  shell  the  homogeneous  precipitation  is  always  repeated  three 
times. 
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Figure  106.  SEM  micrographs  of  p-Eu:La2C>3  synthesized  (a)  in  the  presence  and  (b) 
in  the  absence  of  citric  acid.  Respectively,  the  salt/citric  acid/urea  ratio  are  1/0.25/40 
and  1/0/40. 


To  further  demonstrate  the  usefulness  of  the  modified  homogeneous  precipitation 
technique,  we  prepare  Tb:La203/Eu:Y203  and  Tb:Y203/Eu:La2C>3  core/shell 
nanoparticles.  Figure  107(a)  and  Figure  107  (d)  show  that  the  p-Tb:Y203  and  p-Eu:La203 
cores  remain  spherical  during  the  formation  of  the  respective  p-Eu:La203  and  p-Tb:Y203 
shells  via  homogeneous  precipitation.  The  core-shell  morphologies  can  be  clearly  seen  in 
the  TEM  micrographs,  see  Figure  107  (b)  and  Figure  107  (e).  The  EDS  spectra  of  the 
Tb:La2C>3/Eu:Y203  and  Tb:Y203/Eu:La2C>3  core/shell  nanoparticles  confirm  the  presence 
of  Y  and  La,  see  Figure  107  (c)  and  Figure  107  (f). 


Figure  107.  SEM  and  TEM  micrographs  of  (a,  b)  p-Eu:La203/p-  Tb:Y203  and  (c,  d)  p- 
Tb:Y203/p-Eu:La203  core/shell  nanoparticles,  and  (c,  f)  corresponding  EDS  spectra. 


Figure  108  (a)  and  Figure  108  (c)  show  the  x-ray  diffraction  patterns  of  heat  treated 
(between  600  °C/3  h  and  900  °C/3  h)  p-Eu:La203/p-Tb:Y2C>3  and  p-Tb:Y203/p-Eu:La2C>3 
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core/shell  nanoparticles,  respectively.  The  peaks  can  be  assigned  to  the  m-La2C>2C03,  h- 
La(OH)3,  h-La203,  and  cubic  Y2O3  phases,  compare  with  Figure  104.  The  peaks  due  to  the 
La203  related  phases  are  more  prominent  than  the  Y2O3  phase  for  p-Eu:La203/p-Tb:Y203 
nanoparticles.  The  opposite  is  true  for  the  p-Tb:Y203/p-Eu:La203  core/shell 
nanoparticles.  These  results  indicate  the  dominance  of  the  core. 

As  can  be  seen  from  Figure  108  (b)  and  Figure  108  (d),  the  fluorescence  (Aex=488  nm) 
from  the  uncalcined  p-Eu:La20s/p-Tb:Y203  and  from  the  uncalcined  p-Tb:Y20s/p- 
Eu:La203  core/shell  nanoparticles  is  dominated  by  emission  from  Tb3+.  Simultaneous 
emission  due  to  Tb3+  and  Eu3+  is  seen  for  the  p-  Eu:La203/p-Tb:Y203  and  the  p- 
Tb:Y203/p-Eu:La203  core/shell  nanoparticles  heat  treated  to  (600,  700,  900)  °C/3  h. 
Even  though  the  nanopowder  samples  heat  treated  at  500  °C/3  h  have  a  black 
appearance,  some  Tb3+  emission  is  observed.  While  the  emission  from  the  shell 
components  is  clearly  visible,  the  fluorescence  spectra  of  the  samples  heat  treated  at  900 
°C/3  h  are  dominated  by  the  emission  from  the  core  components,  indicating  that  the 
crystallization  of  the  shell  is  not  as  efficient  as  the  crystallization  of  the  core. 


2  Theta  /  degree 


Figure  108.  (a,  c)  XRD  and  (b,  d)  fluorescence  of  p-Eu:La2C>3/p-Tb:Y203  (a,  b)  and  p- 
Tb:Y203/p-Eu:La203  (c,  d)  heated  to:  (uncalcined,  300,  500,  600,  700,  and  900)  °C/3 
h.  Markers  in  (a)  and  (c)  correspond  to  p-La203  (gray)  and  the  P-Y2O3  phase  (black). 
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Figure  109  shows  the  SEM  and  TEM  micrographs  of  the  polycrystalline 
Eu:La203/Tb:Y203  and  Tb:Y203/Eu:La203  nanoparticles  heat  treated  between  600  °C/3 
h  and  900  °C/3  h.  In  all  cases,  the  polycrystalline  Eu:La203/Tb:Y203  nanoparticles 
remain  segregated,  whereas  the  Tb:Y203/Eu:La203  nanoparticles  show  neckings  and 
sintering  following  900  °C/3  h  heat  treatment.  This  is  reminiscent  of  the  heat  treated  p- 
La2C>3  and  P-Y2O3,  compare  with  Figure  105.  That  is,  the  crystalline  phases  in  the  p-La203 
form  much  earlier  than  the  crystalline  phase  in  the  P-Y2O3,  see  also  see  Figure  104. 


Figure  109.  TEM  micrographs  of  the  polycrystalline  Eu:La203/Tb:Y203  (a,  b,  c)  and 
Tb:Y203/Eu:La203  (d,  e,  f)  nanoparticles  heat  treated  at  600  °C/3  h  (left),  700  °C/3  h 
(center),  and  900  °C/3  h  (right). 


Homogeneous  precipitation  of  nanophase  p-Eu:La203  using  La(N03)3-6H20  as  the  salt 
precursor  resulted  in  diverse  morphologies,  including  sheets  (2-D),  rods  (l-D),  and 
spheres  (o-D).  Homogeneous  precipitation  in  the  presence  of  citric  acid  results  in 
monodisperse  spherical  p-Eu:La203  nanoparticles.  The  optimum  salt/citric  acid/urea 
ratio  is  1/0.25/40  with  5  mM  salt  concentration.  We  have  also  demonstrated  that  the 
method  can  be  used  to  synthesize  core/shell  morphologies  with  p-Eu:La203  as  either  the 
core  (e.g.  p-Eu:La20s/p-Tb:Y203)  or  the  shell  (e.g.  Si02/p-Eu:La203  and  p-Tb:Y203/p- 
Eu:La203).  Heat  treatments  for  3  h  for  various  temperatures  leads  to  the  formation  of 
oxide  phases.  XRD  and  fluorescence  data  for  the  core/shell  samples  heat  treated  at  900 
°C/3  h  indicate  the  dominance  of  the  core.  The  reported  method  of  synthesizing 
monodisperse  spherical  shaped  nanophosphor  of  p-Eu:La2C>3  is  simple  and  low-cost.  It 
requires  a  low  heating  temperature  (90  °C)  for  3  h  in  an  aqueous  solution  in  the  presence 
of  citric  acid  and  urea.  The  h-La203  phase  is  obtained  following  heat  treatment  of  900 
°C/3  h.  Due  to  the  hygroscopic  nature  of  the  h-La203  phase,  it  readily  transforms  into  the 
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h-La(0H)3  phase.  Exposure  to  moisture  should  be  minimized  to  preserve  the  h-La203 
phase.  The  variety  of  different  phases  make  this  material  attractive  as  another  sensor 
material. 

B.  Eu:Ti02  -  based  materials 

Initial  experiments  were  performed  using  Eu:TiC>2.  Figure  no  shows  the  fluorescence 
spectra  excited  at  464  nm  and  normalized  for  the  emission  peak  at  592  nm  and  the  peak 
intensity  ratios  as  a  function  of  temperature.  The  emission  intensity  at  612.5  nm  decreases 
with  increasing  temperature  indicating  the  increasing  symmetry  of  the  Eu  site.  This  trend 
extends  to  temperatures  at  least  as  low  as  80  °C,  making  this  material  promising  as  a 
wide-temperature  range  sensor. 


C.  Magnetic  Core  Sensors 

Previous  testing  of  our  nanosensors  has  highlighted  a  potential  problem  related  to  the 
presence  of  debris.  In  some  cases,  a  large  amount  (several  grams)  of  powder  was  returned 
that  contained  only  small  amounts  of  sensor  material,  see  Figure  ill.  While  we  were 
usually  able  to  extract  a  fluorescent  signal  from  these  samples,  the  signal  has  been  very 
weak.  In  order  to  enhance  the  S/N,  we  have  developed  an  initial  design  that  contains  a 
magnetic  core,  Fe304.  Having  such  a  magnetic  core,  would  allow  us  to  use  a  magnet  to 
extract  our  sensor  materials  from  the  debris.  Figure  112  shows  an  SEM  of  Fe304/Y203 
core/shell  nanoparticles. 
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Figure  ill.  Debris  from  explosion  tests. 
Only  a  very  small  percentage  of  the  debris  is 
actual  sensor  material. 
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Figure  112.  SEM  of  Fe304/Y203 
core/ shell  nanoparticles. 


The  magnetic  properties  of  Fe304/Eu:Y203  core/shell  nanoparticles  are  demonstrated  in 
Figure  113.  All  samples  are  magnetic,  except  the  one  heated  to  950  °C.  Further  analysis  is 
required  to  determine  the  exact  cause. 


Figure  113.  Images  of  Fe304/Eu:Y203  core/shell  nanoparticles  in  the  presence  of  a 
magnet.  Precursor  (left),  heated  to  600  °C  and  700  °C  (center),  and  heated  to  800  °C 
and  950  °C  (right).  All  samples  are  magnetic,  except  the  one  heated  to  950  °C. 


Figure  114  shows  the  fluorescence  spectra  of  Fe304/Eu:Y203  core/shell  nanoparticles 
heated  to  various  temperatures  for  5  s. 
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Figure  114.  Fluorescence  spectra  of  Fe304/Eu:Y203 
core/shell  nanoparticles  heated  to  various  temperatures 
for  5  s. 


X.  LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


c-Ho:Y203 

crystalline  Ho-doped  yttria 

C02 

Carbon  dioxide 

DTRA 

Defense  Threat  Reduction  Agency 

DSC 

Differential  Scanning  Calorimetry 

Dy 

Dysprosium 

EDS 

Energy  Dispersion  Spectroscopy 

Er 

Erbium 

Eu 

Europium 

Eu:Y203 

Eu-doped  yttria  (1%  unless  otherwise  indicated) 

Eu:Zr02 

Eu-doped  zirconia  (1%  unless  otherwise  indicated) 

FTIR 

Fourier  Transform  Infrared 

FWHM 

Full  Width  Half  Maximum 

Ho 

Holmium 

Ho:Y203 

Ho-doped  yttria 

IH 

Indian  Head 

IHEODTD 

Indian  Head  Explosive  Ordnance  Disposal  Technology  Division 

La203 

Lanthanum  oxide  or  lanthana 

La(OH)3 

Lanthanum  hydroxide 

Ln 

Lanthanide 

M203 

Generic  form  for  metal  oxides  (sesquioxides) 

Nd 

Neodymium 

NSWC 

Naval  Surface  Warfare  Center 

p-Eu:Y203 

precursor  of  Eu-doped  yttria 

p-Eu:Zr02 

precursor  of  Eu-doped  zirconia 

Pr 

Praseodynium 

RE 

Rare  earth 

SEM 

Scanning  Electron  Microscopy 

Sm 

Samarium 

Tb 

Terbium 

TC 

Thermocouple 

TEM 

Transmission  Electron  Microscopy 
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TGA 

Tm 

XRD 

Yb 

Y2(C03)3-nH20 

Y203 

Zr02 

Zr(OH)4 


Thermogravimetric  Analysis 
Thulium 
X-ray  diffraction 
Ytterbium 
Yttrium  carbonate 
Yttrium  oxide  or  yttria 
Zirconium  dioxide  or  zirconia 
Zirconium  hydroxide 
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